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Abstract 

We elucidate the close connection between the repulsive lattice gas in equi- 
librium statistical mechanics and the Lovasz local lemma in probabilistic com- 
binatorics. We show that the conclusion of the Lovasz local lemma holds for 
dependency graph G and probabilities {p x } if and only if the independent-set 
polynomial for G is nonvanishing in the polydisc of radii {p x }- Furthermore, we 
show that the usual proof of the Lovasz local lemma — which provides a suffi- 
cient condition for this to occur — corresponds to a simple inductive argument 
for the nonvanishing of the independent-set polynomial in a polydisc, which 
was discovered implicitly by Shearer ||S| and explicitly by Dobrushin [37J EE] • 
We also present some refinements and extensions of both arguments, including 
a generalization of the Lovasz local lemma that allows for "soft" dependencies. 
In addition, we prove some general properties of the partition function of a 
repulsive lattice gas, most of which are consequences of the alternating-sign 
property for the Mayer coefficients. We conclude with a brief discussion of the 
repulsive lattice gas on countably infinite graphs. 



Key Words: Graph, lattice gas, hard-core interaction, independent- set polynomial, 
polymer expansion, cluster expansion, Mayer expansion, Lovasz local lemma, proba- 
bilistic method. 
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1 Introduction 



The lattice gas with repulsive pair interactions is an important model in equilib- 
rium statistical mechanics pHH EH1 [HZ1 E01 EH El [121 [2IH EDI ED EDI EEII21 E3 
El US]- In the special case of a hard-core self-repulsion and hard-core nearest-neighbor 
exclusion (i.e. no site can be multiply occupied and no pair of adjacent sites can be 
simultaneously occupied), the partition function of the lattice gas coincides with the 
independent- set polynomial in combinatorics (also known as the independence poly- 
nomial or the stable-set polynomial) [H EU E21 EEl EH H21 IS31 EI31 [23 [23 EH EH GUI 
Moreover, the hard-core lattice gas is the universal statistical-mechanical model in the 
sense that any statistical-mechanical model living on a vertex set Vq can be mapped 
onto a gas of nonoverlapping "polymers" on Vq, i.e. a hard-core lattice gas on the 
intersection graph of Vo |ll)Ut Section 5.7]. This construction, which is termed the 
"polymer expansion" or "cluster expansion", is an important tool in mathematical 
statistical mechanics [23 EH El EH 121] , and much effort has been devoted to finding 
complex polydiscs in which the polymer expansion is convergent, i.e. in which the 
(polymer-)lattice-gas partition function is nonvanishing jHHl EH E3 123 EH E01 EH 
HTM E3 EH EH EH EH EH EH EUl CTJ21 HU1- One goal of this paper is to make a 
modest contribution to this line of development. 

The Lovasz local lemma [HH1 EH1 11051 H06| is an important tool in probabilistic 
existence proofs in combinatorics. It provides a lower bound on the probability that 
none of a collection of "bad" events occurs, when those events are subject to a set of 
"local" dependencies, controlled by a "dependency graph". The Lovasz local lemma 
(and its algorithmic versions [HI E]) has found a significant variety of applications, 
especially in graph coloring, Ramsey theory, and related algorithmic problems [21 EH 

E3- 

In this paper we would like to elucidate the close relation between these two 
apparently disparate subjects. Following a seminal (but apparently little-known) 
paper of Shearer jUH], we shall show that the conclusion of the Lovasz local lemma 
holds for dependency graph G and probabilities {p x } if and only if the independent- 
set polynomial for G is nonvanishing in the polydisc of radii {p x }- Moreover, we shall 
show that the usual proof of the Lovasz local lemma — which provides a sufficient 
condition for this to occur — corresponds to a simple inductive argument for the 
nonvanishing of the independent-set polynomial in a polydisc, which was discovered 
implicitly by Shearer [§H] and explicitly by Dobrushin [H3 EH]- Finally, we shall 
present some refinements and extensions of both arguments. 

We have tried hard to make this paper comprehensible to the union (not the 
intersection!) of combinatorialists, probabilists, classical analysts and mathematical 
physicists. We apologize in advance to experts in each of these fields for boring them 
every now and then with overly detailed explanations of elementary facts. 
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1.1 The repulsive lattice gas 

In statistical mechanics, a "grand-canonical gas" is defined by a single-particle 
state space X (here a nonempty finite set), a fugacity vector w = {w x } xe x £ C x , 
and a two-particle Boltzmann factor W: IxI-»C with W(x,y) = W(y,x). The 
(grand) partition function Zw(w) is then defined to be the sum over ways of plac- 
ing n > "particles" on "sites" xi,...,x n G X, with each configuration assigned 
a "Boltzmann weight" given by the product of the corresponding factors w Xi and 
W(xi,Xj): 



= e^ e n 



'1.1a) 



n=0 x\,...,x n £X \i=l / \l<i<j<n 



= E ( n <-^;)"-'^ ) n w{x , y) ^\ (1 , b) 

n \xex x ' / y{x,j/}cx y 

where in fll.lb|) the sum runs over all multi-indices n = {n x } x€ x of nonnegative 
integers, and the product runs over all two-element subsets {x,y} C X (x ^ y). In 
this paper we shall limit attention to the repulsive case in which < W(x, y) < 1 
for all x, y. From this assumption it follows immediately that Ziy(w) is an entire 
analytic function of w satisfying |Zjy(w)| < exp(J2 x ^x \ w x\)- 

If W(x,x) = for all x G X — in statistical mechanics this is called a hard-core 
self-repulsion — then the only nonvanishing terms in (jl.lb|) have n x = or 1 for all x 
(i.e. each site can be occupied by at most one particle), so that Z w (w) can be written 
as a sum over subsets: 




Z w (w) = ^ \[[ w x \\ H W(x,y) | . (1.2) 

X'CX \x£X> J \{x,y}CX> 

In this case Z^(w) is a multiaffine polynomial, i.e. of degree 1 in each w x separately. 
Combinatorially, Zyy(w) is the generating polynomial for induced subgraphs of the 
complete graph, in which each vertex x gets weight w x and each edge xy gets weight 
W(x,y). 

If, in addition to hard-core self-repulsion, we have W(x, y) = or 1 for each pair 
x 7^ y — in statistical mechanics this is called a hard-core pair interaction — then 
we can define a (simple loopless) graph G = (X, E) by setting xy G E whenever 
W(x,y) = and x ^ y, so that Z^(w) is precisely the independent- set polynomial 
for G: 

ZgM = E II w * ( L3 ) 

X' C X x&X 1 
X' independent 

Traditionally the independent-set polynomial is defined as a univariate polynomial 
Zg(w) in which w x is set equal to the same value w for all vertices x jU HU H21 EH EB 
E21 EHl 123 EE EH EZ1 • But one of our main contentions in this paper is that Z G is 
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more naturally understood as a multivariate polynomial; this allows us, in particular, 
to exploit the fact that Zq is multiaffine. 1 

More generally, given any W satisfying < W(x, y) < 1 for all x, y, let us define a 
simple loopless graph G = Gw (the support graph of W) by setting xy G E(G) if and 
only if W(x, y) ^ 1 and x ^ y. 2 The partition function Zy/iyf) can be thought of as a 
"soft" version of the independent-set polynomial for G, in which an edge xy G E(G) 
has "strength" 1 - W(x,y) G (0, 1]. 

In summary, we shall consider in this paper three levels of generality: 

(a) The general repulsive lattice gas in which multiple occupation of a site 
may be permitted. 

(b) The lattice gas with hard-core self- repulsion (jl.2j) . in which multiple occupation 
of a site is forbidden, but in which the interactions between adjacent sites may 
be "soft". 

(c) The lattice gas with hard-core self-repulsion and hard-core pair interaction (jl.3|) . 
which is simply the independent-set polynomial for the graph G. 

Let us now return to the general case of a repulsive lattice gas (II. lj) . Since Zw(w) 
is an entire function of w satisfying Zw{0) = 1, its logarithm is analytic in some 
neighborhood of w = 0, and so can be expanded in a convergent Taylor series: 

logZ^w) = 5> n (W0w n , (1.4) 

n 

where we have used the notation w 11 = n^ex w x x i an d of course Co = 0. In statistical 
mechanics, fll.4|) is called the Mayer expansion |112j , and there is a beautiful formula 
for the coefficients c n (W), whose derivation we will review briefly (Proposition |23J- 
As a corollary of this formula, we will show (Proposition 12.8(1 that the Taylor series 
(jl.4() has alternating signs whenever the lattice gas is repulsive: 

(-lp- l c n (W) > (1.5) 

for < W < 1. (Here |n| = J2xex n v) Moreover, in this interval, |cn(W)| = 
(— l)l n l _1 c n (Vr) is a decreasing function of each W(x, y), i.e. an increasing function 
of the "interaction strength" 1 — W(x,y). 

1 More generally, let G = (X, E) be a graph, and suppose that for x ^ y we have W(x, y) = or 
1 according as xy £ E or xy E; but let < W(x, x) < 1 be arbitrary. Then we have the amusing 
identity Zwiyf) = ^g(w) where 

w * = E n\ W *- 

n—l 

2 Strictly speaking, Gw ought to be called the (simple) support graph of 1 — W, since xy £ E(G) 
if and only if 1 — W(x, y) ^ and x ^ y. In particular, in the case of a hard-core self-repulsion and 
hard-core pair interaction, W is the adjacency matrix of the complementary graph G. Please note 
also that Gw does not contain any information about the diagonal weights W(x, x). 
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As a simple consequence of the alternating-sign property for log Zw, we will then 
prove (Theorem I2.10JI the equivalence of a number of conditions for the nonvanish- 
ing of Z w in a closed polydisc D R = {w: \w x \ < R x for all x}. These equivalent 
conditions will play a central role in our study of the Lovasz local lemma (Section HJ). 

Finally, we will provide (Section E} some sufficient conditions for the nonvan- 
ishing of Zw in a closed polydisc Z)r, based on "local" properties of the interac- 
tion W (or of the graph G). Results of this type have traditionally been proven 
EB1 EH EH 121 EDI E21 EH Ell EH by explicitly bounding the terms in the Mayer ex- 
pansion (jl.4j) : this requires some rather nontrivial combinatorics (for example, Propo- 
sition 1231 below together with the counting of trees). Once this is done, an immediate 
consequence is that Zw is nonvanishing in any polydisc where the series for log Zw is 
convergent. Dobrushin's brilliant idea p??l EIHj was to prove these two results in the 
opposite order. First one proves, by an elementary induction on the cardinality of 
the state space, that Z w is nonvanishing in a suitable polydisc ( Theorem 15. it then 
follows immediately that log Zw is analytic in that polydisc, and hence that its Taylor 
series (jl.4|) is convergent there. In Section El we will prove some refinements of this 
result and investigate their sharpness; and in Sections El and [7| we will provide some 
complementary results that give additional insight into the nature of these bounds. 
Let us remark that the Dobrushin-Shearer inductive method employed in Section El 
is limited, at present, to models with hard-core self-repulsion (jl.2j) . for which Zw is 
a multiaffine polynomial. It is an interesting open question to know whether this 
approach can be made to work without the assumption of hard-core self-repulsion. 

1.2 The Lovasz local lemma 

In combinatorics we are frequently faced with a finite family (A x ) x ex of "bad" 
events in some probability space, having probabilities p x = F(A X ), and we want to 
show that there is a positive probability that none of the events A x occurs. Under 
what conditions can we do this? In general all we can say is that PCQrgx Ax) > 
1 — J2 x Pxi since in the worst case the events A x could be disjoint; so we would need 
^2 x Px < 1 to ensure that Ftflzex ^ x ) > 0. On the other hand, if the events (A x ) xeX 
are independent, then P(f\ex Ac) = TlxO- ~ P x )i which is positive as soon as p x < 1 
for all x. This suggests that if the (A x ) x£ x are in some sense "not too strongly 
dependent" , then it might be possible to prove P(rLex Ax) > ^ under relatively mild 
conditions on the {p x }. This is the situation addressed by the Lovasz local lemma: 
we allow strong dependence among some subsets of the (A x ) xe x, but insist that most 
of these events are independent. Specifically, the local lemma applies to collections of 
events in which the dependencies are controlled by a dependency graph, so that each 
event is independent from the events that are not adjacent to it. 

More precisely, let G be a graph with vertex set X. We say that G is a dependency 
graph for the family (A x ) x ex if, for each x G X, the event A x is independent from 
the cr-algebra generated by the events {A y : y G X \ T*(x)}. [Here we have used the 
notation T*(x) = T(x) U {x}, where T(x) is the set of vertices of G adjacent to x.] 
Erdos and Lovasz proved the following fundamental result: 
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Theorem 1.1 (Lovasz local lemma) Let G be a dependency graph for the family 
of events {A x ) xe x, and suppose that (r x ) xe x are real numbers in [0, 1) such that, for 
each x, 

P(4) < r x H (l-r„) . (1.6) 

y&(x) 

Then P(a e x^) > IL*=x(l ~r m ) > 0. 

Erdos and Spencer (50] (see also [21 ESI) later noted that the same conclusion 
holds even if A x and o~(A y : y G X \ T*(x)) are not independent, provided that the 
"harmful" conditional probabilities are suitably bounded. More precisely: 

Theorem 1.2 (Lopsided Lovasz local lemma) Let (A x ) x€ x be a family of events 
on some probability space, and let G be a graph with vertex set X . Suppose that 
(r x ) x( zx are real numbers in [0, 1) such that, for each x and each Y C X \ T*(x), we 
have 

Then P(a e x^) > n, e x(l " O > 0. 

In fact, the arguments of [HHllini (see also |105[ Hn6j) show that in Theorems ll.il 
and 11.21 a slightly stronger conclusion holds: for all pairs Y, Z of subsets of X we 
have 

F{f]A x \f]A x ) > J] (l-r a ). (1.8) 

x& xGZ x£Y\Z 

In this paper we shall (following Shearer j98 ) approach the problem by dividing 
our discussion into two parts: 

1) A best-possible condition to have P(f\ex ^) > 0, * n t erms of the independent- 
set polynomial Z G {— p); and 

2) A sufficient condition to have Z G (— p) > 0, along the lines of Lovasz, Dobrushin 
and Shearer. 

We shall treat the first problem in Section and the second problem in Sections EH3 
We shall also prove a generalization of the Lovasz local lemma that allows for "soft" 
dependencies (Theorems 14.21 and 15. 4j) . 

Let us remark that we have been able to relate the Lovasz local lemma to a 
combinatorial polynomial (namely, the independent-set polynomial) only in the case 
of an undirected dependency graph G. Although the local lemma can be formulated 
quite naturally for a dependency digraph 123 EHj; we do not know whether the 
digraph Lovasz problem can be related to any combinatorial polynomial. (Clearly 
the independent- set polynomial cannot be the right object in the digraph context, 
since exclusion of simultaneous occupation is manifestly a symmetric condition.) 
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1.3 Plan of this paper 

The plan of this paper is as follows: In Section |2 we prove some general properties 
of the partition function of a repulsive lattice gas. In Section|3]we study the additional 
properties that arise in the case of a hard-core self- repulsion (jl.2|) . for which Zw is 
a multiaffine polynomial. We also calculate some simple examples; in particular, we 
show that when the support graph G = Gw is a tree, the partition function Z w 
can easily be calculated by working upwards from the leaves. In Section 0] we prove 
our lower bound on ^Hzex A x ) > in terms of the independent-set polynomial 
Zq(— p). We also prove a "soft-core" generalization of this result. In Section |3] 
we give some sufficient conditions for the nonvanishing of Zw in a polydisc, and 
investigate their sharpness. In Section |H1 we show how the bounds of Section |3] can 
be interpreted in terms of a lattice gas on either the "self-avoiding-walk tree" or the 
"pruned self-avoiding-walk tree" of G. In Section Owe show how this tree bound can 
be understood as arising from the repeated application of a single "unfolding" step. 
Finally, in Section |H] we study the repulsive lattice gas on an infinite graph. 

The reader primarily interested in the Lovasz local lemma can read the statement 
of Theorem I2.1UI (skipping the proof) and Definition 12.141 quickly read Section 13.11 
and then jump directly to Section 01 The reader primarily interested in the lattice 
gas can skip Section HJ 

2 The lattice-gas partition function 

In this section we prove some general properties of the partition function of a 
repulsive lattice gas. We begin with some elementary identities (Section 12. In 
Section l2~2l we derive the Mayer expansion and prove some of its properties, notably 
the alternating-sign property (Proposition 12. 8j) . In Section |2~31 we state and prove 
the "fundamental theorem" (Theorem I2.1()j) . which sets forth a number of equivalent 
conditions defining the set that we shall call TZ(W) [Definition I2.14j and that plays 
a central role in the remainder of this paper. In Section 12.41 we derive the principal 
properties of the set 1Z(W). In Section 12*31 we derive some further consequences of 
the alternating-sign property. In Section 12.61 (which is a digression from the main 
thread of the paper and can be omitted on a first reading) we discuss the algebraic 
irreducibility of the multivariate partition function Zy/{w). Finally, in Section |2~7I we 
make a brief remark concerning the convexity of logZvi/(w) at nonnegative fugacity 
w. 

Nearly all of the results in this section are valid for an arbitrary repulsive lattice 
gas (jl.lj) . in which multiple occupation of a site is permitted. A few of the results are 
restricted to the case of a hard-core self- repulsion (jl.2|) . in which multiple occupation 
of a site is forbidden. 
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2.1 Elementary identities 

Consider a general lattice gas (jl.l)) defined on a finite set X. For any subset 
A C X, let us define the restricted partition function 

Zw,a( w ) = Z w (wl A ) (2.1) 

where 

(wl A ), = l Wx if f GA (2.2) 
1 otherwise 

This simply forces the sites in X \ A to be unoccupied. Since W will be fixed through- 
out this subsection, we shall often omit it from the notation and write simply Z A . 

For any x G A, we can expand the partition function (jl.lb|) in terms of the 
occupation number n x : we obtain the fundamental identity 



Z A (w) = ^^W(x,x)^^ a Z ASjB (W(x,') n w) (2.3) 

n=0 

where 

[W(x,-) n w] y = W(x,y) n w y . (2.4) 

In the special case of a hard-core self-repulsion at site x [i.e. W(x,x) = 0], only 
the terms n — 0,1 appear in this sum; this case will play a key role in this paper 
starting in Section 13.11 (Analogous identities can be derived where we expand in 
the occupation numbers at an arbitrary set of sites S C A; for example, a two-site 
identity for the hard-core case will play a major role in Section 0) 
From ()2.3|) we can deduce the differentiation identity 



OZs(w) ir"- 1 

n=l 



dw x 



E r^TTT W ^ x) n{n ~ 1)/2 Za\ x (W(x, w) (2.5a) 
z — ' In — 1)! 

?1=1 v 7 

00 n 

= y"^W(x,x) n{n+1)/2 Z A \ x (W(x,-) n+1 w) (2.5b) 
z — ' n! 

n=0 

= E f™' 1 " Wfo x)"^ 1 )/ 2 Za^CWCx, -)"W(x, -)w) (2.5c) 

n=0 

= Z A (W(x,». (2.5d) 
Repeated application of ()2.5|) yields the multiple differentiation identity 
d n Z A (w) 



dw xi ■ ■ ■ dw Xn 



11 W(xi, x 3 ) j Z A (f[ W(x h •) w\ 



(2.6) 



where (n W^,, ■) w) x = w x Y[ W(x h x) 
t=i i=i 
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2.2 The Mayer expansion 

We begin by reviewing the derivation of the Mayer expansion ()1.4j) . The first step 
is to trivially rewrite the partition function (jl.laj) as 

oo / n \ 

z ^) = H- [ E IK E W F( x ii x j) > (2.7) 

n=0 ' n,..,j!„eX \i=l / GeGn ijeE(G) 

where is the set of all (simple loopless undirected) graphs on the vertex set 
{1, . . . , n}, and 

F(x,y) = W(x,y)-1 (2.8) 
is called the two-particle Mayer factor. This is of the form 

n=0 ' GeGn 

where the weights 

= e (ik) n p ^ x i) ( 2 - iQ ) 

xi,...,x n eX \i=l J ijeE(G) 

satisfy 

(a) W(0) = 1; 

(b) W(G) = W(G") whenever G = G" (i.e. whenever G and G" differ only by a 
relabelling of vertices); and 

(c) W{G) = W(G'i) W^Gy whenever G is isomorphic to the disjoint union of G\ 
and G2. 

It then follows from the exponential formula |112[ 11151 E3 H08j that 

00 1 

iog^(w) = E^E W ( G ) ( 2 - lla ) 

n=0 ' GeC„ 

00 / n \ 

= E^ E IK E n (2.nb) 

n=0 n,..,x„eJC \i=l / GeC„ijeE(G) 

at least in the sense of formal power series in w, where C n C Q n is the set of connected 
graphs on {1, . . . , n}. Therefore: 
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Proposition 2.1 The coefficients c n (W) of the Mayer expansion \l-4\ i are given by 

Cn (w) = i yi En F ^ x i) ( 2 - i2a ) 

n,...,i B ex Gec n ij£E(G) 

#{i: Xi = x} = n x \/x 

E E II F(x,y) eG{Sx ' Sy) ]jF(x,x) eG ^ (2.12b) 

(S x ) xex GeC n {x,y}CX xGX 

= {l,...,n} 

[5x1 = tij Vie 

where n = \n\, t/ie /irst sum zn \2.12b\) runs over all partitions of {1, . . . ,n} into 
disjoint subsets (S x ) xe x with the specified cardinalities, and ec(S x ,S y ) [resp. ec{S x )] 
denotes the number of edges of G connecting S x to S y [resp. within S x ]. 

Remark. Even in the simple case X = {x} and W(x,x) = 0, for which logZw(w) = 
log(l + w) = Y^nLi — wn i the Mayer expansion (|2.11|) implies the nontrivial 
identity 

J2(-1)\ E ^\ = (-l) n_1 (n - 1)! (2.13) 

Gec„ 

for the generating function of connected spanning subgraphs of the complete graph 
K n . We leave it as an exercise for the reader to find a direct combinatorial proof of 

In order to analyze the Mayer coefficients c n (W), it is convenient to proceed in 
a bit more generality. So let H = (V, E) be a graph — we allow loops and multiple 
edges — and let z = {z e } e£ E be a family of complex edge weights for H. [Later we 
will specialize to H = K n .] Define the generating function of connected spanning 
subgraphs of H (or "connected sum" for short), 

c h(*) = e n * ( 2 - 14 ) 

E'CE e€E' 
(V, E') connected 

It is easy to see that Ch satisfies the deletion-contraction relation 

C H (z) = C H \e{l^e) + Z e C H /e{z?e) (2-15) 

for any edge e £ E: here H \ e is the graph H with edge e deleted, H/e is the graph 
H with edge e contracted (note that we do not delete any loops or multiple edges 
that may be formed), and z^ e denotes the family {^/}/g£;\e- Please note that if e is 
a loop, then H/e — H \ e and hence 

C H (z) = C H \ e {^ e ) + Z e C H/e {-L^ e ) = (1 + Z e ) C H \e(z^ e ) . (2.16) 

Let C (resp. T) be the set of subsets E' C E such that (V, E') is connected (resp. 
is a tree). Clearly C is an increasing family of subsets of E with respect to set- 
theoretic inclusion, and the minimal elements of C are precisely those of T (i.e. the 
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spanning trees). It is a nontrivial but well-known fact fTf\ Sections 7.2 and 7.3] |117| 
Section 8.3] that the (anti-) complex C is partitionable: that is, there exists a map 
R: T — > C such that R(T) D T for all T G T and C = i+J Te7 -[T, R(T)] (disjoint union), 
where [Ei,E 2 ] denotes the Boolean interval {E'\ E\ C E' C E^}- m fact, many 
alternative choices of R are available [T71 Sections 7.2 and 7.3] Sections 2 and 6] 
|15| Proposition 13.7 et seq.] |101| Proposition 4.1], and most of our arguments will 
not depend on any specific choice of R. However, at one point (Proposition 12.6(1 we 
shall need the existence of an R with the following additional property: 

Lemma 2.2 Let H = (V,E) be a connected graph, and fix any Tq G T. Then there 
exists a map R: T — > C such that 

(a) R(T) D T for all T G T; 

f&j C is the disjoint union of the Boolean intervals [T, R(T)], T G T; and 
fcj R(T )=T . 

Proof. If H has one vertex and no edges, then T = C = {0} and the result holds 
trivially; so assume henceforth that E ^ 0. Assign arbitrary weights w e > chosen 
so that no two spanning trees have equal weight (for example, one can choose the 
w e to be linearly independent over the rationals). For each E' G C, let S(E') be the 
(unique) minimum-weight spanning tree contained in E'. (This can be constructed 
by a greedy algorithm, i.e. start from and keep adding the lowest-weight edge in 
E' that does not create a cycle. See, for instance, [22 Section 1.2].) We then define 
R(T) to be the union of all E' that have S(E') = T. To verify that this works, we 
need to show that if S(Ei) = S(E 2 ) = T, then S(Ei U E 2 ) = T; but this follows easily 
from the validity of the greedy algorithm. [Note that this construction includes, as a 
special case, the lexicographically minimum spanning tree for any ordering of E: it 
suffices to take w n = 2 n .] 

By choosing the weights so that w e > Wf whenever e G T and / G E \ T , we can 
ensure that for E' ^ Tq the first edge chosen by the greedy algorithm will not belong 
to T . Therefore, R(T ) = T . | 

Remark. The greedy algorithm works, more generally, for an arbitrary matroid [HE1 
Section 1.8]. It follows that the independent-set complex of a matroid (or, dually, the 
spanning-set anti-complex) is partitionable f7| Sections 7.2 and 7.3]. 

Given the existence of R, we have the following simple but fundamental identity: 

Proposition 2.3 (partitionability identity) Let R: T — > C be any map satisfying 
R(T) D T for allT G T andC is the disjoint union of the Boolean intervals [T, R(T)], 
for T G T '. Then 

c H (z) = £ n*« e n * 

TCE eGT TCE'CR(T) e£E'\T 
(V, T) tree 
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- e n*» n ( 2 - 17 ) 

TCE eeT eeR(T)\T 
(V, T) tree 

This identity (for one specific choice of R) is due originally to Penrose |88j . 

One immediate consequence of the identity (J2.17j) is the following inequality valid 
in the "complex repulsive" regime 1 1 + z e | < 1 : 

Proposition 2.4 (Penrose |88j) Let H = (V,E) be a finite undirected graph equipped 
with complex edge weights {z e } e( zE satisfying |1 + z e \ < 1 for all e. Then 

1^)1 < e 

tce eer 

(V,T) tree 

This bound plays a major role in several traditional (pre-Dobrushin) proofs of con- 
vergence of the Mayer expansion jHHl HUH E3 12H H0U| EH H04j , as well as in a recent 
bound on the zeros of chromatic polynomials |lUlj . 
Let us now define the "generalized connected sum" 

c H (\-z) = Yi xc(E,) ^ 2 - i9a ) 

E' C E ee-E' 
(V, E 1 ) connected 

= \-^-^C H (\z) (2.19b) 

where c(E') = \E'\ — \V\ + 1 is the cyclomatic number of the connected subgraph 
(V, E') [more generally, for a subgraph (V, E') with k(E') components, the cyclomatic 
number is c(E') = \E'\ - \V\+k(E')). Of course, ()2.19bjl shows that C H (\] z) contains 
no more information than C^(z); it is just a convenient way of scaling all the variables 
z e simultaneously. The function C#(A; z) interpolates between the tree sum (A = 0) 
and the connected sum (A = 1); and Proposition 12.41 can be rephrased as saying that 
|C#(l;z)| < Ch(0; |z|). One of us has conjectured that the absolute-value signs can 
in fact be put outside the sum, i.e. |C#(1; z)| < |C#(0; z)| jlUH Remark 2 in Section 
4.1]. In general this is still an open problem; but if the z e are real and lie in the 
interval [—1, 0] - which corresponds to the physical repulsive regime < W < 1 - 
then a vastly stronger result is true: 

Proposition 2.5 (Sokal [101J) Let H = (V,E) be a finite undirected graph equipped 
with real edge weights {z e } eG E satisfying — 1 < z e < for all e. Then 

(-l)^y\-^C H (X;z) > (2.20) 
on < A < 1, for all integers k > 0. In particular, setting k = and A = 1 we have 

(-l)M-itf H (z) > . (2.21) 



(2.18) 
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Proof. By (l2~TTI) and (|2~T9EJ) . we have 



C H (X;z) = Yl H z - II (! + A ^) ( 2 - 22 ) 

T C E eST e£R(T)\T 
(V, T) tree 

and hence 

^C H (\;z) = k\ E II* II + (2-23) 

TCE T C R(T) \ T eeTUT e€R(T)\(TUT) 
(V, T) tree |f | = k 

which has the claimed sign whenever < A < 1 and — 1 < z e < for all e. 3 ■ 

We shall also need a variant of Proposition 12.51 with strict inequality for the case 
k = 0: 

Proposition 2.6 Let H = (V, E) be a finite undirected graph equipped with real edge 
weights {z e } eG E satisfying — 1 < z e < for all e. Assume furthermore that the 
subgraph (V, suppz) is connected [here suppz = {e E E: z e ^ 0}]. Then 

{-l)\ y \- l C H {z) > 0. (2.24) 



Proof. Choose T C suppz such that (V, To) is a spanning tree, and use Lemma 
to choose R so that R(Tb) = To. Then, in the sum (12.22)1 . the term T = To is nonzero 
with sign (— l)' y ' _1 , and all the other terms either have the sign (— l)'^'" 1 or else are 
zero. ■ 

We can also prove some inequalities on the partial derivatives of Ch{z) with 
respect to individual weights z e . Note first that the deletion-contraction identity 
(I2.15J) implies that 

dC H (z) 

—q- z — = C H/e (z^ e ) . (2.25) 

The graph H/e has \V\ vertices if e is a loop, and \V\ — 1 vertices if e is not a loop. 
Repeated application of these facts together with (|2.21|) yields the following result: 

Proposition 2.7 Let H = (V, E) be a finite undirected graph equipped with real edge 
weights {z e } e( zE satisfying — 1 < z e < for all e. Let k > 1 and e±, . . . , e k G E. 

(a) If ex, . . . , efc are not all distinct, we have 



dz ei ■ ■ ■ dz ek 



(2.26) 



3 Equation (4.7) of |1(J1| inadvertently omitted the prefactor kl in (I2.23|l . 
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(b) If e\, . . . , Ck are all distinct and form a subgraph with cyclomatic number c, we 
have 

(-l)[r[-*+c-i 9Ch{z } > o. (2.27) 
dz ei ■ ■ ■ dz ek 

Let us now specialize these results to the Mayer expansion (|2.12a|) by taking 
H = K n (where n — |n| and K n denotes the complete graph on n vertices), Zij = 
F(xi, Xj) = W(xi, Xj) — 1 and then summing over x\, . . . , x n G X with the specified 
cardinalities: we get 



c n {W) = i C kM*)) (2-28) 



xi, . . . ,x„ € X 

Xi = x} = n x Vx 

where x = (xi, . . . , x n ) and z(x)y = F(x{, Xj). Since any subgraph of K n with 0, 1 or 
2 edges has cyclomatic number 0, we can use Propositions 12.51 and 12.71 to deduce the 
following: 

Proposition 2.8 (signs of Mayer coefficients) Suppose that the lattice gas is re- 
pulsive, i.e. < W(x, y) < 1 for all x,y G X. Then, for all x,y,x',y' G X, the 
Mayer coefficients c n (W) satisfy 

(-l) |nM c n (W0 > (2.29) 
(-l)H-i < o (2.30) 

/ lV n|-l 9 2 C n (W) > 

v ; dW(x,y)dW{x',y') ~ K ' ' 

PROOF. (j2.29j) is an immediate consequence of (|2.28j) and (|2.2ip . For (I2.3U|) . note 
first that for any fixed X\, . . . , x n G X, we have 



E 



dW(x,y) ^ dza 

V ' {i,j}eE(x,y;xi,...,x n ) 3 



(2.32) 

z=z(x) 



where E(x, y;xi,..., x n ) is the set of unordered pairs {i, j} {i ^ j) such that {x^ Xj} = 
{x, y}. Since any subgraph of K n with one edge has cyclomatic number 0, the inequal- 
ity ()2.3()j) now follows from ()2.27|) . Similarly, the left-hand side of (12.31)1 gives rise to 
a double sum over pairs {i, j} G E(x, y;xi, . . . , x n ) and {i 1 , j'} G E(x', y'\ x%, . . . , x n ). 
The diagonal terms (if any) in this sum vanish by ()2.26|) . and the other terms have 
the claimed sign by f)2.27j) . I 

Please note that ()2.31|) can be nonzero even when x = x' and y = y', because 
distinct edges of K n could correspond to the same pair x, y. Note also that the sign 
of the third and higher derivatives cannot be controlled by this method, because a 
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subgraph of K n (n > 4) containing three or more edges could have cyclomatic number 
of either parity. Indeed, for the simple case X = {x} and W(x,x) = W we have 

~ w n W n(n-l)/2 

Zw(w) = (2.33) 

n! 

n=0 

and hence 

, ry , S F 2 3F 2 +F 3 3 

log Z w {w) = w + —w H w 

2 6 

+ — w 4 + 0{w 5 ) (2.34a) 

W-l 2 (W + 2)(W -I) 2 o 

= w + — - — + - ^ -w 3 

2 6 

+ + + + n -Lw 4 + 0(w 5 ) , (2.34b) 

has no fixed sign on < W < 1. 

We can also determine the precise conditions under which the inequality (|2.29|) is 
strict: 



so that 



Proposition 2.9 (condition for Mayer coefficients to be nonzero) Suppose that 
the lattice gas is repulsive, i.e. < W(x,y) < 1 for all x, y £ X. Let n = (n x ) xe x be 
a multi-index with support suppn = {x G X: n x ^ 0}. Then (— l)l n l _1 c n (W) > if 
and only if one of the following is true: 

(a) |n| = 1; 

(b) |n| > 2, supp n = {x} and W(x, x) ^ 1; or 

(c) | suppn| > 2 and the induced subgraph GV [suppn] i s connected. 



Proof. In case (a), c n (W) = 1 > 0; in case (b) or (c), for each term in ()2.28|) [actually 
they are all the same!] the subgraph of K n with edge set suppz(x) is connected, so 
that (-l)l n l- 1 C A -„(z(x)) > by Proposition E3 
There are only two other possibilities: 

(d) n = 0; 

(e) | suppn | > 2 and GV [suppn] is disconnected. 
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In case (d), clearly c n (W) = 0; and in case (e), for each term in ()2.28|) the subgraph 
of K n with edge set suppz(x) is disconnected, so that (— l)l n l-! Cjf n (z(x)) = 0. I 

Historical remark. The alternating-sign property ()2.29|) for the Mayer coefficients 
of a repulsive gas has been known in the physics literature for over 40 years: see 
Groeneveld [48] for a brief sketch of the proof, which uses methods quite different 
from ours. 4 Nevertheless, this result does not seem to be as well known as it should 
be. To the best of our knowledge, the inequalities ()2.30j) and ()2.31|) are new. We think 
that the Mayer coefficients c n (W), and more generally the "connected sums" C#(z), 
merit further study from a combinatorial point of view; we would not be surprised if 
new identities or inequalities were waiting to be discovered. 

2.3 The fundamental theorem 

Let us now state the principal result of this section: 

Theorem 2.10 (The fundamental theorem) Consider any repulsive lattice gas, 
and let R = {R x } x( zx > 0. Then the following are equivalent: 

(a) There exists a connected set C C (— oo,0] x that contains both and — R, such 
that Z^(w) > for all w G C . [Equivalently, — R belongs to the connected 
component of Z^(0, oo) R (— oo, 0} x containing 0.] 

(b) Zyy(w) > for all w satisfying — R < w < 0. 

(c) Z w (w) ^ for all w satisfying |w| < R. 5 

(d) The Taylor series for \ogZy/{w) around is convergent at w = — R. 

(e) The Taylor series for log Z\y(w) around is absolutely convergent for |w| < R. 

Moreover, when these conditions hold, we have \Z\y(w)\ > Zyj{— R) > for all w 
satisfying |w| < R. 

In the case of hard-core self-repulsion, (a)-(e) are also equivalent to 

(V) Z w (-R l s ) > for all SCI, where 

(Rls)x = {?* if * eS (2.36) 
l. otherwise 

4 Groeneveld 0Hj writes that "a more detailed account of this work, containing also the corre- 
sponding results for the pressure and the distribution functions will be published in Physica" , but 
to our knowledge that more detailed paper never appeared. See also Penrose |F7| eq. (8.1)] for a 
related result. 

5 Here we use the notation |w| = {Iw^H^gx. This conflicts slightly with our notation |n| = 
J2xgx nx f° r multi-indices, but we trust that it will not lead to any confusion. 
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(f) Z w (-R) > 0, and (-l)\ s \Z w (-R; S) > for all S C X, wh 



ere 



Z w (w;S) = Yl [H w *)\ II W ^y)\ ■ ( 2 - 37 ) 

SCX'CX \xeX' J \{x,y}CX' J 

(g) There exists a probability measure P on 2 X satisfying P{0) > and 



E p ( T ) = n w(x,y)\ (2.38) 

TDS \x€S / \{x,j/}C5 / 

for all S C X. [This probability measure is unique and is given by P{S) = 
(-1)^Z W (-R;S). In particular, P(0) = Z w (-R) > O.J 

Remarks. 1. The conditions (b'), (f) and (g) are inspired in part by Shearer P%1 
Theorem 1]. 

2. Suppose that the univariate entire function Zw{w), defined by setting w x = w 
for all x, is strictly positive whenever —R < w < 0. Then in fact Zyy(w) > 
whenever —R < w x < for all x: this follows from (a) =^> (b) by taking C to be the 
segment [-R, 0] of the diagonal. 

The proof of Theorem I2.1UI will hinge on the alternating-sign property (J2.29)) for 
the Taylor coefficients of log Zw- In preparation for this proof, let us recall the 
Vivanti-Pringsheim theorem in the theory of analytic functions of a single complex 
variable Theorem 5.7.1]: if a power series f(z) = X]^=o a ™ 2?1 with nonnegative 
coefficients has a finite nonzero radius of convergence, then the point of the circle 
of convergence lying on the positive real axis is a singular point of the function 
/. Otherwise put, if / is a function whose Taylor series at has all nonnegative 
coefficients and / is analytic on some complex neighborhood of the real interval 
[0,-R), then / is in fact analytic on the open disc of radius R centered at the origin 
and its Taylor series is absolutely convergent there. Here we will need the following 
multidimensional generalization of the Vivanti-Pringsheim theorem: 

Proposition 2.11 (multidimensional Vivanti-Pringsheim theorem) LetC be 
a connected subset of [0, oo) n containing ; let U be an open neighborhood of C in 
C n , and let f be a function analytic on U whose Taylor series around has all non- 
negative coefficients. Then the Taylor series of f around converges absolutely on 

the set hull(C) = [J Z)r, where Z)r denotes the closed polydisc {w 6 C n : \wi\ < 
Rec 

Ri for all i), and it defines a function that is continuous on hull(C) and analytic on 
its interior. 



Remark. We could equally well start from an open neighborhood U of in C n , and 
define C to be the connected component of U H [0, oo) n containing 0. This is the 
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maximal set C compatible with the given U; note that it is open in [0, oo) n , and that 
hull(C) is then open in C n . 

Proof. Let /(z) = Y2 n anZn ^ e ^ ne Taylor series of / around 0. And let us define 

S = {R G [0, oo) n : X] n a n z n converges absolutely at z = R} (2.39a) 
= {R G [0, oo) n : J2n a n zn converges absolutely for z G D n } (2.39b) 

Note that S is a down-set (that is, R G S and < R' < R imply R' G S) and that 
G S] we shall show that C C S. 

Consider any point z G 5 fl C (here 5 denotes the closure of S). Choose e > 
such that the closed polydisc of radius e in each direction around z — call it D(z , e) 
- is contained in U. Then choose zi G S fl D(z ,e/5), and choose z 2 G D(zi,e/5) 
such that < Z2 < zi and (z 2 )j < (zi)j for all coordinates i with (zi)j > 0. It follows 
that _D(z 2 ,3e/5) C U. Now, since zi G S and z 2 is strictly below zi in all nonzero 
coordinates, the Taylor series around for / and all its derivatives converge absolutely 
at z 2 . And the Taylor series for / around z 2 converges absolutely in D(z 2 ,3e/5). If 
Z3 G D(z2, 3e/5) with Z3 — z 2 > 0, we can write 

A»s) = E^l^( z 3- Z 2) k (2.40a) 

k 

k n>k ^ ' 

= E Q(z 3 -Z 2 ) k zr k (2.40C) 

n 0<k<n ^ ' 

= E fl nZ3» (2-40d) 
n 

where the rearrangements are justified because all the terms are nonnegative, so that 
the convergence of the iterated sum (j2.40b|) implies the absolute convergence of the 
corresponding double sum. It follows that Z3 G S. In particular we can choose 
z 3 = z 2 + (3e/5, . . . , 3e/5). Since S is a down-set, and since the the real points of the 
polydisc -D(z , e/5) all lie below z 3 , we conclude that D(z , e/5) fl [0, 00)" c S. Hence 
D(z , e/5) fl C C S D C. Since z is an arbitrary point of 5 fl C, it follows that S C\C 
is both closed and open in C (and nonempty since G S fl C). Since C is connected, 
we conclude that S fl C = C. I 

For the proof of Theorem 12.101 we need an elementary topological lemma: 

Lemma 2.12 Let K be a convex set in M. n , let C be an open connected subset of K (in 
the relative topology), let U be an open neighborhood of C in M. n , and let Xi,x 2 G C. 
Then there exists a finite polygonal path P C C running from x\ to x 2 and a simply 
connected open set U' in W 1 satisfying P C U' C U . 
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PROOF. It is a well-known fact that any two points in C can be connected by a finite 
polygonal path P lying in C. [Sketch of proof: Define an equivalence relation ~ on 
C by setting x ~ y iff there exists a finite polygonal path in C connecting x to y. 
Because K is convex and C is open in K, the equivalence classes of ~ are (relatively) 
open subsets of C. Since C is connected, there is just one equivalence class.] By 
removing loops, we can assume that P is non-self-intersecting. 

Then we can take U' to be a sufficiently small tube centered on P. [Let 5 be 
the minimal distance between non-adjacent segments of P, and let 5' be the distance 
from P to the complement of U; we have 5' > by compactness of P. Then let U' 
be the set of all points whose distance from P is less than min(5/3, S'/2). It is not 
hard to see that U' is simply connected. 6 ] ■ 

We shall also make use of the following elementary result: 

Lemma 2.13 Let F be a function on 2 X , and define 



Then 



F -( s ) = E F ( X ') 

x'cs 

X'DS 

F -(s) = £(-i) |y| F + CO 

YCS C 



(2.41) 
(2.42) 

(2.43) 



where S c = X \ S . 

Proof. This is a straightforward application of inclusion-exclusion (see e.g. [107, 
Section 2.1] or Chapter IV]). We have 

F(X') = (-l) inX V + (F) . (2.44) 

Hence 

F -(S) = E = EE (-l) inX V + (F) (2.45a) 

X'CS X'CSY^X' 

= E F +( y ) E (- 1 ) inx ' 1 ( 2 - 45b ) 

Y X'CYnS 

6 It suffices to use repeatedly the following fact: If U\ and Ui are simply connected open subsets 
of R ra , and U\C\Ui is nonempty and connected, then Xl\ U U2 is simply connected. For a proof of this 
lemma, see e.g. [Ill pp. 46-47]. This is a special case of the Seifert-van Kampen theorem, which 
gives a recipe for computing the fundamental group (first homotopy group) of U\ U U2 in terms of 
the fundamental groups of U\, U2 and U\ PI U2 [7^1 Chapter 4]. 
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^F + (F)(-l)l ynSC l (-l) l(ynSAX '' (2.45c) 

Y X'CYnS 

J2F + (Y)(-l) lYnSCl 5 Y ns,0 (2.45d) 

Y 

(-l) |y| ^+On • (2.45e) 



yc5 c 



Proof of Theorem 12. 101 (c) (b) (a) is trivial. 

(e) ==>- (d) is trivial, while (d) ==>- (e) follows from the alternating-sign property 
(Ell). 

(e) implies that the sum of the Taylor series for logZv^(w) defines an analytic 
function on the open polydisc D R and a continuous function on the closed polydisc 
Dr. Its exponential equals Z w (w) on _D R and hence by continuity also on D R . 
Therefore (e) =>■ (c). 

Finally, assume (a). Since Z w is continuous on C x (and has real coefficients), we 
can find an open connected neighborhood C of C in (— oo, 0] x on which Zw > 0, 
and an open neighborhood U of C in C x ~ R 2 I X on which Zw ^ 0. Applying 
Lemma f2. 121 (with n = 2|X|), we can find a finite polygonal path P C C from to 
— R and a simply connected open set U' in C x satisfying P C U' C U. Then log Z^y 
is a well-defined single- valued analytic function on U', once we specify log Zw(0) = 0. 
Applying Proposition 12 . 1 II to log Zw on P and U' [using the alternating-sign property 
([2.29)1 ] . we conclude that the Taylor series for log Z w around is absolutely convergent 
on Z) R . Therefore (a) =^> (e). 

The bound |Zw(w)| > Z w {— R) for |w| < R, which is equivalent to Re log Z w {w) > 
\ogZw{— R), is an immediate consequence of the alternating-sign property ()2.29|) . 

Now consider the special case of a hard-core self-repulsion, (b) (b') is trivial, 
and (b') =>■ (b) follows from the fact that Zw is multiafnne (i.e. of degree < 1 in each 
w x separately) because the value of Z w at any point w in the rectangle — R < w < 
is a convex combination of the values at the vertices. 

To show that (b) =^> (f), note that 

Z w (w;S) = (n>-) II W(x,y)\z w (W(S,-)w) (2.46) 
where we have defined 

[W(S,-)w] y = (]JW(x,y))w y (2.47) 

\xeS J 

(note in particular that this vanishes whenever y e S). Hence 
(-lf\Z w (-R;S) = [HrJ I J] W(x,yUZ w (-W(S,-)R) > (2.48) 

\x£S J \{x,y}CS J 
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since — R < —W(S,-)R < 0, with strict inequality when \S\ = or 1 [since the 
product over W(x, y) is in that case empty]. 

To show that (f) =>- (b'), use Lemma 12. 131 applied to the set function 

f(s) = \T[-R*) ( n w &y)) ■ ( 2 - 49 ) 

XxeS J \{x,y}CS J 

We have 

F.(S) = Zw(-Kl s ) (2.50) 

F+(S) = Z W (-R;S) (2.51) 

so that Lemma f2 . 1 31 asserts the identity 

Z W (-R1 S ) = ^(-l)l y l^(-R ; r) . (2.52) 

YCS C 

By (f), the Y = term is > and the other terms are > 0, so R I5) > for 

all S. 

Finally, let us show that (f) <^=^ (g). By inclusion-exclusion, there are unique 
numbers P(T) satisfying (|2.38p . namely P(T) = (— l)' T 'Ziy(— R; T). Moreover, tak- 
ing S = in (J2.38J1 we see that ^j>P{T) = 1. Therefore, P is a probability mea- 
sure if and only if (— l)' T 'Ziy(— R; T) > for all T; and P(0) > if and only if 
Z w (-R; 0) = Z w (-R) > 0. ■ 



2.4 Properties of the set 1Z(W) 

The following definition plays a central role in the remainder of this paper: 

Definition 2.14 (definition of 1Z(W)) We define 1Z(W) to be the set of all vectors 
R > satisfying the equivalent conditions (a)-(e) of Theorem \2.1(A When W is the 
hard-core pair interaction for a graph G, 

xrr / \ JO if x = y or xy E E(G) (r> 
W[X ' V) ~ \l if x ^ y and xy i E{G) 

we also write 7Z(G). 

Proposition 2.15 (elementary properties of TZ(W)) For any repulsive lattice gas, 
the set K(W) is 

(a) open in [0, oo) x 

(b) a down-set [i.e. R e TZ{W) and < R < R imply R G TZ(W)J 
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(c) logarithmically convex [i.e. R,R' G K(W) and < A < 1 imply R A (R') 1_A G 

n(wy. 

Moreover, the set 1Z{W) is an increasing function of each W(x,y) on < W < 1. 

Proof. It is obvious from Theorem 12. lUf b) or (c) that TZ(W) is open in [0, oo) x and 
is a down- set. 

Given any formal power series /(z) = X] n a n zn with complex coefficients, let us 
define its set of radii of absolute convergence by 

cvg(/) = {R > 0: l«n|R n < oo} . (2.54) 

n 

It then follows immediately from Holder's inequality for infinite series that cvg(/) 
is logarithmically convex. 7 By Theorem I2.10f e) we have TZ(W) = cvg(log Zw), so 
7Z(W) is logarithmically convex. 

Finally, from (|2.29J) / (|2.30jl we see that |c n (W)| is a decreasing function of each 
W(x,y) on < W < 1, so TZ(W) = cvg(logZvy) is an increasing function of each 
W{x,y) on < W < 1. ■ 

See Corollary 12.191 below for a strict monotonicity that strengthens Proposi- 
tion ^. 15T b): and see Corollary 12 . 291 for an additional convexity property in the special 
case of a hard-core self-repulsion. 

Remark. If A is any nonempty subset of X, we can obviously define a lattice gas on 
A with interaction W \ A, and there will be a corresponding set 1Z(W \ A) C [0, oo) A . 
Now, for any w G C A , we trivially have Zyy^(w) = Z\y(w, 0) where e C X ^ A ; hence 
R G K(W \ A) if and only if (R, 0) G K{W). So the sets K(W \ A) are the sections 
through of K(W). 

Next let us show that the sets 1Z(W) are bounded in a suitable sense. In the 
case of hard-core self-repulsion this is immediate, because 1Z{W) C [0, l) x . [This 
containment is trivial if \X\ = 1, and otherwise follows from Proposition I2.15r b). 
Alternatively, it is trivial if W(x, y) = 1 for all x ^ y, and otherwise follows from 
the monotonicity statement in the last sentence of Proposition 12.151 ] However, in 
the general case the set 1Z(W) can be unbounded, as is shown by the following two 
examples: 

(a) Let X = {x} and W(x, x) = 1; then Z w (w) = e w and K{W) = [0, oo). 

(b) Let X = {x,y}, W(x,x) = W(y,y) = 1 and W(x, y) = 0; then Zw(w x ,w y ) = 
e w* + e ™ v _ i anc i n{W) = {(R x , R y ): R x ,R y >0 and e~ R * + e~ R y > 1}. This 
set is unbounded, since R x can go to oo as R y — * (and vice versa). [When 
R y = this reduces to example (a).] 

7 See HO Chapters B and G], 7:'. Sections 2.3 and 3.4.3], |S2 Section II.3.8] and [HI PP- 116- 
122] for further discussion of related questions in the theory of analytic functions of several complex 
variables. 
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We now claim that these are essentially the only ways in which TZ(W) can be un- 
bounded: 

Proposition 2.16 (boundedness of 1Z(W)) Consider any repulsive lattice gas, and 
let x G X. 

(a) IfW(x, x) < I, then there exists C x < oo such that TZ(W) C [0, C x ) x [0, oo) x ^. 
[In particular, ifW(x,x) < 1 for all x G X , then 1Z(W) is bounded.] 

(b) Suppose that W(x, x) = 1, and let a G [0, oo) x ^ x . If there exists at least one y G 
X such that W(x, y) < 1 and a y > 0, then the section {R x : (R x ,sl) G IZ(W)} 
is bounded. Moreover, the converse is true provided that a G IZ(W \ (X \ x)). 



Remark. Since IZ(W) is a down-set, this proposition can alternatively be formulated 
as follows: Let a, b G [0, oo) x ; then the set {A > 0: a + Ab G 1Z(W)} is bounded 
if there exists x G suppb and y G suppa U suppb such that W(x, y) < 1; and the 
converse is true provided that a G IZ(W). 

Proof, (a) Using the fact that 1Z(W) is a down-set, it suffices to prove the claim for 
the single-site partition function 

°° «pn(n-l)/2 

Z w (w) = "j wn ( 2 - 55 ) 

n=0 

where < W < 1. If W = 0, we have Z w (w) = 1 + w, so that K(W) = [0, 1). If 
< W < 1, then Zw is a nonpolynomial entire function of order [7U1 Theorem 1.2]; 
and by the Hadamard factorization theorem (713 Theorem 1.13], any such function 
must have infinitely many zeros. In particular, TZ(W) = [0, a) where a is the smallest 
absolute value of a zero of Zw- 8 

(b) If there does not exist y G X such that W(x, y) < 1 and a y > 0, then 
Z w (w) = e Wx Z WtX \ x (w^ x ) whenever suppw C suppaU {x}. It follows that the 
section {R x : (R x , a) G TZ(W)} is either empty [in case a ^ 1Z(W \ (X \ x))] or all of 
[0, oo) [in case a G 1Z(W \ (X \ x))]. 

Now suppose that there does exist such a y. Using the fact that IZ(W) is a down- 
set, it suffices to prove the claim for the two-site partition function with A = {x, y}. 
Moreover, by the monotonicity statement in the last sentence of Proposition 12.151 we 
may assume that W(y,y) = 1. Writing W = W(x,y) G [0, 1), we need to treat 



— —w n x ^ (2.56a) 



n x ,n y =0 y 



n x =0 



^ nj 



8 It turns out (though we do not need this fact here) that all the zeros of Zw are negative real 
numbers |85l 166) . 
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If Wy 7^ 0, this is an entire function of order 1 in w x that is not of the form e aWx P(w x ) 
for any a G C and polynomial P. 9 It again follows from the Hadamard factorization 
theorem that Z w ( ■ , w y ) has infinitely many zeros. Choosing any w y G C with = 
a y (which by hypothesis is nonzero), we conclude from Theorem 12.10( c) that the 
section {R x : (R x ,a y ) G 1Z(W)} is bounded. I 



Let TZ(W) be the closure of K(W), and let &R(W) = K(W) \ TZ(W). [Note that 
this is the boundary of TZ(W) in [0, oo) x , not in M x .} Since 1Z(W) is a down-set, it 
is obvious that 1Z(W) is a down-set as well. 

Proposition 2.17 (properties of closure and boundary of 1Z(W)) For any re- 
pulsive lattice gas: 



(a) If-w G TZ(W), then 

d n Z w {w) 



dw Xl ■ ■ ■ dw Xn 



> (2.57) 



for all n > and all X\, . . . ,x n G X . 
(b) If-w G &R(W), then Z w (w) = 0. 

Proo f, (a) W e have Z w (w) > for -w G 7l(W), so by continuity Z w (w) > for 
— w G 1Z(W). The general claim ()2.57|) then follows from the multiple differentiation 
identity ()2.6|) . using the facts that < W < 1 and that 1Z(W) is a down-set. 

(b) If Zy/iyv) were strictly positive, then —1Z(W) U {w} would be a connected 
subset of (— oo, 0} x on which Zw is strictly positive, so by Theorem l2.10f a) we would 
have -w G K(W). I 

9 Proof. If P{z) = Y.k=o a kZ k , then 

e az P(z) = y \ao + n^+n(n- 1)^| + ... +n(n-l) ■■•(n- K + 1)^ 

71=0 

= q"Q(n) 
~ ^ n! 

n=0 

where Q is a polynomial. But F(n) = exp(w y W n ) is not of the form a n Q{n) for any a € C and 
polynomial Q- 
Alternate proof. We have 



Zw{w Xl w v ) - e w » = ^NK«^)-1] 



Since Wj, is fixed and < W < 1, we have | exp(wj,W / ™ a: ) — 1| < CW n * for some constant C < oo, 
from which it follows that 

\Z w {w x ,w y ) -e w *\ < Ce w ^\ . 

If Z\y(w x , w y ) were equal to e aw * P(w x ), then by taking w x —>■ +oo we conclude that we would have 
to have a = 1 and P = 1 (since W < 1); but is not in fact of this form when w y ^ 0. 
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Under an appropriate connectivity condition, we can prove that the inequality 
()2.57|) is strict for n = 1: 



Proposition 2.18 Consider any repulsive lattice gas. Let R G 1Z(W), and suppose 
that the induced subgraph Gr^suppR] is connected. Then 

<JZ,r(w) > o for all x G supp R . (2.5? 

w=-R 



Here the connectedness hypothesis is essential, as is shown by taking G to be the 
edgeless graph K n , for which Z G (w) = Y]£ =1 (l + Wi), and taking R = (1, . . . , 1). It 
is also essential to consider the induced subgraph Gyy [supp R] , since a site can be 
effectively eliminated by setting its fugacity to zero: for example, when G is the path 
123, we have -Zg(w) = (1 + u>i)(l + w^) + u>2, and a counterexample to (|2.58|) is 
obtained by taking R 2 = and R± = R% = 1. 

Before proving Proposition 12. 18l let us state and prove some corollaries. 

Corollary 2.19 (strict monotonicity) Consider any repulsive lattice gas. Let R G 
72.(W) ; and suppose that the induced subgraph GV[suppR] zs connected. If < R' < 
R wit/i i?^. < i? x /or at /east one x ; taen R' G 7£(W). 



Proof. Since Z W (-~R) > and dZ w (w)/dw x \ w= _ R > by Propositions I2.17f a) 
and I2.18( we have Zw(— (R — eS x )) > for all sufficiently small e > [here 6 X is the 
vector with xth coordinate 1 and all other coordinates 0]. Since TZ(W) is a down-set, 
we have R — eS x G 1Z(W); but by Proposition 12. 17f b) . R — eS x cannot lie in d1Z(W), 
so we must have R — eS x G 1Z(W). Now pick e small enough so that R' < R — ed x 
and use the fact that TZ(W) is a down-set. I 



Corollary 2.20 Consider any repulsive lattice gas. 

(a) If R G c?72.(VK) , inen every neighborhood of — R m (— oo,0] x contains points 
w 7^ — R wift suppw = suppR where Z w > 0, where Z w = and where 
Z w < 0. 

f6j Re 7£(W) «/ and on/y z/ Z W /(w) > for all w satisfying — R < w < 0. 

Proof, (a) Let Xi, . . . , X k be the vertex sets of the components of GV[supp R]. For 
any w with suppw C supp R, we have Z w (w) = Yii=i ^w( w a vertex 

in each JQ, and let w± = — R ± e<5 Xl + e^i=2^^r ^ Proposition 12.181 for all 
sufficiently small e > we have ±Zvy(w lxj > and Zvy(w lxj > for 2 < z < k. 
This proves the existence of points near — R with Zw > and < 0. The existence 
of points w 7^ — R with Z\y = then follows by the intermediate value theorem. 

(b) "Only if" is obvious from Proposition ^ . 1 7f a) together with the fact that 1Z(W) 
is a down-set. To prove "if", suppose that Z w (w) > for all w satisfying — R < 
w < 0. Consider the line segment {AR}o<a<i, and let A max = sup{A G [0, 1]: AR G 
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K(W)}. If A max = 1, then R G TZ{W). If A max < 1, then A max R G dTZ(W), so by 
part (a) we can choose w arbitrarily close to — A max R with supp w = suppR (and in 
particular satisfying — R < w < 0) such that Z w {yt) < — a contradiction. | 

Remark. Corollary 12.20T b) is an analogue of Theorem I2.1()f b) for the closure of 
1Z(W). It is worth noting that the analogue of Theorem 12. l()f a) is false: for example, 
for G = K 2 we have Z G (w) = (l + w 1 )(l + w 2 ), so that the set {R > 0: Z G (-R) > 0} 
is connected and equals [0, l] 2 U [1, oo) 2 , while TZ(W) = [0, l] 2 only. 

Proof of Proposition l2~T£l If R g TZ(W), then W(x,-)R e TZ{W) as well 
[since < W < 1 and 1Z{W) is a down-set], so (j2.58|) follows from the differentiation 
identity ()2.5|) . So we can assume henceforth that R G d7Z(W). By deleting sites 
where R x = 0, we can also assume without loss of generality that suppR = X. Let 
us define 

= 0} (2.59) 

w=-R 

S' = {x G X: Z w (w x , -R^.) = for all w x G C} (2.60) 

Clearly S' C S C X. The Proposition will then be an immediate consequence of 
Lemmas l2~2~D and l2~2~2l below. I 

Lemma 2.21 Let R G dTZ(W), and let S, S' be defined as in Wlty / WM) . If 
supp R = X and G\y is connected, then either S' = S = or else S' — S — X. 

Lemma 2.22 Suppose R G 1Z(W) and for all x G X and all w x G C, we have 
Zw{w x , — R^x) = 0. Then GVfsuppR] is disconnected. 

Proof of Lemma 12". 2 11 Suppose S ^ 0, and consider any x e S. Then Z w (— R) = 
and Zw(— W(x, -)R) = dZ\y(w) /dw x \ w= _ R = [using the differentiation identity 
(|2.5|) ]. By monotonicity of Z w in 1Z(W), we have Zw(— R') = whenever — R < 
— R' < — W(x, -)R. In particular, if y is adjacent to x in [i.e. W^x, y) < 1], then 
Zw(u) y , — R^y) = for all G [— -R^, — W(x, y)R y ]. By analyticity it follows that 
Zw(u) y , —H^y) = for all w y G C, i.e. y G S' C 5. Since G^y is connected, it easily 
follows that 5' = 5 = X. |' 

Proof of Lemma 12.221 As before, we may assume without loss of generality 
that suppR = X. The proof is by induction on \X\. If \X\ = 1, the hypothesis is 
impossible, since -Zvy(O) = 1 ^ 0. So let \X\ > 1 and suppose that Gw is connected. 
Choose a vertex a; G X such that Gw \ x is connected [e.g. let a; be any endvertex of 
a spanning tree in Gw}- By the fundamental identity ()2.Hj) . we have 

^(wjj, -R^c) = iV(0, -Kjt x ) + w x Z w (0, -W(x, -)R^) 

+ terms of order w 2 x and higher (2-61) 



S = {x G X: 



dZ w {vj) 

dw x 
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Since Zw{', — R^r) is identically zero, all coefficients in the Taylor expansion (j2.61|) 
must vanish; in particular we must have 



Z w (0,-R^ x ) = (2.62a) 
Z w (0,-W(x,-)R^ x ) = (2.62b) 

Pick y ^ x with W(x, y) < 1 [which is possible since x is not an isolated vertex of 
Gw\- Since R y > and Zw(0,w y , — R^ x , y ) is monotonic in w y for w y G [— .R^O], it 
must be vanishing on the nontrivial interval [—R y ,—W(x,y)R y ]. By analyticity of 
Z w , it follows that Zyy(0, w y , — R^ X) y) = 101 a ll w y £ C, so in particular 

Z w (0,-Ry,-R^ Xty ) = (2.63a) 
0^(0, w y ,-R^ y ) _ Q (2 63b) 



9w y 



Notice that R = (i^, R^ x ) G 7£(W) implies (0, R^) G 7£(W) [since ft(W) is a 
down-set], or in other words R^ x . G ^(W 7 f (X \x)). Moreover, by (j2.63a|) R^ x must 
lie in d1Z(W \ (X \ x)). So we can apply Lemma 12.211 to X \ x to conclude that 
S' = X \ x; therefore, by the inductive hypothesis, Gw \x is disconnected, which 
contradicts the choice of x. I 



Remark. One might think that Lemma 12.221 would hold not only for R G TZ(W) 
but for arbitrary w(°) G C x . However, this is false: Consider the star K 14: with 
center x and endvertices yx, . . . , with interactions W G (0, 1) along the edges [and 
W{x,x) = W{y uyi ) = 0]. Then 

4 4 

Z w (w x ,w yi , . . . ,w y4 ) = i\p- + w yi ) + w x JJ(1 + Ww yi ) . (2.64) 

i=i i=i 

Setting w y °i = wf^ = — 1 and w y ^ = w y } = —1/W [here w x ^ is arbitrary], we see 
that Zw{w z , w^l) = for every vertex z and every w z G C. 



2.5 Further consequences of the alternating-sign property 

Let us now exploit systematically the alternating-sign property (J2.29|) for the 
Taylor coefficients of logZw- The general context is the following: 

Definition 2.23 (absolute monotonicity) Let U C C" be a union of open poly- 
discs centered at (a "complete Reinhardt domain"). We say that f is absolutely 
monotone in U if it is analytic in U and all the Taylor coefficients of f at are 
nonnegative. 



We will use the following essentially trivial result: 
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Lemma 2.24 (elementary consequences of absolute monotonicity) Let f be 

absolutely monotone in a union U of open poly discs centered at G C n . Then: 

(a) All the derivatives D a f (where ex. is a multi-index) are absolutely monotone in 
U. 



(b) For every A G [0, l] n , the function 

A(z) = /(z) - /(Az) 



(2.65) 



is absolutely monotone in U. [Here Az is the pointwise product, so that (Az) x = 

(c) For all multi-indices a, we have 

\(D<*f)(z)\ < (D<*f)(R) (2.66) 
whenever |z| < R G U. In particular, (£)"/) (R) > whenever < R G U. 
Now, the alternating-sign property ()2.29|) is precisely the statement that 

/(z) = -logZ w (-z) (2.67) 
is absolutely monotone in the domain D-ji^ = |J D R = |J D R . Applying 
Lemma f2.24f c) to /, we obtain: 

Proposition 2.25 For any repulsive lattice gas: 
(a) When — w G 7Z(W), we have 



! <9"logZ^(w) v () 



9u; xi • • • dw Xri 



(2.68) 



/or all n > and all X\, . . . ,x n G X. 
WTien |w| < R G 7£(W), we nave 



<9 n logZ w (w) 



< - 



dw Xl ■ ■ ■ dw Xn 
for all n > and all X\, . . . , x n G X. 



<9u>' ■ • • 

Xl ^"^2 



(2.69) 



w'=-R 



Proof. This is an immediate consequence of the alternating-sign property (|2.29|) 
together with Lemma f2.24f c) and Theorem 12. lUf d) . ■ 



Let us now define 



Mw;A) 



Z w (\w) 



(2.70) 



29 



By Lemma f2.24f b). the function 



fx(z) = loglV(-z;A) = log' 



Z w (-z) 



(2.71) 



is absolutely monotone on Dft/w) whenever A G [0, l] x . Applying Lemma l2.24f c) to 
fx, we obtain: 

Proposition 2.26 For any repulsive lattice gas and any A G [0, l] x : 
(a) When — w G 7Z(W), we have 



d"loglV(w;A) ^ Q 
■ ■ ■ dw Xn 

for all n > and all X\, . . . , x n G X. 
(b) When |w| < R G TZ(W), we have 



(2.72) 



d n logY w (w; A) 



• • • 8w Xri 



< 



d n logY w (w';X) 
dw' Xl ---dw' Xn w , = _ R 

for all n > and xj, . . . , x n G X. In particular, setting n = we have 



2.73) 



Z w (Xw) 



Z w (w) 



< Zw{— AR) 
R) 



(2.74) 



/or a// w G Dr. T/iai is, the quantity \Zw(Xw) /Zw(w)\ takes its maxi- 
mum on the polydisc Dr at the point w = — R, and this maximum value 
Z w (— XR)/Z W (— R) is therefore an increasing function o/R G 72.(W). 



Let us define the function Zw'- [0, oo) x — > [0, oo) by 
Z W (R) = inf max[Z w (w),0] 

— R<w<0 

= inf |Z w (w)| 



Z w (-R) for R G 7fc(W) 
for R £ 7fc(W) 



(2.75a) 
(2.75b) 

(2.75c) 



where the equivalence of these three expressions is an immediate consequence of 
Theorem 12.101 We then have: 

Corollary 2.27 For any repulsive lattice gas, the function Zw is 

(a) continuous 

(b) decreasing 
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(c) satisfies Ziy(AiA 2 R) Z W (R) < Zw^R) Z W (\ 2 R) for all X U X 2 G [0, 1} X 
and R G [0, oo) x 

(d) log submodular [i.e. Z w (Ri A R 2 ) Z W (R\ V R 2 ) < Z w {Ri) Z W (R 2 ), where A 
(resp. y) denotes the elementwise min (resp. max) of vectors] 

Proof, (a) The continuity of Zyy follows easily from ()2.75a|) or (j2.75b|) and the 

continuity of Zjy. [For any r < oo, the function Z\y is uniformly continuous on the 

compact ball B = {w G C: |w| < r}. If R, R G B n [0, oo) x with |R- R| < 5, then 

|Zyi/(R) — Zi4/(R')| < sup |Zyp(w) — Z w (w')\. Uniform continuity of Z w on 

w, w' e b 

|W - W'| < 5 

B therefore implies uniform continuity of Zw on B D [0, oo) x ] 

(b) The decreasing property is an immediate consequence of (J2.75a|) or (|2.75b|) . 

(c) This is trivial if R £ K(W), so assume R G ft(W). Then 3?V(-R; A 2 ) = 
Zvy(A 2 R)/Zvy(R), while 1V(— AiR; A 2 ) = Z^(AiA 2 R)/Zvi/(AiR), so that the claim 
follows from Proposition 12.261 with n — 1. 

(d) The log submodularity follows immediately from (c) by setting R = Ri V R 2 
and Aj = Rj/R for i — 1,2 (setting 0/0 = 1 where needed). [Alternatively, for 
Ri V R 2 G 1Z(W) it follows by integrating (j2.68J) with n = 2 and x\ ^ x 2 ] ■ 

We can also prove a strict version of the cases n = of Propositions I2.25f a) and 
Ell a) as well as Corollary fTTft d. Note first that 

-logZ w (-K) = ^(-f^c^jR" (2.76) 

n 

-\ogZ w (-R) + \ogZ w (-XR) = ^(-l^-^WOa-A^R" (2.77) 

n 

-logZ w (-R) + logZwr(-AiR) + \ogZ w (-X 2 R) - logZ w (-AiA 2 R) 

= ^(-l) |nhl c n (W) (1 - A?)(l - A 2 ) R n (2.78) 

n 

By Proposition I2.8[ all three quantities are nonnegative whenever R G ^(W 7 ) and 
A, Ai, A 2 G [0, l] x . We can determine when they are strictly positive: 

Proposition 2.28 Consider any repulsive lattice gas, and let R G 1Z{W) and A, Ai, A 2 G 
[0, l] x . Then: 

(a) — log Z w (—R) > if and only if 0. 

(b) — log Zw(— R) + log Zw(— AR) > if and only if supp R PI supp(l — A) ^ 0. 

(c) -logZ VF (-R) + logZ l4 /(-AiR)+logZ w (-A 2 R)-logZ W /(-AiA 2 R) > if and 
only there exists a component of Gw [supp R] that meets both supp(l — Ai) and 
supp(l - A 2 ). 
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PROOF. We examine the terms in ()2.76j) - ()2.78j) and determine (for some of them) 
when they are nonzero. 

(a) "Only if" is trivial; for "if" , just consider the term n = S x for some x G supp R, 
and use Proposition 12. 9( a). 

(b) "Only if" is again trivial; for "if" , just consider n = 8 X for some x G supp R PI 
supp(l — A), and again use Proposition 12.9( a). 

(c) For "if" , let S be the vertex set of a component of Gw [supp R] that meets both 
supp(l — Ai) and supp(l — A 2 ), take n = l s , and use Proposition 12.9( a) or 12.9( c). 
For "only if", suppose that c n (W) ^ 0. Then, by Proposition 12.91 either n = S x 
for some x G X, or n = kS x for some k > 2 and some x G X having W(x,x) ^ 1, 
or | supp n | > 2 and GV [suppn] is connected. In all three cases, GV [suppn] is 
connected. Furthermore, to make a nonzero contribution to ()2.78|) . suppn has to 
meet both supp(l — Ai) and supp(l — A 2 ) and be contained in suppR. Therefore, 
there is a component of GVfsuppR] (namely, the one containing suppn) that meets 
both supp(l — Ai) and supp(l — X 2 ). I 

For a repulsive lattice gas with hard-core self-repulsion [i.e. W(x, x) — for all 
x G X], the set TZ(W) has an additional convexity property: 

Corollary 2.29 (convexity of two-dimensional sections) For any repulsive lat- 
tice gas with hard-core self-repulsion, the two-dimensional sections ofTZ(W) parallel 
to the coordinate axes are convex and (when nonempty) are bounded either by a hy- 
perbola or by a straight line. 



Proof. For x, y G X (x ^ y) and a nonnegative vector a = (a z ) ze x\{x,y}, consider 
the section 

P x , y , a n(W) = {(R x ,R y ): (R x ,R y ,a) G K{W)} . (2.79) 

If (0,0, a) ^ 1Z(W), then clearly the section is empty; so let us assume henceforth 
that (0, 0, a) G TZ(W). Let Zw(w x , w y ; —a) be the polynomial obtained from Zy^(w) 
by setting w z = —a z for z ^ x,y. Then, as Zw is multiaffme, we have 

Zw(w x , w y ] —a.) = A + Bw x + Cw y + Dw x w y (2.80) 

where 



A = Z w (0,0;-a) (2.81a) 

B = Z w {0, 0; -W{x, •)<*) (2.81b) 

C = Z w (0, 0; -W(y, -)a) (2.81c) 

D = W(x,y)Z w (0,0;-W(x,.)W(y,-)a) (2.81d) 

as a consequence of the multiple differentiation identity (|2.fij) . Since (0, 0, a) G 1Z(W), 
< W(x, y) < 1 for all x, y, and 1Z(W) is a down-set, we have A, B, C > and D > 
[with D = if and only if W(x, y) = 0]. Moreover, by Corollary 12.27( c) we have 
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AD < BC. W It follows from Theorem EHUl a) that P x>y>oc Tl(W) is the component of 
[0, oo) 2 fl {(R x ,Ry): Z w (—R x ,—R y ;—et) > 0} containing (0,0) [here we have used 
the fact that (0,0, a) E TZ(W)]. Thus, 

P x , y , a n{W) = {(R x , R y ): R x , R y > and BR X + CR y < A} (2.82) 

in case D = 0, and 

Px, y ,<Ji{W) = {(R x , R y ): 0<R X < C/D and < R y < B/D and 

(C/D - R X )(B/D - R y ) > (BC - AD)/D 2 } (2.83) 

in case D > 0. I 

Remarks. 1. The proof shows that the boundary of P x ,y,ctR,(W) is a straight line 
when W(x, y) = and is a hyperbola when < W(x, y) < 1. 

2. Consider w = (w x ,w y ) and w' = (w' x ,w' y ), and let < A < 1; then 

Z w (Xw + (1 - A)w'; -a) - [XZ w (w; -a) + (1 - X)Z w (w'; -a)} 

= _ jDA (l - X)(w x - w' x ){w y - w' y ) , (2.84) 

so that Zw( ■ ; —a) is concave along line segments of negative slope in the (w x ,w y )- 
plane. This provides an alternate proof that the projected section P Xj y )C J^(W) is 
convex [because 1Z(W) is a down-set, only line segments of negative slope need be 
considered] . 

3. The convexity of two-dimensional sections is false without the hypothesis of 
hard-core self-repulsion: consider X = {1,2}, W(l, 1) = W(2,2) = 1 and W(l,2) = 
0. Then Z w (w) = e Wl + e W2 - 1, so that K{W) = {(R u R 2 ): Ri,R 2 >0 and e - Rl + 
e -R-2 > l} ) which is not convex. 

4. Even in the case of hard-core self-repulsion, TZ(W) itself is not in general 
convex. To see this, let X = {1, 2, 3} and let G be the path 123; its independent-set 
polynomial is Zq(w) = W2 + (1 + wi)(l + "U^)- The diagonal section of 7Z(G) given 
by R\ = i?3 = x and R 2 = y is the region of [0, l] 2 where y < (1 — x) 2 , which is not 
convex. 

2.6 Algebraic irreducibility of Z^(w) 

In this subsection — which is a digression from the main thread of the paper and 
can be omitted on a first reading — we discuss the algebraic irreducibility of the 
multivariate partition function Z^(w). We restrict attention to the case of hard-core 
self-repulsion, in which Zw is a multiaffine polynomial. 

Let R be a commutative ring with identity, and let W: X x X — > R be symmetric 
and satisfy W(x,x) = for all x e X. Define the support graph Gw by setting 

10 In the case of a hard-core pair interaction, this follows alternatively from Corollary 12 . 27( d) . 



33 



xy G E(G) if and only if W(x, y) 7^ 1 and x 7^ y. Define the polynomial Zw G R[w] 
by 




(2.85) 



Proposition 2.30 Suppose that R is an integral domain. Then Zw is irreducible 
over R if and only if Gw is connected. 

The proof is based on an easy but crucial lemma f3H Lemma 4.7]: 

Lemma 2.31 Let Pi and P2 be nonzero polynomials in the indeterminates {w x } x ex 
with coefficients in an integral domain R. Suppose that P1P2 is multiaffine (i.e. of 
degree or 1 in each variable separately). Then: 

(a) There exist disjoint subsets Xi,X 2 C X such that Pi uses only the variables 
{w x }x£Xi (i = l,2). 

(b) Pi and P 2 are both multiaffine. 

Proof. Suppose there exists x G X such that both Pi and P 2 use the variable w x . 
For % = 1, 2, let dj > 1 be the degree of Pj in the variable w x , and let Qi 7^ be the 
coefficient of in Pj, considered as an element of the polynomial ring P[w^J. Then 
Q1Q2 7^ because P[w^J is an integral domain |J51 Theorem III. 5.1 and Corollary 
III. 5. 7]. But this shows that the coefficient of w^ 1+d2 in P1P2 is nonzero, contradicting 
the hypothesis that P1P2 is multiaffine (since di + d 2 > 2). This proves (a); and (b) 
is an easy consequence. I 

Proof of Proposition 12.301 If G w is disconnected, it is obvious that Z w is 
reducible. To prove the converse, suppose that Zw is reducible, i.e. Zw = P1P2 
where Pi and P 2 are nonconstant polynomials over R. It follows from Lemma 12.311 
that there exist disjoint subsets Xi,X 2 C X such that Pj uses only the variables 
{w x } x€ Xi (we can assume without loss of generality that Xi U X 2 = X); moreover, 
neither Xi nor X 2 can be empty (since Pi and P 2 are nonconstant). Let 6 R be 
the constant term of Pj [i = 1, 2). Since Zw has constant term 1, we have aia 2 = 1. 
Replacing Pi by a 2 Pi and P 2 by CI1P2, we can assume without loss of generality that 
ai = a 2 = 1. For x G Xi (resp. x G X 2 ), let b x be the coefficient of the linear term w x 
in Pi (resp. in P 2 ). Since the linear term w x has coefficient 1 in Zw and ai = a 2 = 1, 
we have b x = 1 for all x G X. It follows that, for x E Xi and y G X 2 , the term w x w y 
in P1P2 has coefficient 1; since this term has coefficient W(x,y) in Zw = P1P2, we 
have W(x, y) = 1 whenever x G X\ and y G X 2 . Hence Gw is disconnected. I 

Remarks. 1. The univariate polynomial Zw{w) is obviously reducible over C 
whenever \X\ > 2; and it is in many cases reducible over ]R as well. Indeed, it 
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can be reducible over the integers: e.g. for G = P4 (the 4-vertex path) we have 
Z G {w) = 1 + Aw + 3w 2 = (3w + l)(w+ 1). 

2. Results analogous to Proposition 12.31)1 hold for some other combinatorial poly- 
nomials. For example, the irreducibility over C for the (bivariate) Tutte polynomial 
Zm{x,d) of a connected matroid M was proven recently by Merino, de Mier and 
Noy [80J. Likewise, the irreducibility over any integral domain R for the multivariate 
basis generating polynomial of a connected matroid M — or more generally, for any 
multiaffine polynomial whose support is exactly the collection of bases of M — was 
proven recently by Choe, Oxley, Sokal and Wagner 

3. It would be interesting to know whether analogous results can be obtained 
without the hypothesis of hard-core self-repulsion. The difficulty is to find the right 
ring in which to work. Irreducibility in the ring i2[[w]] of formal power series is too 
trivial, because every formal power series with constant term 1 is invertible. On the 
other hand, it might conceivably be possible to prove irreducibility in the ring of 
entire functions (or entire functions of order 1) on <C X — but only for \X\ > 2, since 
it is manifestly false if \X\ = 1 and < W(x, x) < 1 (a nonpolynomial entire function 
of order has infinitely many zeros [THJ Theorem 1.13]). 

2.7 Convexity of log Z at nonnegative fugacity 

In this paper we are primarily concerned with the behavior of Zw(w) for complex 
fugacities w. However, the regime of nonnegative fugacities w is of particular interest 
to probabilists and statistical mechanicians, since the Boltzmann weights [cf. (jl.lj) ] 
are there nonnegative and so can be interpreted, after normalization, as a probability 
measure on lattice-gas configurations. In this "probabilistic regime", the logarithm 
of the partition function is, in very great generality, a convex function of all the 
interaction energies jHZl P- 12]. Let us prove the specialization of this statement to 
our model: 

Let < W,W < 1 and w,w' > 0. Then, for < A < 1, define W x (x,y) = 
W(x, y) x W'(x, yY~ x and (w\) x = w x {w' x ) l ~ x . By Holder's inequality applied to (|l.laj) 
or (jl.lbj) . we obtain: 

Lemma 2.32 Suppose that < W, W < 1 and w, w' > 0. Then, for < A < I, 

ZwA^x) < Zwiw^Zw^w') 1 -* . (2.86) 

That is, log Zw x (w\) is a convex function of A; or in other words, log-ZV(w) is a 
convex function of the vector ({log W(x, y)} x ,y€Xi {log w x } x€ x)- 
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3 The lattice gas with hard-core self-repulsion 



We now restrict attention to the case of a repulsive lattice gas with hard-core 
self-repulsion, i.e. < W(x, y) < 1 for all x, y and W(x, x) = for all x. This means, 
as noted in the Introduction, that Zvk(w) can be written as a sum over subsets: 

Z w (w) = 

IT " < IT w ( x >y)- (3.i) 

X'QXxGX' {x,y}CX' 

Here we shall exploit this special structure, which implies that Z\y is a multiaffine 
polynomial , i.e. a polynomial of degree 1 in each w x separately. Since W will be fixed 
throughout, we shall henceforth often omit it from the notation. 

3.1 The fundamental identity 

Let us define, for each subset A C X, the restricted partition function 

Z A (w) = E IT " II W(x,y). (3.2) 

X'CA x£X> {x,y}CX' 

Of course this notation is redundant, since the same effect can be obtained by setting 
w x = for x G X \ A, but it is useful for the purpose of inductive computations and 
proofs. We have, for any x G A, the fundamental identity 

Z A (w) = Z A \ x (w) + w x Z A \ x {W(x, -)w) (3.3) 

where 

[W(x,-)w] y = W(x,y)w y ; (3.4) 

here the first term on the right-hand side of (J3.3)) covers the summands in (|3.2|) with 
X' $ x, while the second covers X' 3 x. [Note that (|3.3|) is a special case of (|2.3j) .] 
In the special case of a hard-core interaction (= independent-set polynomial) for a 
graph G, ()3.3|) reduces to 

Z A (w) = Z A \ x (w) + w x Z A \ r * {x) (w) , (3.5) 

where we have used the notation T*(x) = T(x) U {x} [here T(x) denotes the set 
of vertices of G adjacent to x\. The fundamental identity ()3.3|) / ()3.5|) will play an 
important role both in the inductive proof of the Lovasz local lemma [cf. (J4.1(J|) and 
(I4.22|) ] and in the Dobrushin-Shearer inductive argument for the nonvanishing of Zy/ 
in a polydisc (Section EJ). 

Let us also remark that because Z A \ X is multiaffine, the last term in (|3.3j) can be 
expanded to rewrite Z A \ x (W(x, -)w) as a linear combination of values at the vertices 
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of a rectangle: 11 

z A (w) = z AV (w) + Wx ( II ^ y))\ II t 1 - w ( x > ^)] I ■ 

(3.6) 

In the hard-core case, the only term with a nonzero coefficient is Y = A \ T*(x), so 
that ()3.6|) reduces to ()3.5|K 

Remark. Repeated use of (|3.3|) obviously gives an algorithm to compute Z w (w). 
But this algorithm takes in general exponential time. In fact, calculating Zc(w) for 
general graphs G (or even for cubic planar graphs) is NP-hard (as noted by Shearer 
|98j ) . since even calculating the degree of Zq(w) — that is, the maximum size of an 
independent set — is NP-hard |4*H pp. 194-195]. Therefore, if P 7^ NP it is impossible 
to calculate Zg(w) for general graphs in polynomial time. 

A key role in (j3.3J) is manifestly played by the rational function 

Z AV (W(x,-)w) 
Za\ x {w) 

Note that JC X; \(w) depends only on {w y } y< =A\ x ', m particular, it does not depend on 
w x . We shall sometimes write it as K. Xt \(w^ x ) to emphasize this fact. Please note 
also that Z^(w) is nonvanishing in the region \w x \ < 1/|/C Xj a( w ^x)| and vanishes at 
w x = -1//C I; a( w ^)- We have 

dlogZ A (w) = JC x ,a{^x) / 3g x 
dw x 1 + JC X} a(w^ x )w x ' 

Since /C Xj a( w ) * s a special case of the function 1V( W ; A) defined in ()2.70|) . we obtain 
as an immediate corollary of Proposition 12. 26f b) : 

Proposition 3.1 Consider any repulsive lattice gas with hard-core self-repulsion. If 
ReK(W), then 

\JC xA M\ < /C X , A (-R) < 00 (3.9) 

for all w e Dr. That is, the maximum of |/C Xj a(w)| over Z)r is attained atw = — R ; 
and this maximum value )C X ^(— R) is an increasing function of HE 1Z(W). 

n If F: C x — > C is multiaffme, then F can be reconstructed from its values at the corners of any 
rectangle with positive volume. For example, if the rectangle is [0, l] x , we have 



f ( z ) = e (n*j ( n a-^po 

YCX \i£Y ) \i(iX\Y J 



This is clear when z g {0, 1} , and equality for all z then follows because both sides are multiaffme. 
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The function K. X) \(— R) will play an important role in Section 

We can translate the fundamental identity (|3.3|) into a recursion for the rational 
functions JC x ^(w). Let us order arbitrarily the sites A\x = {yi, . . . ,yk}- Then we 
write (jSZZD telescoping product: 

where the vectors w^' are defined by 

(w (i) ) = { W ^ Xl V ^ Wy iiy = Vj f ° r S ° me J ~ 1 (3 11) 

1^2/ otherwise 

[Note that it is irrelevant how we define (w''') 9 for y = x, since we set this weight to 
zero anyway by considering Z\\ x .] Applying the fundamental identity ()3.3|) to y^, we 
obtain 

^a\x(w«) 1 + W(x }yi )w y X mAV (w^) 
^(wC-D) " ' 1 + ^^(w^i)) 1 • j 

and hence 

Uw) ^ LljM^) , (3,3, 

This identity will play a central role at the end of Sectional In the special case of a 
hard-core pair interaction, ()3.13|) becomes 

' 1 
/C,,a(w) = ft— (3.14) 



1=1 



where {yi, . . . ,yi} is an ordering of A n T(x). 

Note, finally, that the partition functions Z\(w) can be reconstructed from the 

rational functions JC x ^(w). For, by the fundamental identity ()3.3|) and the definition 
([3.7)1 . we have 

IW) 

„ V \ = 1 + w x JC xA (w^ x ) ; (3.15) 
and hence, setting A = {x\, . . . ,x n } (in arbitrary order) we can write 

Za(w) = Jl^ilW = I] [1 + <r,X,.,,, ,,!wj; . (3.16) 

i=i *{*i,~Pi-in w ) i=1 



3.2 Examples 

Let us take a moment to compute a few examples of partition functions Z^(w). 
We shall assume hard-core self-repulsion [i.e. W(x,x) = for all x] throughout, 
and shall mostly restrict attention to the case of a hard-core pair interaction, i.e. 
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the independent- set polynomial Zg(w) for a graph G. In some cases we shall, for 
simplicity, compute only the univariate polynomial Z G {w) obtained by setting w x = w 
for all vertices x. 

Example 3.1. The complete graph K n . Clearly Zx n (w) = 1 + Wi + . . . + w n . In 
particular, 1Z(K n ) = {R: R\ + . . . + R n < 1}. 

Example 3.2. The n-vertex path P n . We limit attention to the univariate 
independent-set polynomial. Applying the fundamental identity ()3.5|) to an end- 
vertex, we obtain the recursion relation 

Z Pn (w) = Z Pn _ x [w) + wZ Pn _ 2 (w) , (3.17) 

which is valid for all n > if we define Zp = Zp _ 1 = 1 and Zp_ 2 = 0. The solution 
of (jgTZD is 

Z Pn (w) = -=L= (Xl^ - A^ +2 ) (3.18) 
VI + 4w 

where 



(An alternative method of obtaining this result is given in Example 3.2' below.) The 
zeros of Zp n are located at 

w = 'I T^k- ( 3 - 2 °) 

4cos 2 ^ 

for k — 1,2, ... , L^y^J • The zero nearest the origin converges from below to w — — 1/4 
as n —>■ oo. Elementary counting arguments give also the explicit formula 



fc=0 



Example 3.3. The n-vertex cycle C n . We again limit attention to the univariate 
independent-set polynomial. Applying the fundamental identity ()3.5|) . we find 

Z Cn (w) = Z Pn _Sw) + wZ Pn _ 3 (w) , (3.22) 

valid for n > 2. Inserting ()3.18|) / (|3.19J) . we obtain 

Z Cn (w) = X n + + X n _ . (3.23) 

(An alternative method of obtaining this result is given in Example 3.3' below.) The 
zeros of Zc n are located at 

w = — Wy (3 - 24) 
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for k = 0, 1, ... , [f J — 1- The zero nearest the origin converges from below to w = 
— 1/4 as n — > oo. Elementary counting arguments [or ()3.21j) / ()3.22j) ] give also the 
explicit formula 

~" J ,., 

w k . (3.25) 

n — i: \ I;' I 

k=0 



L-J 

„ . , n fn — k\ 

= Esn( k ) 

u—n \ / 



Remark. There is a one-to-one correspondence between matchings on a graph G 
and independent sets on the line graph L(G). Since L(P n ) = P„_i and L{C n ) = C n , 
these formulae for Zp n and Zc n can also be found in the work of Heilmann and Lieb 
on matching polynomials [5B1 pp. 196-197]. Heilmann and Lieb also noted that Zp n 
and Zc n can be written as Chebyshev polynomials. 

Example 3.3'. The "soft" n-vertex cycle C n . Let us now generalize Example 3.3 
by considering the cycle C n with edge weights 

jxrt \ \W if xy is an edge of C n fo 
W(x,y) = < y . 6 n (3.26) 

i- 1 otherwise 

where W is a constant in [0, 1] (the case W = corresponds to the independent- 
set polynomial of C n , while the case W = 1 corresponds to the independent-set 
polynomial of the edgeless graph on n vertices). Then the univariate polynomial 
Zw(w) is given by a "transfer matrix" Section 2.1] [TU Sections 2.2-2.4] |107| 
Section 4.7]: 

= tr (I JwJ < 3 ' 27 > 

where the first (resp. second) row or column corresponds to an empty (resp. occupied) 
site, and the fugacity w is attributed to the row only. It follows that 

Z w (w) = A™ + \ n _ , (3.28) 

where 

A± = i + ww± v /(i-wwy + 4w (3 2g) 

are the eigenvalues of the transfer matrix. Simple algebra then shows that all of the 
zeros of Z\y are real and negative (when < W < 1). 

Example 3.2'. The "soft" n-vertex path P n . Let us now generalize Example 3.2 
analogously. The univariate polynomial Zyy(w) is again given by a transfer matrix: 

i i \"~Vi 



Z W ( W ) = (1 1) ^ wW j j . (3.30) 

After some algebra we find 

Z w (w) = A + X\ +2 + A_A™ +2 , (3.31) 
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with A± given by ()3.29|) and 

A ± = 2 = . (3.32) 

±(2 - W + wW 2 )y/(l-wW) 2 + 4w + W[(l - wW) 2 + Aw] 

Example 3.4. The star K\^ T . Let x be the center vertex of the star and let 
yi, . . . ,y r be the leaves. The independent-set polynomial is 

r 

i=l 

More generally, suppose we have a pair interaction matrix W(x, yi) [with W(yi, yj) = 1 
for all i 7^ j and W(x, x) = W(yi, yi) = 0]. Then 

r r 

Z w (w) = ]J(l + w yi ) + w x Y[[l+W(x, yi )w yi \. (3.34) 

i=i i=i 

Example 3.5. The complete bipartite graph K m>n . Let X = {x\, . . . ,x m } and 
Y = {yi, . . . , y n } be the bipartition. Then the independent-set polynomial is 

m n 
i=l j=l 

Example 3.6. The complete r-ary rooted tree |94[ I98j. Let Tn be the complete 
rooted tree with branching factor r and depth n. We limit attention to the univariate 
independent- set polynomial. Fix r > 1; and to lighten the notation, let us write Z n 
as a shorthand for Z _( r ) . Applying the fundamental identity (13.51) to the root vertex, 
we obtain the nonlinear recursion 

Z n (w) = Z n _i(w) r + wZ n ^ 2 (w) r2 , (3.36) 

which is valid for all n > if we set Z-\ = Z-i = 1. By defining 

Y n (w) = ^) {337) 
Z n _ x [w) r 

we can convert the second-order recursion ()3.36|) to a first-order recursion 

w 

y n _i(w) r 

with initial condition YLi = 1. The polynomials Z n (w) can be reconstructed from 
the rational functions F„(w) by 

n 

z„h = n^H r ""- ( 3 - 39 ) 

fc=0 
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Let w n < be the negative real root of Z n of smallest magnitude (set w n = — oo 
if Z n has no negative real root). Note that w_\ = — oo and Wq = —1. Let us prove 
by induction that w n _i < w n for n > 0. It is true for n = 0. For n > 1 we have 

2 

Z n (w n _i) = Z„_ 1 (w n „i) r + w„-i^„- 2 K-i) r < (3.40) 

since Z„_ 1 (w n _ 1 ) = 0, w n _i < and Z n _ 2 (u> n _i) > by the inductive hypothesis. 
Therefore Z n vanishes somewhere between u>„_i and 0. 

It follows that the w n increase to a limit Woo < as n — > oo. Let us show, following 
Shearer [§E], that 

W ~ = - ( r + l)r + l ^ 3 ' 41 ) 

by proving the two inequalities: 

Proof of >: If w G [wooiO), we have Z n {w) > for all n and hence also 
Y n (w) > for all n. Since Ylj > F , it follows from the monotonicity of (|3.38|) that 
{Y n (w)} n >o is a strictly decreasing sequence of positive numbers, hence converges 
to a limit > satisfying the fixed-point equation = 1 + w/yl, or equivalently 
w = yl +1 — yl- Elementary calculus then shows that w > —r r /{r + l) r+1 ; taking 
w = Woo we obtain > — r r / (r + l)' r+1 . 

Proof of <: If — r r / (r+ l) r+1 < w < 0, the equation w = — yl has a unique 
solution y* G [r/(r + 1), 1). It then follows by induction [using ()3.38|) and the initial 
condition K_i = 1] that 1 = Y^w) > Y (w) > . . . > K„_i(u>) > Y n {w) > . . . > y* for 
all n > 0. In particular, Y n {w) > for all n, so that w n < w for all n. This shows 
that Woo < -r r /(r + l) r+1 . 

Let us conclude by observing that (|3.38|) defines a degree-r rational map R w : y i— > 
l+u>/V" parametrized by w G C\0. Moreover, the zeros of Zy/{w) correspond to those 
values w for which R w has a (superattractive) orbit 0i— > oo i- > 1+w t— >0 
of period + 3 (or some submultiple of n + 3). As n — > oo, these points accumulate 
on a "Mandelbrot-like" set in the complex w-plane. For further information on the 
maps y i-> 1 + w/|/ r , see jHH El UHl CSl 123 EH • 

3.3 Reduction formulae 

Recall that, given W, we have defined a simple loopless graph G = Gw ("the 
support graph of W") by setting xy G -E'(G') if and only if W(x, y) ^ 1 and x ^ y. 
When the support graph G has a very simple structure, the polynomial Z\y(w) can 
be simplified: 

1) Disconnected graph. Suppose that G can be written as the disjoint union of G\ 
and Set A, = V{Gi). Then clearly Zw factorizes: 

^Axua 2 (w) = Z Al (w) Z A2 (w) . (3.42) 

Note that Za^w) depends only on {w x } x& \ i . 

2) Cut vertex. Suppose that G = G x U G 2 where V{G X ) n V(G 2 ) = {^}. [If 
and G2 have at least two vertices, this means that a; is a cut vertex of G. But the 
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formulae below hold also in the trivial cases V{G\) = {x} or V(G 2 ) = {x} or both.] 
Again set A, = V(Gi). The fundamental identity ()3.3|) asserts that 

^Aiua 2 (w) = Z (AlUA2) \ x (w) + w x Z {AlUA2) \ x (W(x,-)w) . (3.43) 

But since the graphs Gi\x and G2 \ x are disjoint 12 , we can apply the factorization 
(13.42(1 to the right-hand side of ()3.43|) . yielding 

^AiUA 2 (w) = Z Al \ x (w) Z A2 \ x (w) + w x Z Al \ x (W(x, -)w)Z M \ x (W(x, -)w) 

(3.44a) 

= Z AAx (w)[Z A2V (w) + wf^Z A2 \ x (W(x,-)w)] (3.44b) 
= Z Al \ x (w)Z A2 (w^ x ,wf^) (3.44c) 
where we have defined 

eff,Gi Z Al \ x (W(x, -)w) 
w x ,U1 = w x — = w x lC xM {w) . (3.45) 

^AxVHW) 

This can be interpreted as "integrating out" the variables in V(G\) \ x, leaving an 
"effective fugacity" w e x ,Gl for the vertex x G V(G 2 ). 

Example 3.4 revisited. Integrating out the leaves y±, . . . , y r of the star Ki >r , we 
find from fl3~M|) : 

1 + i y(x,^ )m 
+ w 



= II k W • (3-46) 

i=1 



3.4 When G is a tree . . . 

When the support graph G = Gw is a tree, every vertex of G is either a cut vertex 
or a leaf, so the partition function Z w (w) can be calculated by repeated use of the 
reduction formula (|3.44(l / (|3.45(1 . This can be done in many ways, but the simplest is 
probably to "roll up" the tree "from the leaves up" , as follows: 

We say that a vertex x G V(G) is a near-leaf of G if it is adjacent to at least 
one leaf of G. If x is a near-leaf of G, then we can write G = G\ U G2 where 
V(Gi) fl V(G 2 ) = {x}, G\ is a star with center x consisting of x and some or all of 
the leaves of G adjacent to x, and G 2 is a tree. Integrating out G%\x using (|3.46(1 . 
we obtain 

yEV(Gi\x) y 

and of course 

^aa.(w) = n + ( 3 - 48 ) 
yeV(Gi\x) 

12 We recall that G\x denotes the graph obtained from G by deleting the vertex x and all edges 
incident with it. 



J] t^l WW (3.47) 

1 ~\~ VJy 
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Applying this process repeatedly, we can obtain a linear-time algorithm for evaluating 
Z w (w) at any point w G C x whenever the support graph Gw is a tree. (It is also 
an algorithm for computing Zyy(w) as a polynomial in w, but in this context it may 
no longer be linear-time, because of the need to multiply rational functions of w.) 
The simplest approach is to orient the tree by choosing arbitrarily one vertex of the 
tree as the root, and then to work upwards from the leaves. We obtain the following 
algorithm: 

Algorithm T. Given numbers {w x } x( zx and {W(x,y)} Xiye x for which 
the support graph G = Gw is a tree, we compute Z\y{w) as follows: 

1) Pick (arbitrarily) a root Xq G X. Define 

depth[x) = dist(x, x ) 

D = max depth[x) 

children(x) = {y: xy G E(G) and depth(y) = depth[x) + 1} 

2) For d = D, D - 1, . . . , 0, do: 

For all vertices x of depth d, set 

1 + W(x,y)w eS 



y£children(x) * 



[Note that all the needed Wy S have been set at the previous 
iteration. If Wy S = —1 for any child y of x, the algorithm is 
declared to fail.] 

3) Output 

Mw) = n ( i+ < ff )- ( 3 - 5 °) 

xdV{G) 

The correctness of this algorithm (when it succeeds) follows from (j3.44bj) . (|3.44cjl and 
f|3.48|) . Indeed, for any vertex x, the w^f produced by this algorithm is w^ ,Gx where 
G x is the subtree consisting of x and all its descendants. 

This also provides an algorithm for testing whether a given vector R lies in 1Z(W): 

Theorem 3.2 Consider any repulsive lattice gas with hard-core self-repulsion for 
which the support graph Gw is a tree. Then a vector R > lies in TZ(W) if and 
only if Algorithm T with w = — R produces —1 < w^r < for all vertices x. 



Proof. If Algorithm T with w = — R produces — 1 < u^ ff < for all x, then by 
the monotonicity of (J3.49|) in each Wy S it also does so when w = — R' for any vector 
< R < R. By (13301) this means that Z w (-R!) > for all such R, hence that 
R G K{W). 
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Conversely, let x be a vertex such that w x s (—1,0] but Wy S G (—1,0] for all 
descendants y of x. Since Wy S G (—1,0] for all children y of x, and w x = — R x < 0, it 
follows from ((3.49)1 that w x s < 0, so we must have w x G < — 1. Now set = R z for z 
in the subtree consisting of x and its descendants, and R' z = otherwise. Applying 
Algorithm T to w' = -R', it follows from (13301 that Z w (w') < 0. Thus R' G" K(W); 
and since TZ(W) is a down-set, we have R G" 7£(W). I 

Example 3.6 revisited. Consider again the tree Tn \ Comparison of ()3.38|) with 
(13.49(1 shows that in the multivariate case w x s = Y n (w) — 1 for all vertices x of height 
n above the leaves. Then ()3.39() is equivalent to ()3.50j) . 

3.5 Upper bounds on 1Z{W) when Gw is a tree 

We can use Theorem 13.21 to prove upper bounds on the set 1Z{W) whenever the 
support graph Gw is a tree. Let us begin by considering the special case of a hard-core 
pair interaction, i.e. the independent-set polynomial for a tree G on the vertex set X. 
As before, pick (arbitrarily) a root x G X, and define Aj = {x G X: depth(x) = i} 
for % = 0, . . . , D. Given any vector R > 0, define Ri to be the geometric mean of R x 
over all vertices x of depth i: 




(3.51) 



Now apply Algorithm T to w = — R, let p x 



—w x} and define 
\ V\Xi\ 




(3.52) 



(we assume here that w x < for all x G Xj). By ((3.49(1 we have 



n (i-jv)' 



(3.53) 



so that 



Ri 



(3.54) 





II (!-^) ^ l-T^lEft 



(3.55a) 





(3.55b) 
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by two applications of the arithmetic-geometric mean inequality. It follows that 

^ - • (3 - 56) 

So, suppose we fix the numbers {Ri)f =Q and define (j>i)fLo by the recursion 

R- 

with initial condition Pd+i = (which is well-defined as long as Pi+i remains < 1). 
Then it follows immediately from the monotonicity of the right-hand side of (|3.56|) 
that pi > pi as long as the former is well-defined. In particular, if for some level % we 
have Pi > 1, then pi is either > 1 or else ill-defined, so that there exists x G Xi with 
wl s ^ (—1,0]. It then follows from Theorem 13.21 that R ^ 1Z{W). Moreover, since 
this calculation depends only on the numbers (Ri)fL , the same conclusion holds for 
any R with the given geometric means on levels. 

Let us now generalize this argument to the case of a soft-core pair interaction 
< W(x, y) < 1. It is easily verified that 1 — Wp > (1 — p) w whenever < p < 1 
and < W < 1. Therefore, defining p x = —w^ as before and applying Algorithm T 
to w = — R, we have 

Px - n a-Py) 1 -^ ' (3,58) 

y£children(x) 

Now let us assign to each vertex x a weight 

d 

a x = Hil-Wix^Xj)} (3.59) 
j'=i 

where xo,x±, . . . , xj, = x is the unique path in Gw connecting x to the root xq. [In 
particular, a xo = 1 and if y is a child of x, then a y = [1 — 2/)]a x .] Setting 



aj = 



5^ « x , (3.60) 



let us define Ri and pj to be the weighted geometric means 



(\ 1 / a i 
n Ra * ( 3 - ei ) 

n ( 3 - 62 ) 
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The preceding argument can then be repeated verbatim (using the weighted arithmetic- 
geometric mean inequality), yielding 

fi - (i-jww* ' (3 ' 63) 

which can be analyzed analogously to ()3.56j) . We have therefore proven: 

Proposition 3.3 Consider any repulsive lattice gas with hard-core self-repulsion for 
which the support graph Gw is a tree. Define weights (a x )xex o^nd (a.i)f =Q by / K3.6ty) . 
and define numbers {j>i)f =0 by the recursion 

Pi = T\ ~ R \a^la ( 3 - 64 ) 

(as long as % + \ remains < 1), with initial condition Pd+i = 0. Suppose that there 
exists an i for which % > 1. Then R ^ 1Z(W) for all vectors R > satisfying 

]1 Rx x/ ^ > Ri for all i. 

xeXi 

Just as Proposition 13.31 "homogenizes" each level of a tree, we can also "homoge- 
nize" between levels. Let us again begin by considering the special case of a hard-core 
pair interaction. We first define branching factors hi = \X i+ i\/\Xi\. Choose any 
number b > and define 

7* = \Xi\/V . (3.65) 
Fix an integer k > 1 (the number of levels to be averaged together) and define 

i+k-l 

li,k = ^2 Ij ■ (3-66) 

j=i 

Then define the weighted inter-level geometric means 

i+fc-l \ VTi.fc 



Vi.k 



I [ f/ (3-67) 



Ri,k = ( n ^7 ( 3 - 68 ) 



for % — 0, . . . , D + 1 — k, and pi t k = R^k = for % > D + 2 — k. It follows from (|3.56|) 
that 

~ . Ri,k 
PiJk > 



n (i-p j+ i)W7i,* 
j=i 



^ . (3.69) 
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Now 

i+k-l i+k-l 

n (i-p j+ iv j+i/ii+i ' k < i - e — Pi+i 

AA z — ' 7i+i,fc 

i+fc-1 

fTj+i/7i+i,fc 



s i - n pjs 



= 1 - % +1 (3.70) 
by two applications of the arithmetic-geometric mean inequality, so that 

R- 

hk ~ (T^^w^' (3 ' 71) 

for i — 0, . . . , D + 1 — k. This holds true for an arbitrary choice of the numbers b and 
k. 

In the special case in which the sequence {b , bi, . . . , &d-i} is periodic of period k 
(where D > k) , we can choose b to be the geometric mean of the bi over one period: 

'jfe-i \V* 

II M • (3-72) 

vi=0 / 

Then the sequence {70, 71, ... , 7d} is also periodic of period k, so that 70^ = 7i,& = 
. . . = 7D + i_fc,fc, and ()3.71|) simplifies to 



R 



Pi,k > 7" ^ 7r (3.73) 



(1 -p 4 +l,fc) b 

for i = 0, . . . , D + 1 — k [note that for i = D + 1 — k we have Pi+i,* = 0, so the value 
of the exponent &7i+i,fc/7i,fc is in this case irrelevant]. 

So we can argue as before: suppose we fix the numbers (Ri,k)i^o l ~ k and define 
(Pi,k)i^Q l ~ k by the recursion 

^ = (l-^)^W7 !lfc (3 ' 74) 

(as long as Pi+\ y k remains < 1), with initial condition PD+2-k,k = 0. Then it follows 
immediately from the monotonicity of the right-hand side of (|3.71j) that Pi t k > Pi t k as 
long as the former is well-defined. In particular, if for some level i we have %^ > 1, 
then is either > 1 or else ill-defined, so that there exists x E Ujt^" 1 -^0 with 
w^ ff ^ (—1,0]. It then follows from Theorem 13.21 that R ^ 1Z{W). Moreover, since 
this calculation depends only on the numbers {Ri,k)f=Q~ h \ the same conclusion holds 
for any R with the given geometric means on-and-between levels. 

The same argument works with soft-core pair interaction if we replace \Xj\ by 5j. 
We have therefore proven: 
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Proposition 3.4 Consider any repulsive lattice gas with hard-core self-repulsion for 

which the support graph G\y is a tree. Define weights (a x ) x ex and ifiii)f =Q by 

Fix an integer k > 1 and a real number b > 0, and define 7, = oti/b 1 . Then define 

{li,k)?=o~ k , (Pi,k)F=o 1 ~ k and (R l , k )^=o 1 ~ k h WM> ~ WM> - Finally, define numbers 
(pi,k)iLo by the recursion 

(as long as remains < 1), with initial condition pc+2-k,k = 0. Suppose that 

there exists an i for which %^ > 1. Then R ^ TZ(W) for all vectors R > satisfying 
i+k-i ( _ \ 7i / 7i ' fe 

II II R* x/ " 3 I > Ri,k for alli = 0,...,D + l-k 

j=i 

We will return to these ideas in Section 



4 The Lovasz local lemma 

4.1 Hard-core version 

Let (A x ) x( zx be a finite family of events on some probability space, and let G 
be a graph with vertex set X. We say that G is a dependency graph for the family 
(A x ) x( zx if, for each the event A x is independent from the cr-algebra c(A y : y G 

X \ T*(a;)). [Here we have used the notation = U {x}, where T(x) is the 

set of vertices of G adjacent to x.] Note that this is much stronger than requiring 
merely that A x be independent of each such A y separately. 

A family of events typically has many possible dependency graphs: for instance, 
if G is a dependency graph for events (A x ) x< zx, then any graph obtained by adding 
edges to G is also a dependency graph. In particular, if the events A x are independent, 
then any graph on X is a dependency graph. Nor must there be a unique minimal 
dependency graph. Consider, for instance, the set of binary strings of length n with 
odd digit sum (giving each such string equal probability), and let Ai be the event 
that the ith digit is 1. Any graph without isolated vertices is a dependency graph for 
this collection of events. 

There is also a stronger notion of a dependency graph G for a collection of events 
(A x ) x( zx, where we demand that if Y and Z are disjoint subsets of X such that G 
contains no edges between Y and Z, then the a-algebras o~(A y : y gF) and u(A z : z G 
Z) are independent. In this case we shall refer to G as a strong dependency graph for 
the events (A x ) x£ x- 13 Alternatively, the dependency-graph hypothesis can be replaced 
by a weaker hypothesis concerning conditional probabilities, as in the lopsided Lovasz 

13 For instance, this situation arises in any statistical-mechanical model with variables living on 
the set X and pair interactions only on the edges of G, where each A x depends only on the variable 
at x. 
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local lemma (Theorem I1.2JI . It will follow from Theorem 14.11 below that all three 
hypotheses lead to the same lower bounds on B^D^ex Ar)- 

The various forms of the Lovasz local lemma (e.g. Theorems II. II and 11.2)) provide 
a sufficient (but not necessary) condition to have PdH^ex Ar) > ®> ^ n terms of the 
existence of numbers (r x ) x ex satisfying p x < r x Y[ ye j<r x ){l — r y ) for suitable proba- 
bilities or conditional probabilities p = (p x ) x ex- As discussed in the Introduction, 
we shall approach the problem by dividing our analysis into two parts. First, in this 
section, we examine a best-possible condition to have P(flxex Ac) > ^> m terms of the 
independent- set polynomial Zg{— p); then, in the next section, we discuss a sufficient 
condition to have Zg{— p) > 0, in terms of the existence of such numbers (r x ) xe x or 
generalizations thereof. 

The following result is a development of Shearer |9H1 Theorem 1]. 

Theorem 4.1 Let (A x ) xe x be a family of events on some probability space, and let 
G be a graph with vertex set X . Suppose that (p x )xex o,re real numbers in [0, 1] such 
that, for each x and each Y C X \ Y*(x), we have 



(a) Ifpe TZ{G), then 
and more generally 



V{A X \ p| A y ) < p x . (4.1) 

yeY 

F(f]A x ) > Z G (-p) > (4.2) 



x£X 



for any subsets Y,Z C X. Moreover, this lower bound is best possible in the 
sense that there exists a probability space on which there can be constructed a 
family of events (B x ) xeX with probabilities F(B X ) = p x and strong dependency 
graph G, such that W(0 xeX B x ) = Zc(—p)- 

(b) If p TZ(G), then there exists a probability space on which there can be con- 
structed: 

(i) A family of events (B x ) xe x with probabilities F(B X ) = p x and strong de- 
pendency graph G, satisfying P(flxex ^x) — 0/ and 

(ii) A family of events (B x ) xe x with probabilities f(B' x ) = p' x < p x and strong 
dependency graph G, satisfying F(B' X fl B') = for all xy G E{G) and 

p(a e x^i)=o. 

Remarks. 1. Please note that G is here an arbitrary graph with vertex set X; it 
need not be a dependency graph for the events (A x ) xe x- Rather, given G, we can 
regard p as defined by 

P * = rc^(/^n^) (4-4) 

yeY 
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(this is clearly the minimal choice). There is then a tradeoff in the choice of G: adding 
more edges reduces p x (since there are fewer conditional probabilities to control) but 
also shrinks the set 7Z{G) (by the last sentence of Proposition I2.15|) . 

2. Though (|4.1j) is the weak hypothesis of the lopsided Lovasz local lemma (Theo- 
rem ll.2|) . we will prove in (a) and (b) that the extremal families (B x ) xe x and (B x ) xe x 
have G as a strong dependency graph. Therefore, all three dependency hypotheses 
lead to the same optimal lower bound on W(f} xe x ^)- 

3. The proofs given here of Theorems 14.11 and 14.21 are logically independent of 
nearly all of Theorem 12. 101 More precisely, if we define 1Z(G) and 1Z{W) by condition 
(b) of Theorem 12.101 then the only part of Theorem 12.101 that is used in the proofs 
of Theorems 14. II and 14.21 is the (relatively easy) implication (b) ==>- (f). 

Proof. For p G 71(G), we wish to define a family of events (B x ) x( zx [on a new prob- 
ability space] such that the hypotheses of the theorem are satisfied and W(f) xeX Bx) 
is as small as possible. An intuitively reasonable way to do this is to make the events 
B x as disjoint as possible, consistent with the condition (|4.1|) [or with either of the 
two stronger notions of dependency graph]. With this in mind, for A C X let us 
define 

Y[ Px if A is independent in G 
*eA ' (4.5) 

xeA ^ otherwise 

This defines a signed measure on the cr-algebra generated by (B x ) xe x] indeed, inclusion- 
exclusion gives 

F(f]B x nf]B x ) = ^(-l^l-Wptfl^) (4.6a) 

a;eA x^A IDA xel 

e (-i) mHAi n^ ( 4 - 6b ) 

IDA, I independent xEl 

= (-l)l A 'Z G (-p;A), (4.6c) 

where Zq(— p; A) is defined as in (|2.H7j) . In particular, taking A = 0, we have 
^(flxex Ba) — Zg{~ p)- Theorem 12.1 Of f) implies that ()4.6c|) is nonnegative for all A, 
so that (|4.5jl defines a probability measure on a(B x : x G X). [This is the probability 
measure defined in Theorem 12. lOf g) .] 

If Y and Z are disjoint subsets of X such that G contains no edges between Y 
and Z, it follows from (|4.5|) that for Y C Y and Z C Z the events flxeyo ^ x an< ^ 
HigZo ^ x are ^dependent. This implies (see, for instance, |116| Theorem 4.2] or [TBI 
Theorem 4.2]) that u(B x : x G Y) and u(B x : x G Z) are independent, and so G is a 
strong dependency graph. 

We next show that (B x ) x& x is a family minimizing P^xex B x ). For A C X, we 
define 

P A = P(f|A E ) (4.7) 

zeA 

QA = F{f\B x ). (4.8) 

x€A 
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Let us now prove by induction on |A| that Pa/Qa is monotone increasing in A. Note 
first that by inclusion-exclusion, 

Q A = ^(-1)1^(1^) (4.9a) 

I(ZA x&I 

e (-i) i7i n^ ( 4 - 9b ) 

/CA, / independent xdl 

= Z G (-pl A ). (4.9c) 
Thus Qa > for all A, since p G 71(G) and 71(G) is a down-set. Furthermore, for 

i"CAU{j/}, / independent xEl 

= qa-vv £ (-i) |7i n^ ( 4 - iob ) 

7CA\r(j/), I independent x£l 
= QA ~ PyQ.A\T{y) • (4.10c) 

[Note that this is just the fundamental identity (|3.5|) applied to Z G (—p 1a)-] On the 
other hand, 

Pau{v} = Pa - W(A y n P| A x ) (4.11a) 

xGA 

> Pa- P(A,n f| A,) (4.11b) 

xeA\r(y) 

> P A - PyP A \r(y) (4.11c) 

by the hypothesis (|4.1|) . Now we want to show that P A u{y}/QAu{y} > Pa/Qa, or 
equivalently that -Pau{j/}Qa — <5au{j/}-Pa > 0. By ()4.1(J|) and (|4.11|) we have 

Pau{ v }Qa- Qau{ v }P\ > [Pa -P y PA\T{y)]QA - [Qa - P y Q a\t( v )]Pa (4.12a) 

= Py [PAQA\T(y) ~ QaPa\T(ij)] (4.12b) 

> (4.12c) 

since 

ft > £™. (4.13) 

VA VA\r(2/) 

by the inductive hypothesis. 

Since Pa/Qa is monotone increasing in A, we have Px/Qx > P0/Q0 = 1, which 
proves ()4.2j) . More generally, for any subsets Y, Z C X, we have Pyuz/Qyuz > 
Pz/Qz and hence Pyuz/Pz > Qyuz/Qz, which gives (jOji . 

For p ^ 7Z(G), choose a minimal vector p' < p such that p' > and Zq(— p') = 
[such a p' is in general nonunique]. Then the family of events (B' x ) x& x defined by 
fl4~5jl with p x replaced by p' x satisfies F(f) x&x B~' x ) = Z G (-p') = [by (EM with 
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A = 0}. Since p' is in the closure of 7Z(G), it follows by the minimality of p' and the 
continuity of Zq that this is a well-defined probability measure; note that if x and y 
are adjacent then F(B' X fl B' ) = by (|4.5|) . Thus we have constructed a collection of 
events satisfying part (ii) of the Theorem. 

To construct a collection of events satisfying part (i), let (C x ) xe x be an (indepen- 
dent) collection of independent events satisfying 

[1 - ¥(B' X )} [1 - P(Cy] = 1- Px . (4.14) 
Then the events B x = B' X U C x satisfy ¥(B X ) = p x and f(f]B x ) < P(f]X) = °- ■ 

Remarks. 1. If (A, x ) x& x is a family of events satisfying ()4.2j) with equality, then we 
have Px = Qx in the foregoing proof; and since P = Q = 1, the monotonicity 
of Pa/Qa implies that we have Pa = Qa for every A C X. Thus, if (A x ) xGX is 
an extremal family, the probabilities of all events in a(A x : x G X) are completely 
determined and are given by ()4.5j) / ()4.6j) . 

2. The Lovasz local lemma can be formulated more generally for families of 
events with a dependency digraph: each event A x is independent from the u-algebra 
cr(A y : y G X\r^(x)), where F*,(x) = T + (x)U{x} and T + (x) is the out-neighborhood 
of x. See e.g. Lemma 5.1.1] or |2*3*1 Theorem 1.17]. It would be interesting to have 
a digraph analogue of Theorem 14. 1[ but we do not know how to do this. 

3. There are other probabilistic inequalities that are expressed in terms of a 
dependency graph (see for instance Suen |lU9j or Janson 568 ); it would be interesting 
to know if any of these have counterparts in statistical mechanics. One obstacle 
here is the need for a counterpart of Theorem 14.11 However, even without such a 
result, there may be scope for proving further inequalities in the presence of weak 
dependency conditions of the form discussed in the next subsection. 

4.2 Soft-core version 

Let us now consider how to extend Theorem 14.11 to the more general case of a 
soft-core pair interaction, i.e. to allow "soft edges" xy of strength 1 — W(x, y) G [0, 1]. 
The first step here is to replace the hard-core dependency condition (|4.1|) by an 
appropriate soft-core version. 

Let W: X x X — > [0, 1] be symmetric and satisfy W(x, x) = for all x G X; and 
let (A x ) x& x be a collection of events in some probability space. For each x G X, let 
S x be a random subset of X, independent of the a-algebra a(A x : x G X), defined by 
the probabilities 

F(yeS x ) = W(x,y) (4.15) 

independently for each y G X. [Thus in the case of a hard-core pair interaction, we 
have S x = X \ T*(x) with probability 1.] Let {p x ) x( zx be real numbers in [0, 1]. We 
say that (A x ) x& x satisfies the weak dependency conditions with interaction W and 
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probabilities (p x ) X £x if, for each i6l and each Y C X \ x we have 

EMP^n f) A y U < Px e(f( p| I,)) . (4.16) 

V j/eYns, / \ y eYnS x J 

[Note that in the special case of a hard-core pair interaction, we have YC\S X = Y\T*(x) 
with probability 1, so that (|4.16|) reduces to (|4.1|) .] Of course, the reference here to 
a random subset S x can be replaced by an explicit expression for the probabilities 
P(y n S x = Y'), so that KWh is equivalent to 



E 11^^) ( II [l-W(x,y)}\F(A x n P) Ay) < 

Y'CY \y£Y> / \y£Y\Y' J y&> 

Px J2 ( II W ( x > y) J ( II l±-W(x,y)}\nf)A y ). (4.17) 

Y'CY \<,eY' / \?/eY\Y' / y&' 

We can now state a soft-core version of Theorem 14.11 

Theorem 4.2 Let (A x ) xe x be a family of events in some probability space, and let 
W: X x X — > [0, 1] be symmetric and satisfy W(x,x) = for all x G X . Suppose 
that (A x ) x& x satisfies the weak dependency conditions \4-lb^ / \4-ll\) with interaction 
W and probabilities {p x ) x ^x- 

(a) Ifpe 1Z(W), then 

¥(f]A x ) > Z w (-p) > (4.18) 

and more generally 

nr\A x \r\A x ) > ^ ( 7 plyu f } > o (4.19) 

1 ' 1 ' " Z w (-pl z ) K } 

for any subsets Y, Z C X . Furthermore, this bound is best possible in the sense 
that there exists a family {B x ) x& x with probabilities F(B X ) = p x that satisfies the 
weak dependency conditions J^.i6| )/ J^.i71 ) with interaction W and probabilities 
{Px) X £X, has strong dependency graph Gw, and has F(f] x€X ^x) — Zw(~P)- 

(b) If p ^ 1Z{W), then there exists a probability space on which there can be con- 
structed: 

(i) A family of events (B x ) x£ x with probabilities F(B X ) = p x and satisfying the 
weak dependency conditions d^.i6| )/ d^.i7| ) with interaction W , such that 
nf] x& xB x )=0; and 

(ii) A family of events (B x ) xe x with probabilities F(B' X ) = p' x < p x and satisfy- 
ing the weak dependency conditions \4 ■ 1 b\j / l4 ■ 1 7j ) with interaction W , such 
that F(B' X n B' y ) = W(x,y)¥(B' x )¥(B' y ) for all x,y and W>(f\ xeX ^ x ) = °- 
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Proof. For p e 1Z(W), we define a family of events (B x ) xe x [on a new probability 
space] by 

p (fl^) = [Up-] ( II w ^y)\ ( 4 - 2 °) 

xeA \x£A / \{x,y}CA J 

[here the second product runs over all two-element subsets {x,y} C X (x ^ y)). As 
before, inclusion-exclusion gives 

¥(f]B x nf]B x ) = ^(-^l^-^PCfl^) (4.21a) 

= E(- i ) |Ai ri(-^) n w ^y) ( 4 - 2ib ) 

IDA x£l {x,y}CI 

= (-l)l A IZ w (-p;A). (4.21c) 

In particular, we have F(f] xeX B x ) = Z w (—p) > 0, and more generally for A C X 
we have 

xeA ica xei 

= (Up*) f II w ^v) 

ICA \x£l J \{x,y}CA 

= Z w (-pl A ) 
> 0. 

We define P A and Q A as in IftTj) / (j4~K|l . Note that Q A = Z w (-pl A ) > 0, as 
= Z w (—p) > and ^(W 7 ) is a down-set. Let us now prove by induction on |A| 
that P\/Qa is monotone increasing in A. For A C X and y G X \ A, we have 

QauM = Z w (-p 1au{^}) (4.22a) 
= Z w (-pl A ) - Py Z w (-W(y, -)pl A ) (4.22b) 

= ZwrC-pU) -PvJ2 \Il w ( x >v)) ( II [i-^(^y)] ]Z w (-pl Y ) 

YcA \ y eY ) \j/eA\y / 

(4.22c) 

= qa - p y ( n ^ y))[ n t 1 - ^ I ^ ( 4 - 22d ) 

FcA \yeY J \yeA\Y J 

where we have used the fundamental identity (j3.3|) / ()3.6|) . On the other hand, 

PauM = PA-P^nf)^) (4.23a) 
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PA-J2(U w ^yn ( II [i-w(x, y )}\nA y nf]A x ) 

YCA \x&Y J \xeA\Y J x£A 



(4.23b) 



^ P ^~J2 [H W ( X >V)) ( II [l-W(x,y)]\F(A y nf)A x ) 

YCA \x£Y J \xeA\Y J x&Y 



(4.23c) 



> P A -PyJ2 \H W &y)) ( II [1-^,2/)] I PCfl^). 

YCA \x6Y / \zGA\Y / x&Y 

(4.23d) 

where the last line uses the weak dependency condition 1)4.16)1 / (|4.17|) . We then have 



Wa 



p a-p,e f n w ^y)) ( n [i-w(x, y )] ] 

YCA Wy / \ z6A\Y / ^ A 



(4.24a) 

< Pa-PvZ) ( n^'^J ( II t 1 - W &V)\ I P ^ ( 4 - 24b ) 

YCA VxeY / \z6A\Y / 

< ^auM, (4.24c) 

where the first inequality uses the inductive hypothesis and the second inequality uses 
(14.23)) . Hence PauM / 'Q 'au{ v } > Pa/Qa as claimed. 

The bounds (14.18)) and (J4.19)) follow as in Theorem 14.11 

Part (b) follows as in Theorem 14.11 For p ^ 1Z(W) choose a minimal vector p' 
with < p' < p and Z w (—p') = 0. The events (B' x ) xeX defined by (I4.2U)) with p x 
replaced by p' x satisfy F(f) xeX B x ) = and F(B' X n B' y ) = W(x, y)F(B' x )F(B' y ) for all 
x, y. A collection of events satisfying (i) is constructed as before. I 

Remark. The uniqueness of a family of events satisfying ()4.18j) with equality follows 
as before: since Px/Qx = P0/Q0 = 1 and Pa/Qa is increasing we have Pa = Q\ for 
all A C X and so the probabilities of events are given by ()4.21)) . 



5 Sufficient conditions for Zw 7^ in a polydisc 

In this section we shall exhibit sufficient conditions on a set of radii R = {R x }xex 
so that the partition function Z^y(w) is nonvanishing in the closed polydisc |w| < R. 
We shall restrict attention to the case of a repulsive lattice gas with hard-core self- 
repulsion, i.e. < W(x, y) < 1 for all x, y and W(x, x) = for all x, and we shall 
use the notation introduced in Section 13.11 Our main tool will be the fundamental 
identity ()3.3)) . applied inductively. 
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5.1 Basic bound 

Our first (and most basic) bound is due to Dobrushin [S3 EHj in the case of a 
hard-core interaction; the generalization to a soft repulsive interaction was proven 
recently by one of us |101j . The method of proof is, however, already implicit (in 
more powerful form) in Shearer [HH Theorem 2]. 

Theorem 5.1 (Dobrushin [371131], Sokal [TTTT]1 Let X be a finite set, andletW 
satisfy 

(a) < W(x, y) < 1 for all x,y e X 

(b) W{x,x) =0 for allx EX 

Let R = {R x }xex > 0. Suppose that there exist constants {K x } xe x satisfying < 
K x < 1/R X and 

l-W{x,y)K y R y 



k x > n 

y^x 



1 - KyRy 



(5.1) 



for all x G X . Then, for each subset A C X , Z\(w) is nonvanishing in the closed 
polydisc Dr = {w G C x : \w x \ < R x for all x} and satisfies there 



d\ogZ A (w) 



<9uu 



< 




for all x G A 



for all x G X \ A 

Moreover, ifw, w' G Z)r and w' x /w x G [0, +oo] for each x G A, then 



log 



Za(w' 



log 



1 - K x \w' 



1 - K x kw, 



(5.2) 



(5.3) 



where on the left-hand side we take the standard branch of the log, i.e. \ Imlog • 



< 7T. 



Remark. It follows from (j5.1j) that K x > 1 and hence that R x < 1. 

PROOF. Note first that (|5.2jl for any given A implies (|5.3j) for the same A, by inte- 
gration. 

The proof is by induction on the cardinality of A. If A = the claims are trivial. 
So let us assume that ()5.2j) [and hence also (J5.3J) ] holds for all sets of cardinality < n, 
and let a set A of cardinality n be given. Let x be any element of A. Let us apply 
the fundamental identity (|3.3jl . and observe that W(x, -)w G Dr since |VF(x, ?/)| < 1. 
Therefore, by the inductive hypothesis we have Z A \ x {w) ^ and Z\\ x (W(x, -)w) ^ 0; 
and from ()3.3|) we have 

d i v ( \ ^x,a(w) 

■logZ A (w) = — — - — -- — (5.4) 



dw x 1 + K, X}A (w)w x 
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where 

Z AV (W(x, -)w) 
^,a(w) = — x — — . (5.5) 

Za\x-(w) 

Now by the inductive hypothesis (|5.3|) for A\x, and using the fact that w' = W(x, -)w 
satisfies w' y /w y = W(x, y) > 0, we have 

| Ki , (w) i < n 1 ~ < n i - vy(»,g)jr si (5 6) 

j/gA\x y y y&X\x y y 

which is < K x by the hypothesis (|5.1|) . This proves (|5.2|1 for A, and hence completes 
the induction. ■ 

It is convenient to rewrite Theorem 15.11 in terms of the new variables r x = K X R X : 

Corollary 5.2 Let X be a finite set, and let W satisfy < W(x,y) < 1 for all 
x, y G X and W(x, x) = for allx G X . Suppose that there exist constants < r x < 1 
satisfying 

Rx < r x T[- )~ Ty (5.7) 

for all x & X. Then, for all w satisfying |w| < R, the partition function Zy/ satisfies 

\z w (v)\ > z W (-r) > n( i - r -) > ° ( 5 - g ) 

and more generally 



Z w (wl: 



x£Y 

In particular, if we define the "maximum weighted degree' 



and write 



A w = max^[l- W(x,y)\ (5.10) 

y+x 



F{Aw) = 2 + Aw - v^Tiy (5 n) 

R(A W ) = F(A w )e^- F ^ (5.12) 



we have 



\ZwM\ > [1-F(A W )] |X| >0 (5.13) 
whenever \w x \ < R(A\y) for all x G X . 
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Proof. Setting r x = K x R Xl we find that (J5~T]) becomes (|B~7|) . and (f5~3j) with A = X 
and w' = becomes ()5.8|) . 

To obtain the last claim, note first that 

' ' ' ' 1 > e -^W)r/(l-r) (514) 



1-l^r l-r + (l-W)r 1 + 

whenever < W < 1 and < r < 1. Therefore, if we set r x = r for all x G X, we 
have 

r * IT i u7, ^ \ ^ re-^liX-r) . (5 . 15) 

We then choose r to maximize the right-hand side of (j5.15|) ; simple calculus yields 
r = F(A W ) and A w r = (1 — r) 2 , so that the right-hand side of ()5.15|) is bounded 
below by R(A W ). It follows that if we define R x = R(A W ) and r x = F(A W ) for all 
x G X then (|5.7|) and so (|5.8|) are satisfied. I 

Remarks. 1. The radius R(Aw) behaves as 

1 - 2A^ 2 + \A W + 0(A% 2 ) as A w 
R( A w) = { , r , /A (5-16) 



i 



as Avk — > oo 



Example El 6 (the r-ary rooted tree) shows that this bound is sharp (to leading order) 

1/2 

as Aw — » oo. At the other extreme, the 1 — const x A w behavior at small Aw 
is also best possible, since the two-site lattice gas with W(x,x) = W(y,y) = 
and W(x,y) = 1 — e has Zw{w) = 1 + 2w + (1 — e)w 2 and hence has a root at 

1/2 ^ ^ 

w = — 1/(1 + y/e). [However, the coefficient 2 rather than 1 in the A^ term of (|5.1fij) 
may not be best possible.] 

2. From Proposition I2.15f c) we know that the set TZ(W) is log-convex; it is 
therefore natural to ask whether the subset of TZ(W) produced by ()5.7|) is also log- 
convex. (If it weren't, then we could improve Corollary 15.21 bv taking the log-convex 
hull.) It turns out that the set of vectors satisfying (J5.7|) is indeed always log-convex: 
in fact, if (r, R) and (r', R') are two pairs of vectors satisfying ()5.7jh and < A < 1, 
then (r A r /1_A , R A R /1_A ) also satisfies (|5.7|) . In the hard-core case this follows from the 
inequality 1 — r A {r' y ) 1 ~ x > (1 — r y ) x {\ — r' y ) 1 ~ x , which is proven by two applications 
of the weighted arithmetic-geometric mean inequality as in (|3.55|) . In the general 
case it can be shown by a similar argument (using four applications of the weighted 
arithmetic-geometric mean inequality!). 

Specializing Corollary 15.21 to the case of a hard-core pair interaction for a graph 

G, 

j j t / \ fO if x = y or xy G E(G) 1P _s 
W(x,y) = ., , , y , (5.17) 

v ' yj \1 itx^y and xy £ E(G) v ; 

we have: 
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Corollary 5.3 Let G be a finite graph with vertex set X , and let R = {R x }xex > 0. 
Suppose that there exist constants < r x < 1 satisfying 

R x < r x H (1 - r y ) (5.18) 

y£T(x) 

for all x G X . Then, for all w satisfying |w| < R ; the independent-set polynomial 
Zq satisfies 

\Z G (w)\ > Z G (-R) > Y[(l-r x ) > (5.19) 

x&X 

and more generally 

x& 

In particular, if G has maximum degree A, then \Zg(w)\ > [A/(A + 1)]'^' > 
whenever \w x \ < A A / (A + 1) A+1 for all x G X . 



Z G (w l Y uz) 
Z G (w l z ) 



Proof. The last claim is obtained by setting r x = 1/(A + 1) for all x G X. 
Remark. The radius A A /(A + 1) A+1 behaves for large A as 

A A J_ I" J_ 7 3 n(A~ 4 ] 

(A + 1) A + X ~ iA [ ~ 2A + 24A2 ~ 16A3 + { ' 



(5.21) 



which agrees with 1)5.16)1 to leading order in 1/ A but is slightly larger (hence better) 
at order 1/A 2 . 

Combining Corollary 15.31 with Theorem 14.1) we immediately obtain the lopsided 
Lovasz local lemma (Theorem 11.2)) . Combining Corollary 15.21 with Theorem 14.2) we 
obtain a "soft-core" version of the lopsided Lovasz local lemma: 

Theorem 5.4 Let (A x ) x( zx be a family of events in some probability space, and let 
W: X x X — > [0, 1] be symmetric and satisfy W(x, x) = for all x G X. Suppose 
that (A x ) x£ x satisfies the weak dependency conditions I[4.1b\ )/ fi4 -Ity with interaction 
W and probabilities (p x ) x ex- Suppose further that {r x ) x( zx are real numbers in [0, 1) 
satisfying 

p x < r x H (l-r„) . (5.22) 

2/ero) 

Then 

X £X X £X 

and more generally for sets Y, Z C X , we have 

P(f| A x \ f| A x ) > H (l-r x ) = 0. (5.24) 

x & x ez xeY\z 
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Defining the weighted degree Aw as in (|5.10p . we obtain the following: 



Lemma 5.5 Let (A x ) xe x satisfy the weak dependency conditions \4-lb\j / \4-ll\) with 
interaction W and probabilities (p x )xex- If Px < ^.^/(Aw + l) A w+i j or ever y 
xeX, then¥(f] x£X A x )>0. 

Proof. As in the proof of Corollary EHH set r x = r = 1/(A W + 1) for all x G X. 
Then check (15. 7j) : 

1 — r„ -n- , 1-W(x,y) 



L J-l-W(x,y)r v ~ \1 K 

y ^x \ y y ^ x 

< r(l — r] 



In the first inequality we have used the fact that 1 — W(x, y)r y < (1 — r y ) w ^ x ' v ' for 
0<W(x,y)<l. ■ 

It would be interesting to see applications of Theorem 15.41 and Lemma 15.51 



5.2 Improved bound 

Let us now attempt to improve the bound of Theorem 15 .11 Note, first of all, that 
we need not insist that the bound (|5.2j) hold with the same constant K x for all A 3 x; 
rather, we can use constants K x \ that depend on A. Inspection of the inductive 
argument shows that we can define the constants K x ^ e [0, +oo] as a function of the 
family {R x } by the recursion 



K 



x.A 



1 — W(x,y)K y A\xRy 

[ [ — ~ if K V) \\ x R y < 1 for all terms in the product 

y e A \ x *■ ~ KyA\xRy 

W(x,y) £ 1 



foo otherwise 

(5.26) 

[Note that K Xj { x y = 1 for all x, because the product ()5.26j) is empty] It follows imme- 
diately by induction that each K Xj \ is an increasing rational function of {Ry} ye A\x up 
to the first pole, and +oo thereafter. More precisely, suppose we define the rational 
functions ^^(R) by the recursion ()5.26|) without the restrictions K y ^\ x R y < 1; then 
it is easy to see that 

Kx A ( R ) = | K xA {R) if K xA {R') < oo for < R' < R (5 _ 2?) 
I +oo otherwise 

Moreover, it is easily proven by induction that all the partial derivatives of K Xy \ are 
nonnegative at R = 0; it then follows from Proposition 12 . 1 ll that the Taylor series of 
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K X) \(H) about converges throughout the region where the first case in ()5.27j) holds, 
and that all the partial derivatives of K xA = K xA are nonnegative there. Finally, it 
is obvious from the definition (j5.26|) ff. that K x> \ is an increasing function of A. 

Let us now define a graph G with vertex set V = {x G X: R x > 0} and edge set 
E = {x,y G V: W(x, y) ^ 1}; and for each A C X, let G\ be the subgraph of G 
induced by A n V. Then only the connected component of G\ containing x plays any 
role in the definition of K x ^\ that is, if Ga has several connected components with 
vertex sets Ai, . . . , A& and x G A«, then K Xj \ = K x ^. 

Let us now call a pair (x, A) "good" if K Xj a < oo and K xA R x < 1. It follows 
immediately from the definition ()5.26|) ff . that if (x, A) is good, then (y, A\x) is also 
good whenever y G A \ x with W(x, y) ^ 1 and R y > 0, i.e. whenever y is a neighbor 
of x in G\. [Indeed, this follows under the weaker hypothesis that K x ^ < oo.] 

We then have: 

Theorem 5.6 (improved Dobrushin— Shearer bound) Let X be a finite set, and 
let W satisfy 

(a) < W(x, y) < 1 for all x, y G X 

(b) W(x,x) = for allxeX 

Let R = {R x } x( zx > 0. Define the constants K Xs a G [0, +oo] as above. Suppose that 
in each connected component of G\ there exists at least one vertex x for which the 
pair (x, A) is good. Then Z\(w) is nonvanishing in the closed polydisc Dr.; and for 
every good pair (x, A) and every w G Dr,, we have 



<91ogZ A (w) 



dw 7 



< 



K 



x.A 



1 - K x a\w x 



(5.28) 



Moreover, if w,w' G -Dr and w' x /w x G [0, +oo] for each x G A, and in addition the 
pair (x, A) is good whenever w' x ^ w x , then 



log 



Za(W) 



Z A (w) 



< 



E 

x £ A 



log 



1 - K x a\w' x \ 



1 - K x a\w x \ 



where on the left-hand side we take the standard branch of the log, i.e. | Imlog • 



(5.29) 



< 7T. 



PROOF. Note first that (J5.28|) for any given A implies (J5.29)) for the same A, by 
integration. 

The proof is by induction on the cardinality of A. If A = the claims are trivial. 
So let us assume that (j5.28j) [and hence also (|5.29J) ] holds for all sets of cardinality 
< n satisfying the stated hypotheses, and let a set A of cardinality n satisfying these 
hypotheses be given. If G\ is not connected, i.e. has components with vertex sets 
Ai, . . . , Afc (k > 2), then Z\(w) factorizes as 



(w) = \{Z 



(5.30) 



i=l 
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moreover, each connected component has cardinality < n and satisfies the hypotheses 
of the theorem, so it follows immediately from the inductive hypothesis that (|5.28j) 
and (|5.29|) hold also for A. We may therefore assume that G\ is connected. 

Let x be any element of A for which the pair (x, A) is good. Since G\ is connected, 
the vertex x has at least one neighbor y in each connected component of G\\ x . More- 
over, as noted previously, the pair (y, A\x) is good whenever y is a neighbor of x 
in G\. Therefore, the inductive hypothesis is applicable to A \ x. Let us now apply 
the fundamental identity (|3.3|) . and observe that W(x, -)w G Z)r since ?/)| < 1. 

By the inductive hypothesis we have Z A \ x (w) 7^ and Z A \ x (W(x, -)w) 7^ 0; and from 
(J3.3|) we have 

logZ A (w) = — — — — — (5.31) 



dw x 1 + K, x> a(w)w x 

where 

Z A \ x (W(x, » 
£x,a(w) = — K — — — . (5.32) 

Now each y 6 A \ x with w' y 7^ w y necessarily has W(x,y) 7^ 1 and R y > (i.e. 
is a neighbor of x in G\), so that the pair (y,A \ x) is good. We may therefore 
apply the inductive hypothesis ()5.29|) : using the fact that w' = W(x, -)w satisfies 
w ' y / w y = W(x,y) > 0, we have 

it, A( w)i < n 1 - w(*<v)><.^ < KxA (5JB) 



y e A \ x 
W(x,y) ? 1 



since \w y \ < R y . This proves ()5.28|) for (x, A), and hence completes the induction. 



Remark. It is unclear whether the set of vectors R satisfying the hypotheses of 
Theorem 15.61 is log-convex. If it is not, then the conclusion of Theorem 15.61 can be 
improved by taking the log-convex hull. 

As a corollary of Theorem 15 .61 we can deduce a bound due originally (in the Lovasz 
context) to Shearer [§H1 Theorem 2], which improves the last sentence of Corollary 15 .31 
by replacing A by A — 1. Indeed, we can very slightly improve Shearer's bound by 
allowing one vertex xq to have a larger radius R xo : 

Corollary 5.7 Let G = (X, E) be a finite graph of maximum degree A > 2, and fix 
onevertexxo EX. Suppose that \w XQ \ < (A — 1) A /A A and that \w x \ < (A — 1) A ~ 1 /A A 
for all x 7^ x . Then Z G (w) 7^ 0. 



Proof. Since Zq factorizes over connected components, we can assume without loss 
of generality that G is connected. (Indeed, if G is disconnected, then we can allow 
one "xo-like" vertex in each connected component.) Set R xo = (A — 1) A /A A and 
R x = (A - 1) A ~VA A for all x ^ x . 
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We first claim that if xq ^ A, and x G A is a vertex with at least one neighbor in 
X\A, then 



A ^ 



K xA < [-^lj (5 ' 34) 

(note the strict inequality). The proof is by induction on |A|, using the definition 
(|5.26|) : it certainly holds if A = {x}; since every y appearing in the product on the 
right-hand side of ()5.26j) has at least one neighbor outside of A\a; (namely, x itself), 
Ky,A\x satisfies (|5.34j) by the inductive hypothesis and so K y \\ x R y < 1/A; and finally, 
since x has at least one neighbor outside A, there are at most A — 1 factors in the 
product. Thus 



K X)A < 
It then follows that 



A-l / A xA-1 



1/AJ \A 



K xo ,x < [ - K — [ ) , (5.35) 



since the bound (|5.34|) applies to all the terms K y ^\ xo appearing on the right-hand 
side of (|5.26p . We therefore have K XOtX R xo < 1, and so the pair (x ,X) is good. The 
claim then follows from Theorem 15.61 ■ 

Replacing A A / (A + 1) A+1 by (A— 1) A ~ 1 /A A may seem to be a negligible improve- 
ment, since both quantities have the same leading behavior « l/(eA) as A — > oo, 
and differ only at higher order: 



(A-1) A - X 1 

A^ ~ 7K 



1 7 



2A 24A 2 16A 3 



0(A- 4 ) 



(5.36) 



[cf. (j5.21|) ]. 14 But Shearer's bound (A — 1) A ~ 1 /A A has the great merit of being best 
possible: for, as he showed j^Hl, if G is the complete rooted tree with branching 
factor r = A — 1 and depth n, then Zg{w) has negative real zeros that tend to 
w = —(A — 1) A_1 /A A as n — > oo (see Example|Hl6 above). 

We remark that Corollary 15 . 71 does not appear to extend naturally to the soft-core 
case (note that having one neighbor outside A in the argument around ()5.34|) need 
not reduce the weighted degree of a vertex in A by 1). 



5.3 Optimal bound 

Now let us try to further improve Theorem 15.61 Where in the proof did we lose 
equality? This happened only in ()5.33|) . where we applied (|5.29|h which in turn arose 
from integrating ()5.28|) : the point is that the bound ()5.28|) in general improves as 
we pass from w downwards to w' = W(x, -)w, but we failed to take advantage of 
this fact. The solution is to order (arbitrarily) the vertices y\, . . . ,yj. arising in the 

14 The amusing similarity of lj!T2T)) and (|53(3|) arises from the fact that -(-A)" A /(-A + l)" A+1 = 

(A-i) A -yA A 
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product ()5.33|) and to write /C Xj a(w) as a telescoping product; then we will have a 
sharp bound, i.e. one that becomes equality when w = — R. 

Let us show this first in the special case of a hard-core pair interaction: setting 
n A = {yi, y k }, we have 

r f m \ Z A\r*(x)(w) T-r ^A\ g \{yi,..., y< }(w) 
Za\ x M f = \ ^A\x\{ W i I ...,i/ i _ 1 }(w) 

Therefore we can improve ()5.33|) by writing 

k 

|/C„ A (w)| < \[ 7~ 7r- > ( 5 - 38 ) 

i=l 1 A wM{9l.->!/i-l}% 

and so we can replace the definition (|5.2(jj) by 

K °S = A rr^t 1 ^ (5-39) 

»=i j/»,A\x\{j/i,...,i/i_i} 

if -^^AVaiNfyi y < ^ ^ or a ^ terms i n the product, and = +°° otherwise 

(note that K°^^ = 1). Moreover, if we use K° p ^, then all the bounds in the proof 
become equality when w = — R. 

In the general case of a "soft" interaction W(x,y), things become slightly more 
complicated, since the vectors w' = W(x, -)w have their components depressed but 
not set to zero. To handle this case, we need to define the numbers K° P A explicitly as 
functions of a vector R^. We begin by writing JC X) \(w) as a telescoping product as 
in (pTTOjl : 

^,a(w) = n / M ;it-m ( 5 - 4 °) 

where the vectors w^' are defined by 



(w w ) = <! W(k% ' V ^ Wy l{y = Vj f ° r S ° me 3 ~ i (5 41) 

' w y otherwise 



Therefore we can improve (J5.33|) by writing 

\ K ^)\ < 11 — ^ny- (5-42) 

[where R^ is defined by the obvious analogue of (|5.41|) ]. and so we could have 
replaced the definition (|5.26|) by 
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if Ky^^iRytx )Ryi < 1 f° r all terms in the product, and K° P l = +00 otherwise. 

Another way of looking at all this is: We want to choose the constants K XjJ \ so 
that the bound ()5.28j) holds. By ()3.8|) and (j3.9J) . the optimal choice is manifestly 

KTk = sup £ XiA (-R') , (5.44) 

0<R'<R 

where /C Xi a(w) was defined in (|3.7p . We now claim that ()5.44|) is identical to ()5.43|) 
et seq. Indeed, the recursion (I3.13|) shows that the rational function /C Xj a(— R) is 
identical to what one obtains from the recursion 1)5.43)1 if one omits the condition 
following the equation; and the monotonicity of (j5.43j) up to the first pole guarantees 
that implementing the condition following the equation is equivalent to taking the 
supremum over all R' satisfying < R' < R. 

In summary, this second improvement of Theorem 15. II is optimal in the sense that 
it is equivalent to calculating the exact lC x ^{w) and hence [by (|3.16|) ] the exact Za(w). 
So this "optimally improved Dobrushin-Shearer theorem" is not usually going to be 
useful in practice; but it does give insight into what has been lost in Theorems 15.11 
and 15.61 

6 Tree interpretation 

Further insight into Theorem 15 . 61 and its "optimal" improvement a la ()5.43|) / ()5.44|) 
can be obtained by considering the tree structure underlying the recursions ()5.26p and 

6.1 Tree interpretation of Theorem 15.61 

Let us first consider (j5.26J) . which provides a recursive definition of the functions 
K Xj a(R) [or K X) \(TX) if we ignore the conditions K y ^\ x R y < 1 following ()5.26|) ]. This 
recursion can be encoded as a tree: at the top (root) of the tree is the pair (x,A); 
immediately underneath it are all the pairs (y, A \ x) for which W(x, y) 7^ 1 [that is, 
for which xy is an edge of the support graph G = Gw restricted to A]; and so on. 
Another way of saying this is as follows: 

Definition 6.1 (see also [78J) Let G be a simple loopless graph, and letx\ G V(G). 
Then the path-tree or self-avoiding-walk tree (SAW-tree) of G rooted at x± — let us 
call it SAW(G,xi) — is the graph whose vertex set is the set of paths in G starting 
at X\ [i.e. paths x\Xi ■ ■ ■ Xk in G with k > 1] and whose edges connect P to P' 
whenever P' is a one-step extension of P [i.e. whenever P' = x± ■ ■ ■ Xk with k > 2 
and P = Xi ■ ■ ■ Xk-i]. The root of SAW(G, x±) is the zero-step path x~[. [Note that if 
G happens to be a tree, then SAW(G, X\) is isomorphic to G for any x\ e V(G).] 

We then associate a path xi ■ ■ ■ Xk with the pair (x fc , A) where A = V(G)\{xi, . . . , Xk-i}- 
Note that a given pair (xk,A) can correspond to many paths xl ■ ■ ■ Xk, which arise 
from different orderings of the set {X2, ■ ■ ■ , Xk-i} = V(G) \A\x\. 



66 



Given a nonnegative vector R = {R x }x<=v(G), let us now assign to the vertices of 
SAW(G, xi) the fugacities 

and to the edges of SAW(G, x\) the weights 

W(xi ■ ■■x k -i,x 1 ■■■x k ) = W(x k -i, x k ) . (6.2) 

We can then apply Algorithm T from Section EOl to calculate the effective fugacities 
w ^... x on the tree SAW(G, xi), working upwards from the leaves: 

Proposition 6.2 In the foregoing set-up, we have 

= -Rx k K Xk ,v(G)\{x 1 ,...,x k _ 1 }(R) (6.3) 

where the rational functions K Xj \(R) are defined by the recursion \5.2b\) ignoring the 
conditions K y ^\ x R y < 1. 



Proof. The proof runs inductively up from the leaves. 

If xi ■ ■ ■ Xk is a leaf (i.e. a maximal SAW), then w ^... x = w Xl ... Xk = —R Xk by 

Algorithm T and the definition (jfi.lj) : moreover, i^x fc ,v(G)\{xi,... ) a!j ! _i} = 1 because x k 
is an isolated vertex in G \ {xi, . . . , Xk-i} [since x\ ■ ■ ■ x k is a maximal SAW], so that 
the product ()5.26|) is empty. 

Let P = X\ ■ ■ ■ Xk be a non-leaf, and suppose that (|6.3|) holds for all one-step 
extensions P' = x~y ■ ■ ■ XkXk+i- Then Algorithm T gives 



J x 1 ---x k ~ w xi—x k II -. 



n 



x fc+1 eV(G)\{x 1 ,...,x fc } x 1 -x k x h+1 

1 - W(x k ,X k+ i)R Xk+1 K Xk+1 y(G)\{ Xl ,..., Xk } 



= -■R* fc ^ aifc ,V(G)\{xi,...,Xj t _i}, (6-4) 

where the second equality uses the inductive hypothesis, and the last equality follows 
immediately from ([5.26)1 . ■ 

Now let us apply Theorem 13. 2l to the tree SAW(G, X\) with the fugacities w defined 
by ()6.1|) and edge weights W defined by ()6.2j) . We conclude that the vector R defined 
by 

Rxt^xt = Rx k (6.5) 

lies in TZ{W) if and only if R x K XjA < 1 for all A C V(G) and all x E A. [By Q , 
this latter condition is sufficient to have 7£(jy); indeed, only certain sets A C V(G) 
actually arise. Conversely, given any A C V(G) and any x G A, we can choose 
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x\ — x and restrict to the subtree of SAW(G, x\) consisting of paths inside A. Using 
the monotonicity of 1Z(W) to restrict to this subtree, we conclude by Algorithm T 
that G (—1,0] and hence by ()6.3|) that R X K X ^ < 1.] On the other hand, by 
Theorem 15 . 61 this is a sufficient condition to have R G 1Z(W). 

In summary, the bounds produced by Theorem 15.61 for the lattice gas on G with 
edge weights W correspond to solving exactly (via Algorithm T) the lattice gas on 
the tree SAW(G, x-\) with edge weights W [for any Xi G V(G)]. And this produces a 
lower bound on the set TZ(W). 

6.2 Tree interpretation of the "optimal" bound: hard-core 

case 

Next let us consider the "optimal" recursion ()5.43|) in the special case of a hard- 
core pair interaction, i.e. the independent-set polynomial for the graph G. 

For each vertex x\ ■ ■ ■ Xk of the SAW-tree, let us choose (in any way we like) an 
ordering of the children x\ - ■ ■ x^x^+i- We then define the pruned SAW-tree corre- 
sponding to this ordering to be the subtree of SAW(G,Xi) whose vertex set consists 
of those paths x\ ■ ■ ■ Xk satisfying the rule that "you can't use an elder sibling of a 
vertex of G that you have previously used" . More precisely, for each path x\ ■ ■ ■ Xi, 
let us define the set S(xi ■ ■ ■ Xi) of "spurned vertices at step i" to be the set of all 
x G V(G) such that x\ ■ ■ ■ Xi-\x is a path that precedes x\ ■ ■ ■ x~i~yxi in the ordering of 
children of X\ ■ ■ ■ x~i~[. A path xi ■ ■ ■ Xk then belongs to the pruned SAW-tree if and 
only if for all 2 < i < j < k we have Xj ^ S(xi ■ ■ ■ Xi). 

Let us remark that the pruned SAW-tree is in general much smaller than the 
full SAW-tree. For example, if G = K n , then the SAW-tree has (n — 1)! leaves and 
(n — l)!^^ 1 V*' vertices, while the pruned SAW-tree (for any choice of ordering) 
has 2 n ~ 2 leaves and 2 n_1 vertices. 

Having defined the pruned SAW-tree, we then identify a path X\ - ■ ■ Xk in the 
pruned SAW-tree with the pair (xk, A) where 

k 

A = A(x x , ...,x k ) = V(G) \ {x 1 , . . . , x k -i} \ [JSixx ■■■x i ) , (6.6) 

i=2 

that is, V(G) minus the vertices already visited or already spurned. [Note that this 
differs from the definition A = V(G)\{xi, . . . , Xk-i} used in the preceding subsection.] 
Given a vector w = {w x } xe v(G) of fugacities on G, let us now assign to the vertices 
of the pruned SAW-tree the fugacities 

Wxj^ = w Xk (6.7) 

and, as before, assign to the edges the weights 

W(xl ■ --Xk-uxi ■■■x k ) = W{x k - U x k ) (6.8) 

(which, in the hard-core case currently under consideration, takes the value for each 
edge of the pruned SAW-tree) . We can then apply Algorithm T from Section 13.41 to 
the pruned SAW-tree to calculate the effective fugacities w ^... x : 
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Proposition 6.3 For all x\ 



■Xk belonging to the pruned SAW-tree, we have 



= w Xk )C XkMxii ... tXk) (w) . (6.9) 
Proof. The proof runs inductively up from the leaves. 



If x\ • • • Xk is a leaf (i.e. a maximal pruned SAW), then 



by Algorithm T and the definition (j6.7j) : moreover, fC Xkt A( xl ,..., Xk ) = 1 because Xk is 
an isolated vertex in the subgraph of G induced by A(xi, . . . , Xk) [since x\ ■ ■ ■ Xk is a 
maximal pruned SAW], so that all factors in the product (j3.14|) are 1. 

Let P = x\ ■ ■ ■ Xk be a non-leaf, and suppose that (|6.9|) holds for all one-step 
extensions P' = X\ ■ ■ ■ XkXk+i that are pruned SAWs. These one-step extensions are 
children of X\ ■ ■ ■ Xk, hence are ordered: let the ith such walk in order be X\ ■ ■ ■ XkUi 
(1 < i < I for some / > 1). Note that {yi, ...,?/;} = T(xk) fl A(x 1; . . . , Xk) and that 

h.(x 1 ,...,x k ,y i ) = A(xi, . . . ,x k ) \ {x k ,yi, . . . ,Vi-i} ■ (6.10) 

Then Algorithm T gives 



n 



i=1 x l"' x k% 



tl 1 + TU *i^i^(*i.".,Xfc,»i)( W ) 



1 



, I - + Uj W^'yi I A(a:i,...,a; fe )\{a' t .,yi,...,yi_i}( w ) 
= «V r A;,, .a,,,.. ,, v: (\V) (6.11) 

where the second equality uses the inductive hypothesis and the fact that W(xk, Vi) = 
0, the third equality uses (|6.1U|) . and the final equality uses (|3.14|) . I 

Now let us apply Theorem IH.2I to the prunedjSAW-tree with the fugacities w 
defined by ()6.7|) with w = — R, and edge weights W defined by ()6.8j) . We conclude 
that the vector R defined by 

RxT^xT = R Xk (6-12) 

lies in TZ(W) if and only if ^^^(-R) < 1 for all A C V{G) and all x e A; and by 
the discussion surrounding l|5.4H|) / (j5.44|) . this is a necessary and sufficient condition 
to have R G K(W). 

Therefore, we have shown that the "optimal" Dobrushin-Shearer bound a la 
(15. 4.3 J) / (|5.44jl for the independent-set polynomial (= hard-core lattice gas) on G cor- 
responds to computing exactly (via Algorithm T) the independent-set polynomial for 
the pruned SAW-tree [for any x\ € V(G) and any choice of orderings of children]. 
And this produces an exact computation of the set TZ(W). 
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6.3 Tree interpretation of the "optimal" bound: general case 



Finally, let us consider the "optimal" recursion ()5.43j) in the general case of a 
soft-core pair interaction W. Here we work again on the full SAW-tree SAW(G, xi), 
where G = Gw is the support graph of W\ "pruning" will be replaced by a "soft 
suppression of spurned vertices". As before, we begin by choosing (in any way we 
like) an ordering on the children of each vertex in the SAW-tree, and we use this 
ordering to define the set S(x{ ■ ■ ■ x~t) of "spurned vertices at step i". Then, given a 
vector w = {w x } x£ v(G) of fugacities on G, we assign to the vertices of the SAW-tree 
the modified fugacities 

w Xk | J W(xi-i,x k ) . (6.13) 

2 < i < k 
x k G S(xi ■ ■ -Xi) 

[Note that the product in (|6.13|) could equally well be written 2 < i < k, since 
Xk tf: S(xi ■ ■ ■ Xk).} As before, we assign to the edges the weights 

W(xi ■ ■ -x k -i,xi ■■■x k ) = W(x k -i, x k ) . (6.14) 

We can then apply Algorithm T to calculate the effective fugacities w ^... x on the 
tree SAW(G, xi). We obtain: 

Proposition 6.4 In the foregoing set-up, we have 

<^r = >>\;r ri/ ^n.v><:: In ,}(w r ^-" ) (6-15) 

where w[ xl '" Xfe l is defined by 

w [ y V " Xk] = v>v II Wix^y). (6.16) 

2 < i < k 
y £ S(xi ■ ■ - xi) 

Proof. The proof runs inductively up from the leaves. 
If x~i ■ ■ ■ Xk is a leaf (i.e. a maximal SAW), then 

w Xk J] W(x t . 1} x k ) = (6-17) 

2 < i < k 
x k g S(xi ■ ■ ■ Xi) 

where the first equality uses Algorithm T and the subsequent equalities are the defini- 
tions ()6.13|) and ()6.16|) . On the other hand, K, Xk y(G)\{ Xl ,...,x k _ 1 } = 1 because Xi ■ ■ ■ Xk 
is a maximal SAW, so that the product (|3.14j) is empty. 

Let P = xi ■ ■ ■ Xk be a non-leaf, and suppose that (j6.15|) holds for all one-step 
extensions P' = xl ■ ■ ■ XkXk+i- These one-step extensions are children of x\ ■ ■ ■ Xk, 



Xl---Xk 



''V;..,,-. = = 
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hence are ordered: let the ith such walk in order be x\ • • -XkUi (1 < z < Z for some 
I > 1). Note that . . . , yj} = r(x fc ) fl {xi, . . . , x^-i}. Then Algorithm T gives 



ir'T. = <<-,■,...,,,]_} 



•'' l '"' 7 ''' J A 1 1 + w2 s 



Note now that from (jfi.lfij) we have 

#[*i-* fcW ] = J ^(^^ y)w { y V " Xk] if y = %■ for some j < i 
^ Wy otherwise 

Therefore, the product on the right-hand side of (|bM8|) is identical to (|3.13|) if we set 
x = Xk and A = V(G) \ {xi, . . . ,Xk-i} and w is replaced by w^ Xl '" Xk \ for then (|6.19|) 
becomes precisely the vector w^ 1-1 ) defined in (|3.11|) . Therefore, (|6.18p equals the 
right-hand side of (|6.15jh as claimed. I 



7 Unfolding 

In Section 16.21 we showed that the "optimal" Dobrushin-Shearer bound a la 
fl5.43|l / (jf).44j) can be interpreted, in the special case of the independent-set polyno- 
mial (= hard-core lattice gas) for a graph G, as an exact computation of the set 1Z{G) 
based on computing exactly (via Algorithm T) the independent-set polynomial of the 
pruned SAW-tree of G. In this section we would like to show how this tree bound can 
be understood as arising from the repeated application of a single "unfolding" step. 

Let G = {V,E) be a graph, and let us select a pair of adjacent vertices x, y. As 
previously, G \ x (resp. G \y) denotes the graph obtained from G by deleting the 
vertex x (resp. y) and all edges incident with it. We then define G xy to be the graph 
obtained from the disjoint union of G' = G\y and G" = G \ x by adjoining an extra 
edge connecting the vertex x in G' to the vertex y in G" (Figure To each vertex 
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z in G \ x \ y = (G \x) \ y, there corresponds a vertex z' in G" and a vertex z" in 
G"'. Given a vector w with index set V(G), we define the corresponding "diagonal" 
vector w with index set V(G xy ) by 

w x = w x (7.1a) 

Wy = Wy (7-lb) 

w z , = w z (7.1c) 
w z ii = w z (7. Id) 

By considering the occupation of sites x and y, it is easily seen that 

Z G (w) = Z G \ x \ y (w) + w x Z G \ r * {x) (w) + WyZ G \ r * {y) (w) (7.2) 

and 

^g»»(w) = Z G \ x \ y (w) 2 + w x Z G \ x \ y (w)Z G \ r * {x) (w) + WyZ G \ x \y(w)Z G \ r * {y) (w) 

(7.3a) 

= ^G\x\ y (w) Z G (w) . (7.3b) 
We have the following fundamental result: 

Theorem 7.1 R G 7£(G) and on/y z/ R G lZ(G xy ). In other words, 71(G) is 
(isomorphic to) the "diagonal cross section" oflZ(G xy ). 



Proof. Suppose first that R e 71(G). Then we also have R G 7Z(G \x\y) by the 
monotonicity statement at the end of Proposition ^. 151 Therefore, if — R < w < we 
have Z G (w) > and Z G \ x \ y (w) > 0, and hence deduce by ()7.3b|) that Zq xv (w) > 0. 
In particular, Zq xv (— AR) > for < A < 1. Applying Theorem 12.1 Uf a) (b) 
to the line segment connecting to — R, we conclude that Zq xv (w') > whenever 
— R < w' < (whether or not w' lies "on the diagonal"). Hence R G lZ(G xy ). 

Conversely, suppose that R G 7Z(G xy ). If — R < w < we have — R < w < 
and hence Zg xy (w) > 0; by ()7.3b|) this implies Z G (w) ^ 0. I 

Assume now that G is connected. Let G xy be the component of G xy containing x 
and y. Any other component of G xy must be contained either in G' \ x or in G" \ y; 
either way, it corresponds to some component of G \ x \ y. Now, for any component 
H of G \ x \ y, there are three possibilities: either H is adjacent to x in G, or H is 
adjacent to y in G, or both. In the first case, the copy H' of H in G' is contained 
in G xy , while the copy H" of H in G" is disjoint from G xy ; in the second case, the 
reverse holds; in the third case, both copies of H are contained in G xy . Therefore, any 
component of G xy other than G xy has a mirror image contained in G xy , and hence is 
isomorphic to a subgraph of G xy (with the same weights when we are "on-diagonal" ) . 

Given a vector w with index set V(G), let us define the vector w with index set 
G xy by restricting w from V(G xy ) to V(G xy ). 
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Lemma 7.2 Let H be a subgraph of G, and let R = {R x }weV(G) £ 71(G). Then 



Proof. Define first K1 V (h) [with index set V(G)} by 

(R1 V(H) ) X - | Q i{x £ V{H) (7-4) 

Since H1v(h) < R, we have Rly(/ri £ ^-(C) because 7^(G) is a down-set. Equiva- 
lently, R# £ 72.(G[V(H)]), where G[V(if)] is the induced subgraph of G with vertex 
set V(H). But 71(H) D 7l(G[V(H )]) by the monotonicity statement at the end of 
Proposition 12. 151 I 

Corollary 7.3 Assume that G is connected. Then R £ 71(G) if and only if R £ 
7Z(G xy ) . 

Proof. If R £ then R £ TZ(G xy ) by Theorem O Then R £ TZ(G xy ) by 

Lemma \7.2\ because G xy is a subgraph of G xy . 

Conversely, suppose that R £ 7Z(G xy ). Let Hi, . . . , Hk be the components of G 3 ^ 
other than Then 

k 

Zq xv (w) = Z Qxv (w) \\Z Hi (<n \ Hi) (7.5) 

i=l 

If — R < w < 0, then — R < w < and hence Zg xy (w) > 0. But since each 
component of Hi is isomorphic to a subgraph of G xy (with the same weights in the 
mirror copy), Lemma (7.21 implies that Z^(w |~ Hi) > as well. Hence Zq xv (w) > 
by O, so that Z G (w) ^ by (j7~3bT) . This shows that R £ 7^(G). | 

Remark. It would be interesting to know whether there is a version of unfolding in 
the case of soft-core interaction. 

8 Infinite graphs 

In this section we discuss briefly the repulsive lattice gas on a countably infinite 
graph. We begin by deriving some general properties valid on an arbitrary countably 
infinite graph (Section IHHJ). Next we discuss two cases of special interest: trees (Sec- 
tion IHI2I) and regular lattices (Section I8.3J) . Finally we show how quantitative bounds 
on 7Z(W) can be obtained, using as an example the square lattice 1? ( Section l8.4|) . 
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8.1 General properties 

Let X be a countably infinite set, and let W: X x X — > [0, 1] be symmetric. For 
every nonempty finite subset A C X we can consider the partition function Z A (w) 
defined for w £ C A in the obvious way, i.e. by considering the lattice gas on A with 
interaction W \ A. Then, to each such A there corresponds a set 1Z(W \ A) C [0, oo) A . 
From the Remark after Proposition 12.151 we observe that if A' C A and R £ [0, oo) A , 
then 

R £ K{W \ A') ^ (R, 0) £ K(W \ A) (8.1) 

(here has index set A \ A'). 

Now let us define, for each finite Acl, the set 

TZ A (W) = TZ(W \ A) x [0,oo) x \ A (8.2) 

in the infinite-dimensional space [0, oo) x . Otherwise put, a vector R £ [0, oo) x 
belongs to K A {W) if and only if R \ A lies in TZ(W \ A). Note that 

&R A (W) = dK(W \ A) x [0,oo) xvv . (8.3) 

Clearly TZ A (W) is open (in the product topology) and is a down-set. Moreover, it 
follows immediately from ()8.1|) and the fact that 1Z{W \ A) is a down-set that if 
A' C A, then 1Z A '{W) ^ 1Z\(W). In other words, the {7Z A (W)} form a decreasing 
family of sets in [0, oo) x [when A runs over the collection of finite subsets of X ordered 
by inclusion]. We define the limiting set 

K(W) = f]n A (W) . (8.4) 

A 

Note that, because {7Z A (W)} is a decreasing family, we also have 

oo 

n(w) = f)n An (w) (8.5) 

n=l 

for any increasing sequence A^ C A 2 C . . . whose union is all of X. 

Remark. If W(x,x) < 1 for all x £ X, then it follows from Proposition 12.161 that 
there exist constants C x < oo such that 1Z(W) C [X,. gX [0, C x ], so that 1Z{W) is 
relatively compact in the product topology. 

In view of the fact that each set 1Z A {W) is open, it is perhaps surprising that the 
limiting set 1Z{W) is "almost" closed: 

Theorem 8.1 Let X be countably infinite, and let R £ 1Z(W) (where the closure is 
taken in the product topology). If every component of G\y [supp R] is infinite, then 
ReK{W). 

Corollary 8.2 Let X be countably infinite, let r £ [0, oo) x , and suppose that every 
component o/Gvy[suppr] is infinite. Then {A > 0: Ar £ 7Z(W)} is a closed interval 
of [0, oo) (which may reduce to {0} or to all of [0, oo) ). 
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Remark. If W(x, y) < 1 for some x, y G supp r, then it follows from Proposition ^. 161 
that the interval {A > 0: Ar G 1Z{W)} is bounded, i.e. not all of [0, oo). 

Proof of Theorem 18.11 By hypothesis we have 

r g f]n A (w) c f]njw). (8.6) 

A A 

If R G 7Z\(W) for all A, we are done. So assume that R G <97£a(WO for some A, 
i.e. R [ A G <9ft(W C A) and in particular Z A (-R) = by Proposition HZEff b) . Let 
Ai, . . . , Afc be the vertex sets of the components of GVfsupp R n A]. Since Z\(w) = 
Yli=i ^Aj(w) whenever supp w C supp R, there must be at least one index i for which 
Z Ai (-R) = 0, so that R f Aj G &R(W \ Ai) [recalling that R \ Aj G K{W \ Aj) by 
()8.6)) ]. Now choose a vertex x G (suppR) \ Aj that is adjacent to Aj in Gyy [this is 
possible since each component of GV[suppR] is infinite]; and define Aj = A.; U {x}. 

By USED we hav e R [ Aj G TZ(W \ Aj); and hence by Corollary (PH we have (R \ 
Aj, 0) G 7l[W \ Aj) [where corresponds to the xth entry], i.e. R \ Aj G TZ(W \ Aj). 
But this contradicts R \ Aj G [ Aj) since TZ(W \ Aj) is open. I 



8.2 Infinite trees 

In Example El 6 we considered the complete r-ary rooted tree and showed, following 
Shearer pj^, that there are negative real roots of the univariate polynomial Zq(w) 
that tend to Wqq = —r r /{r + l) r+1 as the number of levels tends to infinity. Let us 
now use the "homogenization" ideas of Section to extend this result to trees that 
are not complete. 

Proposition 8.3 Let G be an infinite tree with (arbitrarily chosen) root vertex Xq, 
and let X { = {x: dist(a;,a; ) = i}. Let b = limsup {X^ 1 ^. If R > b b /{b + l) b+1 , then 

/ \VPQI 

any vector R > satisfying I FJ R x J > R for all i does not lie in 1Z(G). In 



particular, Rl 71(G). 



Proof. Let R > be a vector satisfying R, L = dl^gx RxY^ X ^ > R for all i, and 
suppose that R G 1Z{G). Then by definition we have R G IZ(Gd) for every D, where 
Go = G[[Ji<D Xi] is the subtree consisting of the first D + 1 levels of G. 

Let us fix D, and apply Algorithm T of Section to the tree Gd with w = — R. 
For < i < D, define Pi as in (|3.52|) . Let 6j = |Xj +1 |/|Aj| and define % by the 
recursion (|3.57|) . Then it follows from Proposition 13.31 [or the discussion after (|3.57|) ] 
that Pi > Pi for < i < D, since by hypothesis R G H{Gn). By monotonicity of 
(I3.57|) in the i?j, it follows that the sequence (q\ D ^)°^i defined by the recursion 



1 " Q: 



i+l) 
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with initial condition q\ D ^ = for i > D + 1 satisfies < qf ^ < Pi for < % < D. In 
particular, we have < < 1 for every i. 

Now let D — > oo. Monotonicity of ()8.7|) in q\®\ implies that, for each i, q\ D ^ is 
increasing in D. Since g- D ^ is bounded above by 1, taking g« = linio^oo (if " 1 gives a 
sequence (qi)fZo satisfying 

Qi = 7 ^—rr (8-8) 

for every i. Furthermore, < qi < 1 for every and so (j8.8|) implies that < q,[ < 1 
for i > 1. 

The proposition is then a consequence of the following lemma. ■ 
Lemma 8.4 Let (bi)^l l be a sequence of positive real numbers, let b = limsup(oi02 • • • h 

i— >oo 

and let R > 0. Suppose that there exists a sequence (gj)i^i satisfying < Oj < 1 and 

* £ ra* (8 - 9) 

/or alli>\. Then R < b l /{b + 

As preparation for proving Lemma f8 .41 let us prove an analogous result for finite 
sequences. 

/ k \ V* 

Lemma 8.5 Let k > 1, let 6 , 6i, . . . , 6fc_i > ; and define b = I Y[bi ) ■ Let 

\i=0 J 

R > 0. Suppose that there exist q , q x , . . . , q k e [0, 1) satisfying 

R_ 

and 

go < 9fc • (8.11) 

T/ien # < 6 5 /(6 + 

Proof. When k = 1, the lemma is proved by straightforward calculus. So let us 
treat the case k > 1. Define weights 

for i = 0, 1, . . . , k and note that 70 = 7fc = 1- Let T = Ti- Now define 



> - 7-^77 fori = 0,l,...,k-l (8.10) 



go 



l[q?) (8.13a) 



v i=0 



91 = ( 8 - i3b ) 



v i=l 
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The argument at ()3.65|) - ()3.73|) shows that 

qo > j— ~ ■ (8.14) 
(1 - Qi) b 

If Qo S Qki then % < q\\ it then follows from the k = 1 case of the lemma that 
R < b b /{b + l) b+1 . I 

Proof of Lemma 18.41 If R > b b /(b + l) b+1 , let us choose e > and c < b such 
that 

R c e b b 

TT~e > (c+l) 5 +! > (6+l)6+i • t 8 ' 15 ) 

Since limsup i ^ 0O (6i6 2 • • • ^i-i) 1 ^ > c, we can find integers 1 < jx < 3% < ■ ■ ■ with 

jt+i 

[I bt > e ?t+1 - it (8.16) 

i=jt+X 

for t = 1,2, .. . . Since the sequence (qj t )^i lies in the interval [R, 1) and R > 0, 
there must exist £ such that g Jt < (1 + e)qj t+1 . To simplify the notation, let us set 
k = jt+i — jt and % = qj t+i for i = 0, . . . , k. By hypothesis we have 

ft > 7; 5 ^ fori = 0,l,...,fc-l (8.17) 

(1 - g 4+ i) b ' 

and §o < (1 + e)gfc. It follows that if we set S 1 = + e), ^ = qi for i = 1, . . . , k, 
and g = qo/(l + e), we have 

%>- i—- fori = 0,l,...,A:-l (8.18) 

and g < Lemma 1831 together with ()8.16|) then imply that S < c c /(c+l) c+1 , 
which is a contradiction. ■ 



8.3 Regular lattices 

The most important situation in statistical mechanics is that of a model defined 
on a regular lattice (which we take for simplicity to be Z d ) with a translation-invariant 
interaction [lOOj . In this case one expects to be able to prove that the free energy per 
site (or "pressure") Fa = |A| _1 logZA converges to an infinite- volume limit 

Foo = lim F A , (8.19) 

A /"oo 

where A y oo denotes convergence in the F0lner-van Hove sense, i.e. |A| —>■ oo in 
such a way that the surface-to-volume ratio tends to zero (see |113| Section 2.4.1 
and Appendix A. 3.1] for a variety of equivalent conditions). Indeed, there are several 
standard arguments for proving such convergence: 
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1) Almost additivity. A large volume A is subdivided into smaller (but still large) 
cubes Aj separated by wide "corridors" , and | log — £\ log Z Ai | is shown to be 
suitably small jHZl Section 1.2]. This method applies to all models with bounded 
interaction energies, and yields F0lner-van Hove convergence [67, Theorems 
1.2. 3-1.2.5]. 

2) Superadditivity. In some cases it can be shown that log Z\ u \* > log Z\ + log Z^ 
whenever A and A' are disjoint. This occurs, in particular, for systems with 
certain symmetries, for ferromagnets, and for systems with negative interac- 
tion energies. Moreover, arbitrary models with bounded interaction energies 
can be reduced to the case of negative interaction energies simply by adding 
suitable constants to the interaction terms |1UU[ Theorem II. 2. 4 and Examples 
1-4 following it; see also Section II. 4]. However, this method does not yield 
F0lner-van Hove convergence, but only a slightly weaker convergence in which 
|A|/diam(A) d must be bounded below away from zero [1 131 Appendix A. 3. 3]. 
Moreover, the limiting free energy F^ could be +oo; a separate argument is 
needed to exclude this possibility. 

3) Subadditivity . In some cases it can be shown that log^AuA' < log^A + log^A' 
whenever A and A' are disjoint. This occurs, in particular, for systems with pos- 
itive interaction energies. Moreover, arbitrary models with bounded interaction 
energies can be reduced to the case of positive interaction energies by adding 
suitable constants to the interaction terms. Like superadditivity, this method 
does not yield F0lner-van Hove convergence, but only the slightly weaker con- 
vergence in which |A|/diam(A) d is bounded below away from zero. Moreover, 
the limiting free energy F^ could be — oo; a separate argument is needed to 
exclude this possibility. 

Remark. Although we shall concentrate here on the case of a regular lattice 
(namely, Z d ), arguments of the foregoing types can usually be generalized to handle 
arbitrary quasi-transitive amenable infinite graphs |6TH I91j; we expect that the same 
should be true for our results. 

Let us therefore consider the lattice gas on X = Z d with a translation-invariant 
interaction W(x, y) = W(x — y) satisfying < W(x, y) < 1. Let us place the 
same fugacity w > at each site. [More generally, we could consider a periodic 
or quasiperiodic fugacity w = {w x } xeZ d; but let us stick to a constant fugacity for 
simplicity] Under these assumptions we can prove that the infinite-volume limit of 
the free energy exists. Indeed, the subadditivity argument applies immediately to our 
model, thanks to the hypothesis of repulsive interactions; moreover, since Z\ > 1 for 
all A (thanks to the contribution of the empty configuration), it follows that F^ > 
and in particular that F^ 7^—00. But this method does not yield F0lner-van Hove 
convergence. 

To prove F0lner-van Hove convergence, we use the almost-additivity argument. 
The standard theorems [H7J Theorems 1. 2. 3-1. 2. 5] do not apply to our model, because 
the interaction energies can be unbounded, indeed for either of two reasons: 
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(a) the infinite interaction energy when W(x, y) = for some x ^ y; and/or 

(b) the unboundedness of the interaction energies in the absence of hard-core self- 
repulsion [i.e. for W(x, x) 7^ 0], arising from the fact that the occupation number 
n x at a site can be arbitrary large. 

We can nevertheless make slight adaptations in the standard almost-additivity argu- 
ment so as to make the proof go through. We assume that the total interaction of 
any site with the rest of the world is finite, i.e. 



J2ll-W(x,y)} < 00. 



(8.20) 



In particular, there is a number R < 00 such that W(x, y) > whenever \x — y\ > R. 

Let us begin with the case of hard-core self-repulsion, so that we need only deal 
with problem (a). This is handled by the following lemma: 

Lemma 8.6 Consider any repulsive lattice gas with hard-core self-repulsion and fu- 

n 

gacities w x > 0. Suppose that A = {J Aj U A (disjoint union). Then 

i=l 



.i=i 



\ 



1 < i < j < n 

x e Aj 



< Z A < [Y[Z Ai ] ( Y[(l + w x ) ) . (8.21) 

\x€A 



.1=1 



Proof. We make use of the definition 



Za = E Il w * II w ^y)- 

X'CA x£X> {x,y}CX> 



122) 



The lower bound is on Z A is obtained by considering only those configurations in 
which A is empty (i.e. X' n A = 0) and writing 



n w(x lV ) > m n w( X , V ) 

{x,y}QX> \ i=l {x,y}CX'nX t 



\ 



n n w ^ 

l<i<j<n x £ Ai 



, (8.23) 



which is valid since < W(x, y) < 1. The upper bound on Z A is obtained by writing 

n 

n w(x, y ) < n n w( X , V ). (8.24) 

{x,y}CX> i=l {x,y}CX'nX t 



We then have: 
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Theorem 8.7 (infinite-volume limit, hard-core case) Consider a translation-invariant 
repulsive lattice gas onZ d , with hard-core self-repulsion, satisfying J2 y€Zd [l—W(x,y)} < 
oo, with the same fugacity w > at each site. Then F^w) = liniA/oo log Z A (w) 
exists in F0lner-van Hove sense and satisfies < F^ivj) < log(l + w). Moreover, 
-Poo( w ) is an increasing and convex function of log w; in particular, it is continuous 
on < w < oo. 



Proof. Fix integers a, c > and consider the paving of Z d by disjoint cubes of side 
a + c with corners located at (a + c)Z d , i.e. cubes C n = [0, a + c) d + (a + c)n for n G Z d . 
Consider also the subcubes C' n C C n of side a with the same lowermost corner, i.e. 
C' n = [0, a) d + (a + c)n. For any finite subset A C Z d , let C\ be the collection of all 
cubes C n that are contained in A, and let C' A = {C' n : C n G Ca}. By Lemma | 



logZ A - ^2 ^°S z c 



C'&C' A 



< a{c)\A\ + (3 



A\ |J C 



(8.25) 



where a(c) = — logW(a;) and (3 = log(l + w). By translation invariance, 

x&Z d 

\x\ > c 

log Zc> = logZ C ' for all C G C' A . Moreover, \C' A \ = \C\\ by construction. Dividing 
by |A|, we°get 



^log^A 



and hence 



< a(c) + /3 1 



a' 1 


\c 


A 




A 





(8.26) 



1 - 



a d \Cs 



|A| 



/3 + -r|logZ, 



< a(c) +2/3 1 







A 




A 





(8.27a) 
(8.27b) 



since 1 < Zcj, < (1 + w) a . Taking ()8.27bj) for two finite subsets A, A C Z d and 
subtracting, we get 



\Fa~F a \ < 2a(c) + 2/3(1- 



a d \C A 
|A| 



2/3(1 ^ 



|A| 



(8.28) 



Now let (A n ) oo in F0lner-van Hove sense; this implies that the fraction of 
volume of A n contained in the cubes C\ n tends to 1, i.e. lim n ^ 00 (a + c) d |CA,J/|A n | = 1. 
Therefore, 

limsup \F Am — F An | < 2a(c) + 4/? 1 - [ )<1 



m,n—>co 



a + c 



(8.29) 
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Now, given any e > 0, we can [by ()8.20|) ] choose c large enough so that a(c) < e; and 
we can then choose a large enough so that 1 — [a/ (a + c)] d < e. Hence 

limsup|F Am -F An \ < 2e + 4/3e. (8.30) 

m,n— >oo 

Since e is arbitrary, we have lim sup m n ^ co \F\ m —F\ n | = 0, so that (|A n | _1 log Z\ n )^ =1 
is a Cauchy sequence and hence converges. (The limit is the same for all F0lner-van 
Hove sequences, since interleaving two F0lner-van Hove sequences yields another one.) 

The bounds < F OQ (w) < log(l + w) follow immediately from 1 < Z\ < (l+w)' A L 
Finally, each Fx{w) is manifestly increasing, and by Lemma r2.32l it is a convex function 
of logu>; and these properties are preserved under pointwise limits. ■ 

Now let us turn to problem (b), arising from the unboundedness of the occupation 
number n x when W(x,x) ^ 0. We begin with a well-known fact about positive 
correlation of increasing functions on a totally ordered space: 

Lemma 8.8 Let \i ^ be a nonnegative measure on a totally ordered space Q, and 
let f and g be increasing functions on Q. Then (provided the functions concerned are 
integrable) : 

(a) (Jfgd/j) (Jlrf/i) > (ffdn)(fgdn) . 
f fe° 9 d\x 

(b) —p — is an increasing function of a G M. 

Proof, (a) We have 

(ffgdfi) (Jldfx) - (Jfdfi) (fgdfx) = - J [f{x)-f{y)] [g(x) - g(y)]dfx{x) d/i(y) , 

(8-31) 

and the integrand on the right-hand side is nonnegative both when x < y and when 
x > y. 

(b) We have 

d fjfe^d^ _ (J fge^dy){je^dy) - (J fe<* » d^jf ge™ d/i) 



da V / e a 3 dp J (J e a 9 dfx) 2 

Now apply part (a) with ji replaced by e a9 fi. I 

Now let A be a finite set; for each x G A, let Q x be a totally ordered space; and 
let Q = Ylxex ^x- F° r eacri x G A, let ^ ^ be a nonnegative measure on fl x and 
let F x be a nonnegative decreasing function on Q x . Finally, let H be an increasing 
function on Q and let a > 0. Define 

Z { A {F * ] > a > H) = f (]jF x ( Vx ))ex P [-aH( V )]l[d^ x ) . (8.33) 
J \xeA J zeA 
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Lemma 8.9 Under the above hypotheses, 

J Fx dp-. 
f 1 da x 



Proof. Choose one site z e A, and let us study the integral over ip z with </?^ 2 held 
fixed. By Lemma f8.8f b). the quantity 

J F z {ip z ) exp[-aH{if z ,!f^ z )}dn z {ip z ) ^ 
J exp[-<x£f (ip z , ip^ z )] d/i z (<p z ) 

is an increasing function of a 6 1, so its value at a > is bounded below by its value 
at a = 0, which is 

/ F z (ip z ) d[i z (ip z ) 



fdftzivz) 

Therefore, 



(8.36) 



fa* 



J F z ((p z ) exp[-aH((p z ,(p^ z )] d/j, z ((p z ) > (Jj^^J^j x J exp[-aH((p z , <p^ z )] d\i t 

(8.37) 

Now multiply both sides by rizgA\z F x (ip x ) and integrate with respect to rizeA\2 dfi x (ip x )] 
we obtain 

z ({F x}> a,H) > ( J^^A x ^{l-^WO (8 38) 

V J 1 d/x 2 / 

Applying the same argument successively to each site in A, we obtain (J8.34|) . ■ 

Let us now specialize these results to the repulsive lattice gas, by taking Q x = N, 
dfi x (n x ) = w x x /n x \, 

H(n) = Y,[-^gW(x,x)] nx( - nx 2 ~ 1) + Yl [-logW(x,y))n x n y , (8.39) 

x<=A ~ {a;,j/}CA 

and F x (n x ) = l(n x < K x ) for arbitrarily chosen positive constants K x . We then have 

Z (K) 

-A- > JJ[l- 7 (lf, lWjB )], (8.40) 

A xeA 

where Z^"' denotes the sum (jl.lbj) restricted to the configurations satisfying n x < K x 
for all x, and 



oo r 



j(K,w) = e- w J2 — . (8.41) 
We can now prove a variant of Lemma 18.61 



n=K 
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Lemma 8.10 Consider any repulsive lattice gas with fugacities w x > 0, and let 

n 

(K x ) x< z\ be arbitrary positive constants. Suppose that A = |J A, U A (disjoint union). 

8=1 

Then 



( 



n 

v« = l 



/ 



]] [l-l(K x ,w x )} 



l<i<n 



/ 



\ 



n w ^vY 



1 < i < j < n 
x g A 4 



< z A < iim n 



(8.42) 



vrrgA 



Proof. The upper bound is proved exactly as in Lemma f8. 61 For the lower bound, 
we begin from the trivial fact that Z\ > Now consider only those configurations 

satisfying n < K and, in addition, n x = for all x G A ; for such configurations we 
have 



{x,j/}CA\A \i=l {x,y}CX> 



n n w(x, y ) K * K y 

l<i<j<n Ig A, 



It follows that 



zf) > 



(K |"Aj) 



v i=l 



n n w ^yy 



l<i<j<n x 6 A, 
2/ g A, 



/ 

(8.43) 



(8.44) 



Now use (|8.40p for each set Aj. 



Theorem 8.11 (infinite-volume limit, general case) Consider a translation-invariant 
repulsive lattice gas on 7L d , satisfying ^2 y&Z d[^- — W(x,y)] < oo, with the same fug acity 
w > at each site. Then Foo^w) = liniAz-oo | yV | 1 \ogZ\(w) exists in F0lner-van Hove 
sense and satisfies < F^w) < w. Moreover, F^iw) is an increasing and convex 
function of log w; in particular, it is continuous on < w < oo. 

Proof. We begin by defining families of cubes Ca and C' A as in the proof of Theo- 
rem IHIZI By Lemma f8. 101 with K x = K for all x, we have 



logZy\ 



^2 lo £ Z c 



< a(c)K 2 \A\ + w 



A\ |J C 



cec' 



+ \A\\log[l ->y(K,w) 



$.45) 
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where a(c) 



^ logW(x). By the same arguments as before, we have 

\x\ > c 



F a I < a(c)K 2 + 2w 1 - 



|A| 



log[l - 7(fT, w) 



(8.46) 



Taking ()8.46j) for two finite subsets A, A C Z d and subtracting, we get 



Fa~F x \ < 2a(c)K 2 + 2w[l- 



a d \C A 
|A| 



+ 2w\l- 



a d \C~A 



|A| 



+ 2|log[l-7(if,«>) 



(8.47) 



Now let (A n ) y oo in F0lner-van Hove sense, so that lim n _ >00 (a + c) d |CA n |/|A n | 
1. It follows that 



limsup|F Am -F A J < 2a{c)K 2 + Aw 

m,n— >oo 



a + C 



2|log[l- 7 (*»]| . (8.48) 



Now, given any e > 0, we first choose K large enough so that | log[l — j(K, w)]\ < e; 
then we choose c large enough so that a(c)K 2 < e; and finally we choose a large 
enough so that 1 — [a/ (a + c)] d < e. Therefore, 



lim sup \F Am - F An | < 4e + Aew . 
The remainder of the argument is exactly as in Theorem 18.71 



.49) 



Historical remark. Results like Theorem 18 . 1 1 1 have been proven in the vastly more 
general context of "unbounded spin systems" by Lebowitz and Presutti [73], based 
on superstability estimates due to Ruelle [03] • See also [7H E] for related work. We 
think that it is nevertheless useful to give an elementary and self-contained proof for 
our special case. 

Let us now return to the case of hard-core self-repulsion, and consider the conver- 
gence of the finite- volume free energies F\(w) for complex fugacities w. By hypothesis 
dSIZD) we have 

A w = £)[1-W(0,j/)] < oo. (8.50) 

Therefore, by Corollary 15.21 all of the partition functions Z\(w) are nonvanishing 
in the disc \w\ < R(Aw) where R(Aw) is defined by (|5.12|) . Moreover, they sat- 
isfy the trivial bounds |Za(u>)| < e' A " w '. It follows that the free energies F^(w) = 
|A| _1 log Z\(w) are analytic in the disc \w\ < R(A W ) and satisfy there ReF\(w) < 
\w\ < R(Aw). Finally, Theorem 18.111 shows that the F\(w) converge to a limit when 
w lies in the real interval [0, R(Aw)). These three facts are sufficient to imply the 
convergence of F\(w) to an analytic limit F^w) everywhere in the disc \w\ < R(A W ), 
using the following standard result on normal families of analytic functions: 
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Proposition 8.12 (exp log Vitali) Let D be a domain in C, let S be a subset of D 
having at least one accumulation point in D, let M < oo, and let (f n ) ( £ = i be analytic 
functions in D satisfying: 

(a) Hef n (z) < M for all n and all z G D; and 

(b) lim fn(z) exists (and is finite) for all z G S . 

Then there exists an analytic function on D such that f n (z) — > fao(z) uniformly 
for z in compact subsets of D. 

For a proof of Proposition 18. 121 see e -g- |HH1 P- 343]. 

Corollary 8.13 Consider a translation-invariant repulsive lattice gas on Z d , with 
hard-core self -repulsion, with 

Y^[i--W(0,y)} = A w < oo. (8.51) 

Then: 

(a) Each function F\(w) = (A)" 1 log Z\(w) is analytic in the disc \w\ < R(A W ) 
and satisfies there ReF^m) < \w\. 

(b) There exists an analytic function Fooiw) on the disc \w\ < R(A\y) such that 
liniA/oo F\{w) = Foo^w) uniformly for w in compact subsets of \w\ < R(Aw)- 

Remark. Proposition 18.121 is a very special case of a much more general "Vitali- 
Porter-type" theorem for normal families of analytic functions [103, Lemma 3.5], in 
which condition (a) can be weakened to "there exists a nonempty disc A C C such 
that f n (z) A for all z G D" jUHl Example 2.3.9] or even to "there exist Wi,w 2 G C 
with W\ 7^ W2 such that f n (z) ^ {^1,^2} for all z G D" [SHI Section 2.7]. Moreover, 
it is sufficient for this hypothesis to hold locally in D. Detailed accounts of these 
results can be found in l9*o]. 

8.4 Quantitative bounds for the lattice Z d 

By Corollary EH the set {A > 0: Al G U(W)} is a closed interval [0, A c ] provided 
that each component of Gw is infinite. In this section we sketch briefly some methods 
for finding reasonably sharp upper and lower bounds on A c . We shall illustrate these 
examples with reference to the hard-core lattice gas on 7L d with nearest-neighbor 
edges, paying particular attention to the case of the square lattice Z 2 . 
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8.4.1 Upper bounds on A c 

Let G be any finite subgraph of Z d , and let w* = —A* be the negative real root 
of Zq of smallest magnitude. Then, by the definition of 1Z(W) for an infinite graph 
[cf. f!8.4|l ] together with the monotonicity of 7Z{G) in G [Proposition 12.15] . we can 
conclude immediately that A c (Z d ) < A*. Moreover, these bounds converge to the 
exact value A c if we take any increasing sequence G\ C Gi C . . . whose union is all 
of Z d . 

More generally, as a consequence of the discussion in Section IB~2*l we can take any 
finite subtree T of the pruned SAW-tree for Z d , and let w+ = —A* be the negative 
real root of Zt of smallest magnitude; we again have A c (Z d ) < A*. 

Example 8.1. Clearly Z d contains finite paths of all lengths. By Example El 2 
(or alternatively Example El6 with r = 1), these give bounds A* tending to 1/4 as the 
path length tends to infinity. Thus A c (Z d ) < 1/4. 

Example 8.2. Now consider a finite subgraph G C Z d consisting of a long path 
P along one axis (the spine) and, radiating from each vertex of P, 2d — 2 disjoint 
long paths perpendicular to P (the antennae). Let us assign fugacity w < to every 
vertex of G. Using Algorithm T on the antennae, we must first iterate the recursion 
w c s i — ► w /(I +w c fi). If the antennae were infinitely long, we would approach the fixed 



point w' = —\ J r\J wJ r\ [provided that —1/4 < w < 0; otherwise we are outside 
71(G)]. By taking the antennae sufficiently long, we can get as close to this value as 
we wish. We are now left with a "caterpillar" consisting of the spine (with fugacity w 
on each vertex) along with 2d — 2 pendant vertices (each with fugacity w') attached 
to each vertex of the spine. Applying Algorithm T to this graph, we get the recursion 
w cS h-» w/[(l + w e g)(l + w') 2d - 2 ] = io*/(l + ^eff), where 

_ w w 

w * = (TT^p = 7i / P M " 2 ' (8 ' 52) 

We require — | < w+ < in order to stay within 1Z(G). For d = 2 this yields the 
bound A C (Z 2 ) < 4/25. For d -> oo it yields A c (Z d ) < (logd)/(2d) + 0((loglogd)/d). 

Example 8.3. To get the correct bound A c (Z d ) = 0(l/d), one can argue using 
a subtree of the pruned SAW-tree. Note first that any walk using only steps in the 
positive coordinate directions is guaranteed to be self-avoiding. Moreover, if we define 
the pruning such that at each vertex all steps in the positive coordinate directions are 
preferred to all steps in the negative coordinate directions, then every walk using only 
positive coordinate steps appears as a vertex in the pruned SAW-tree of Z d . Thus, 
the complete <i-ary rooted tree (consisting of these walks) is a subtree of the pruned 
SAW-tree of Z d , and so by Example El 6 we have 

d d 1 

^ ' * (ITW* ~ 3 ' (8 ' 53) 
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Example 8.4. Asymptotically correct upper bounds on A c (Z d ) can be obtained by 
using large cylinders, for which A c can be computed by the transfer-matrix method. 15 
Let us illustrate the method for Z 2 . Consider the strip Sl = {(x,y) G Z 2 : < 
x < L}. Since S t C S 2 C . . . C Z 2 , we have A c (Si) > X C (S 2 ) > ... > A C (Z 2 ). 
[In particular, Hni5-»oo X c (Sl) exists.] On the other hand, since each finite subgraph 
of Z 2 is contained (modulo translation) in Sl for all sufficiently large L, we have 
lim^oo X C (S L ) < A C (Z 2 ). It follows that X C (S L ) I A C (Z 2 ) as L -> oo. 

An analogous argument can also be made using strips with periodic boundary 
conditions, which are more convenient for computation [110j. To see this, let Sl be 
the strip Sl with an extra edge added from (L — 1, y) to (0, y) for each y. As before, 
each finite subgraph of Z 2 is contained (modulo translation) in Sl for all sufficiently 
large L, so that limsup^^ X c (Sl) < A C (Z 2 ). On the other hand, the pruned SAW- 
tree of Sl (for any given choice of ordering) is a subtree of the pruned SAW-tree 
of Z 2 (provided we make an appropriately compatible choice of ordering): just map 
each path in Sl to its "universal cover" in Z 2 , i.e. to the path in Z 2 obtained by 
making the same sequence of north, south, east or west steps. Thus X c (Sl) > A C (Z 2 ) 
for every L (though we do not necessarily have monotonicity in L). It follows that 
X c {Sl) — * A C (Z 2 ) as L — ► oo. 

Example 8.5. Let us also remark that A c (Z d ) is strictly decreasing in d; indeed, 
we can derive an inequality bounding A c (Z a!+1 ) above in terms of A c (Z d ). To see this, 
note that Z d+1 contains a copy of the graph Gd obtained from Z d by attaching two 
semi-infinite paths ( "antennae" ) to every vertex. If we place fugacity w < on every 
vertex of Gd and integrate out the antennae as in Example 8.2, we are left with a 
copy of Z d with effective fugacities 

W W /o r A\ 

WeS = (l + w>)* = — \ — r=i ■ ( } 



X c (Z d ) > , (8.55) 



If w e (-A c (Z d+1 ), 0], then w e{i must be in (-A c (Z d ), 0]. It follows that 

A c (Z d+1 ) 

I _ Ac(Z <m) + ^1 _ Ac (Z<w 

or equivalently 

which is the desired bound. 

Using (J8.56j) together with the initial condition A C (Z X ) = 1/4, it is easy to show 
by induction that A c (Z d ) < l/(2d+ 2). However, this bound is less sharp than the 
bound (j8.53j) obtained in Example 8.3. 



is 



See e.g. jllO] for a brief discussion of transfer matrices for the hard-core lattice gas. 
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8.4.2 Lower bounds on A c 

We now turn to proving lower bounds on A c . Corollary 15. HI gives 

Ac(Z } " (A + l)A + i - (2d ~ 2id ' (8 ' 57) 

while Corollary 15.71 gives the slightly better bound 

x /rr2s (A — 1) A_1 (2d -If 4 - 1 1 , . 

A ^ ) * A^ ^ ~ 2* ■ (8 ' 58) 

For d = 2 the latter bound yields A C (Z 2 ) > 27/256. Either of these bounds shows, 
when combined with Example 8.3, that A c (Z d ) = Q(l/d) as d —>■ oo. 

To get quantitatively better bounds, we can take a supertree of the pruned SAW- 
tree of Z d and calculate exactly for it. In particular, if we choose a supertree that is 
eventually periodic, the computation reduces to finding the fixed point of a recursion. 

Example 8.6. The pruned SAW-tree of Z d is obviously contained in the 2d- 
branching tree, for which the computation has been performed in Example El6. This 
yields the bound (|8.57|) . [This argument obviously works for any infinite graph of 
maximum degree A, not just r L d \ 

Example 8.7. A smaller supertree of the pruned SAW-tree is obtained by giving 
the root 2d children and every other vertex 2d — 1 children. Applying Algorithm T, 
the recursion (until we reach children of the root) is w' \— > w/(l + w') . We end up 
with the bound ()8.58|) . [The argument again works for any infinite graph of maximum 
degree A.] 

Example 8.8. Another possibility (not the only one) is to take a large finite 
piece T of the pruned SAW-tree of Z 2 , and then repeat it periodically (attaching to 
each leaf of T a copy of T starting at the root, and repeating this). It is easily seen 
that this gives a supertree of the pruned SAW-tree. We therefore use the recursion 
w e ff, leaves | — ► ^eff.root an d demand that there exist an attractive fixed point in (—1,0]. 
We have various choices about how to order the edges to define the pruning. For 
instance, we can choose an ordering of the edges at each vertex in a translation- 
invariant way, find the first k levels of the pruned SAW-tree, and then repeat these 
periodically. Alternatively, we can order the vertices arbitrarily at the root of the 
SAW-tree, and at all other vertices of the SAW-tree order the edges according to the 
angle the path turns through (e.g. 0, +tt/2, — vr/2) in taking the step. This defines a 
pruned SAW-tree, and we can take the first k levels and repeat periodically as before. 
We conjecture that either method gives bounds converging to A C (Z 2 ) as k grows, but 
we do not know how to prove this. 

Remark. In Examples 8.2 and 8.7, we have proven for Z 2 the rigorous bounds 

27 4 
0.105468... = — < A C (Z 2 ) < — = 0.16. (8.59) 
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By extensions of those arguments and with a little calculation, these bounds can be 
narrowed further; it would be interesting to see how far one can go. It is worth noting 
that Todo [110J has given the extraordinarily precise numerical estimate 

A C (Z 2 ) = 0.119 338 88188(1), (8.60) 

obtained by using transfer matrices and the phenomenological-renormalization method 
(a variant of finite-size scaling). Furthermore, his computations up to L = 38 show 
[TTT] that 16 

X C (S 38 ) = 0.119 365(1) , (8.61) 

which by Example 8.4 provides an upper bound on A C (Z 2 ). See also Guttmann [S3] for 
an earlier and only slightly less precise estimate of A C (Z 2 ), obtained by series analysis. 
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16 What Todo jlllj actually computed is the number A£ for which 

(a) the transfer matrix has a unique eigenvalue of largest modulus for w G (— A£,0], and 

(b) the transfer matrix has two dominant eigenvalues of largest modulus for w = — A£. 

Let us show that A£ = X c (Sl). Let be the cylinder of width L and length n. By monotonicity, 
X C (S^) decreases in n and hence has a limit as n — > oo; it is easy to see that this limit is X c (Sl)- 
Now, the Beraha-Kahane- Weiss theorem ^3 ll(J3j tells us that 

(a') for every e > 0, there exist 5 > and no < oo such that ZWn) has no (real or complex) zeros 

within a distance 5 from the interval [— A£ + e, 0] when n > no; and 
(b') there exist (possibly complex) zeros of Z^ n ) tending to — A£ as n — > oo. 

Since — X C (S^) is the closest zero to the origin of Z^ n ), it follows easily from (a') and (b') that 
A* - X C (S L ). 
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